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Abstract 


Timber  supply  analysis  techniques  for  broad  geographical  levels 
are  summarized  by  land  allocation,  timber  growth  and  yield,  short- 
term  harvest  flows,  and  long-term  timber  investment  modeling  compo- 
nents. Representative  techniques  of  these  major  analytical  compo- 
nents are  summarized  in  tabular  form  for  supply  regions.  Aspects  of 
aggregation  and  uncertainty  in  timber  supply  analyses  are  also 
discussed. 
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Aggregate  Timber  Supply  Analysis 

Ralph  J.  Alig,  Bernard  J.  Lewis,  and  Paul  A.  Morris 


INTRODUCTION 

Aggregate  timber  supply  analyses  examine  long-run 
trends  in  stumpage  supply  and  identify  opportunities  for 
influencing  the  course  of  those  trends  through  manage- 
ment of  the  forest  resource.  These  analyses  support  the 
formulation  and  direction  of  public  and  private  timber 
policies.  In  this  report,  aggregate  timber  supply  analysis 
refers  to  techniques  for  analyzing  timber  resource  sup- 
plies for  geographic  areas  ranging  from  substate  to 
national  levels. 

This  paper  presents  a  broad  overview  of  major 
analytical  techniques  for  aggregate  timber  supply 
analysis,  and  is  intended  as  a  reference  document  for 
analysts  of  natural  resource  supplies.  It  provides  brief 
descriptions  of  major  techniques  and  their  applications, 
and  includes  critiques  from  those  who  have  designed, 
used,  or  examined  these  techniques.  Bibliographical 
references  refer  the  reader  to  specific  documents  con- 
taining detailed  operational  information  (i.e.,  computa- 
tional structure,  operating  costs,  etc.)  for  the  cited 
analytical  techniques. 

Analyses  of  timber  supply  have  been  approached 
from  a  variety  of  perspectives;  the  objectives  of  a  par- 
ticular timber  supply  study  and  data  availability  dictate 
the  specific  analytical  techniques  employed.  As  a  start- 
ing point,  one  might  consider  the  physical  quantity  or 
"supply"  of  timber  in  a  particular  geographic  region.  At 
a  given  point  in  time  a  certain  amount  of  wood  fiber  is 
present,  and  the  growth  and  mortality  within  this  timber 
base  may  be  modeled  and  projected  under  different 
management  assumptions  and  for  a  variety  of  future 
time  periods. 

The  "physical  supply"  of  timber  is  governed  by  both 
physical  and  biological  factors  that  affect  timber  growth 
and  mortality,  while  the  broader  concept  of  economic 
timber  supply  is  based  on  economic  feasibility  of  grow- 
ing and  harvesting  timber.  Economic  timber  supply  is  a 
schedule  of  the  quantities  of  stumpage  at  different 
prices  that  producers  are  willing  to  offer  for  sale  in  a 
specific  market  area  in  a  given  time  period  under  the  ex- 
isting institutional  framework.2  This  definition  of  timber 
supply  is  adopted  for  use  in  this  report.  It  is  not 
synonymous  with  that  for  "physical  supply,"  because  it 
incorporates  the  notion  of  responses  to  different  timber 
prices  based  on  opportunity  costs  related  to  the  manage- 
ment and  availability  of  standing  timber.  As  an  example 
of  this  distinction,  degrees  of  timber  accessibility  are 
more  fully  reflected  by  economic  supply  than  "physical 
supply"  estimates. 

An  additional  influence  on  timber  supply  is  the  form 
of  institutional  and/or  administrative  constraints  which 

intended  meanings  of  certain  terms  are  defined  in  the  Glossary. 


reduce  the  amount  of  timber  from  that  which  would 
otherwise  be  available  from  an  economic  standpoint. 
Parks  and  other  areas  withdrawn  from  the  economically 
available  timber  base  are  a  familiar  example  of  institu- 
tional restrictions  influencing  timber  supply.  From  an 
administrative  perspective,  constraints  on  timber  pro- 
duction may  take  the  form  of  targets  for  other  forest 
outputs— wildlife,  water,  etc.— which  must  be  achieved 
concurrently  with  the  production  of  timber,  thus  reduc- 
ing the  economic  timber  supply. 

The  combination  of  all  of  the  above  factors  com- 
plicates the  analysis  and  description  of  timber  supply. 
Furthermore,  from  the  aggregate  standpoint,  often  these 
distinctions  may  be  lost  in  the  actual  process  of  data 
aggregation.  If  a  given  model  fails  to  incorporate  one  or 
more  of  these  distinctions,  the  validity  and  subsequent 
utility  of  its  representation  of  actual  timber  supply 
becomes  tenuous.  However,  the  necessity  for  accurate 
information  depicting  timber  supplies  at  all  levels  of 
geographic  resolution  remains.  Marty  (1969)  sum- 
marizes the  general  utility  of  modeling  timber  supply 
from  an  aggregate  perspective: 

An  aggregate  timber  supply  model,  using  a 
classification  of  timber  supply  sources  based  on 
physical  and  economic  characteristics  and  an  intra- 
source  projection  technique  which  explicitly  takes 
into  account  the  more  ultimate  determinants  of  tree 
growth  and  removal,  could  be  useful  in  predicting 
the  supply  response  to  various  changes  in  economic 
conditions.  The  model  would  allow  the  analyst  to 
estimate  timber  inventories  under  various  assump- 
tions about  future  prices  and  production  costs.  This 
information  would  provide  an  estimate  of  the 
change  in  total  volume  of  stumpage  that  would 
become  available  under  different  future  stumpage 
price  expectations,  thus  providing  a  traditional 
long-run  supply  curve  for  timber.  Similarly,  it  would 
be  possible  to  estimate  the  effect  on  supply  of 
changes  in  production  cost  factors  like  timber  grow- 
ing technology,  labor  costs,  taxes,  and  so  forth. 

SCOPE 

This  report  describes  the  major  analytical  techniques 
utilized  in  aggregate  timber  supply  analyses  that  are 
available  in  published  or  mimeographed  form;  it  ex- 
cludes numerous  techniques  developed  for  proprietary 
use  in  private  industry  and  other  methods  currently 
being  developed.  Techniques  examined  are  almost  ex- 
clusively those  developed  for  use  in  the  United  States  or 
Canada,  which  are  generally  representative  of  available 
techniques.  Fries  et  al.  (1978)  provide  examples  of  long- 
term  timber  yield  modeling  in  other  countries. 
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Analytical  techniques  for  describing  aggregate  de- 
rived demand  for  timber,  based  on  society's  demand  for 
products  derived  in  part  from  timber,  are  not  consid- 
ered explicitly  in  this  report.  Therefore,  the  topics  of 
price  formation  and  projection  are  also  not  discussed.3 

Because  physical,  economic,  and  institutional  factors 
act  simultaneously  to  influence  aggregate  supply,  their 
representation  in  current  analytical  techniques  must  be 
considered.  However,  a  specific  breakdown  of  models 
along  these  lines  is  impractical,  because  many  modeling 
efforts  incorporate  each  of  these  three  factors  to  vary- 
ing degrees  within  their  analytical  frameworks.  The 
following  approach  for  describing  analytical  techniques 
is  adopted  in  this  report. 

Aggregate  timber  supply  analysis  requires  incorpora- 
tion of  the  following  processes: 

1.  Projecting  changes  in  the  forest  land  base. 

2.  Sampling  and  measuring  the  forest  inventory  ac- 
cording to  various  criteria  and  assembling  inven- 
tory information  to  facilitate  monitoring  of  changes 
in  various  inventory  strata. 

3.  Projecting  the  growth  and  yield  of  the  timber 
resource,  including  its  response  to  varying  degrees 
of  management. 

4.  Projecting  timber  removals. 

5.  Projecting  timber  investment  potentials  and 
strategies  over  time. 

This  report  describes  existing  analytical  methodology 
in  terms  of  four  of  the  above  areas.  The  collection  and 
organization  of  timber  inventory  data  (item  2  above)  is 
not  discussed  in  depth  but  is  considered  whenever  it  is 
relevant  to  the  other  areas  noted  above.  Accurate  inven- 
tory information  is  a  critical  input  to  the  design  and 
effectiveness  of  timber  supply  analyses.  While  the 
distinction  between  assembling  inventory  data  and 
subsequent  analysis  of  the  data  is  not  always  clear,  in 
this  report  inventory  data  will  be  regarded  as  the  infor- 
mation base  which  underlies  supply  analysis  and  allows 
it  to  proceed.  Detailed  discussions  regarding  the  nature 
of  inventory  data  and  other  timber  supply  data  may  be 
found  in  Clawson  and  Stewart  (1965),  Davis  (1976),  Tee- 
guarden  (1977),  Adams  et  al.  (1979),  Evans  (1979),  Eisen- 
man  et  al.  (1980),  and  USDA  Forest  Service  (1981,  1982). 

Institutional  constraints  in  supply  modeling  are  con- 
sidered as  they  impact  the  forest  land  base  and  the 
establishment  of  output  objectives  (e.g.,  timber  and 
water)  from  forest  lands.  However,  techniques  for 
analyzing  supply  in  the  presence  of  particular  institu- 
tional harvest  requirements  (e.g.,  nondeclining  even- 
flow)  are  only  briefly  addressed.  Chappelle  et  al.  (1976), 
Bell  (1977),  Johnson  and  Scheurman  (1977),  and  Hann 
and  Brodie  (1980)  examine  some  of  the  major  techniques 
applied  to  analyze  cases  where  these  constraints  exist. 

^Aggregate  "demand"  studies  have  historically  estimated  con- 
sumption quantities  rather  than  demand  functions,  analogous  to 
production  estimates  versus  supply  functions.  See  Gregory  et  al. 
(1971)  for  a  review  of  "demand"  projection  techniques  and  Adams 
and  Haynes  (1980)  for  an  example  of  combined  supply  and  demand 
modeling  for  price  projection. 


ORGANIZATION 

The  report  first  presents  a  historical  overview  of 
salient  techniques  that  have  been  employed  in  analyzing 
aggregate  timber  supply.  Current  analytical  techniques 
are  then  described  according  to  four  major  component 
areas  which  contribute  to  aggregate  timber  supply 
analysis— land  allocation,  growth  and  yield  changes  of 
the  timber  inventory,  harvest  flows,  and  timber  invest- 
ment strategies.  It  is  recognized  that  these  components 
are  not  distinctly  separate  in  theory  or  in  their  use  in 
supply  modeling.  Such  a  classification  does,  however, 
permit  the  separation  of  areas  in  which  specific  tech- 
niques have  tended  to  develop.  This  in  turn  facilitates 
the  description  of  relative  availability  of  methodologies 
and  permits  monitoring  of  progress  of  research  in  each 
of  these  four  component  areas.  For  each  area,  descrip- 
tions of  specific  models  are  provided  as  examples  of  cur- 
rent methodologies. 

Following  the  discussion  of  modeling  components,  cer- 
tain ramifications  of  data  aggregation  and  uncertainty 
with  respect  to  the  reliability  of  timber  supply  projec- 
tions are  briefly  explored,  reflecting  the  fact  that  these 
aspects  often  seem  to  be  inadequately  addressed  in 
aggregate  timber  supply  studies.  The  report  concludes 
with  a  general  synthesis  of  relative  strengths  and  weak- 
nesses of  existing  timber  supply  methodology. 


HISTORICAL  REVIEW 

Long-range  studies  of  aggregate  timber  supplies  in  the 
United  States  originated  with  Hough's  (1878,  1880) 
"Report  upon  Forestry"  for  the  U.S.  Department  of 
Agriculture.  Major  studies  that  followed  Hough's  work 
are  listed  in  Appendix  A,  which  includes  the  date, 
primary  analytical  techniques  employed,  and  scope  of 
each  study.  Successive  studies  have  in  general  been 
more  comprehensive  and  reliable  because  of  improve- 
ments in  analytical  techniques  and  data  availability 
(USDA  Forest  Service  1982). 

Quantitative  determination  of  "growth  and  drain" 
guided  major  timber  supply  assessments  (e.g.,  USDA 
Forest  Service  1948)  for  most  of  the  last  one-hundred 
years.  Difficulties  inherent  in  rigid  application  of  the 
"drain  ratio"  concept  (ratio  of  current  removals 
(harvest)  to  current  growth)  have  been  discussed  by 
Vaux  and  Zivnuska  (1952)  and  Bentley  and  Davis  (1976),. 
among  others.  For  example,  the  drain  ratio  indicates 
changes  in  total  volume  of  timber  growing  stock,  but 
does  not  provide  a  measure  of  growth  capability. 

Until  the  latter  half  of  this  century,  projections  of  ag- 
gregate timber  supplies  in  published  studies  typically 
focused  on  physical  measures  such  as  growth,  mortality, 
and  removals,  with  little  economic  interpretation.  The 
McSweeney-McNary  Act  of  1928  directed  the  Secretary 
of  Agriculture  to  undertake  comprehensive  surveys  of 
timber  supplies  to  judge  their  capabilities  in  meeting 
timber  "requirements"  (Bentley  and  Davis  1976).  Thus, 
USDA  Forest  Service  studies  under  the  McSweeney- 
McNary  authority  (e.g.,  USDA  Forest  Service  1948) 
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focused  on  a  "balance  of  the  timber  budget  of  the  United 
States,"  with  maintenance  of  certain  relative  price 
trends  of  forest  products. 

One  of  the  earliest  systematic  methodologies  used  for 
the  projection  of  physical  measures  of  aggregate  timber 
supply  was  the  stand  projection  method  (e.g.,  USDA 
Forest  Service  1965).  Availability  of  computer  tech- 
nology and  an  improving  data  base  in  the  form  of  rela- 
tively consistent  timber  inventory  data  facilitated  the 
development  and  application  of  the  Timber  Resource 
Analysis  System  (TRAS)  in  the  1960's  based  on  the  stand 
table  approach  (Larson  and  Goforth  1970,  1974).  Since 
then  several  other  timber  inventory  projection  tech- 
niques (e.g.,  individual  tree  based  systems)  have  become 
candidates  for  use  in  regional  timber  studies;  the  major 
ones  will  be  examined  later  in  the  section  concerning 
timber  inventory  projection  methodologies. 

Economic  measures  of  alternative  timber  supply 
levels  were  first  used  in  the  periodic  USDA  Forest  Serv- 
ice appraisals  in  the  1970  Outlook  (USDA  Forest  Service 
1973).  In  this  analysis,  supply  projections  were  based  on 
balancing  net  growth  and  removals  in  the  year  2000; 
utilizing  this  criterion,  economically  available  supplies 
of  softwood  sawtimber  were  then  estimated.  The 
estimates  of  economic  availability  depended  to  a  major 
degree  on  judgment  rather  than  upon  analytically  de- 
rived relationships.  Little  attention  was  given  to  the 
timber  management  activities  or  total  output  level  im- 
plied in  a  sustainable  "net  growth-removals"  balance.4 

The  notion  of  supply  and  demand  as  schedules  or 
functions  (of  quantities  versus  prices)  was  more 
thoroughly  introduced  into  the  decennial  Forest  Service 
assessments  of  timber  supply  in  the  1980  Resources 
Planning  Act  (RPA)  Assessment  (USDA  Forest  Service 
1982).  Statistical  analysis  of  the  response  of  private 
stumpage  supply  to  changes  in  prices  and  timber  inven- 
tory levels  was  a  basis  for  short-term  harvest  projec- 
tions (Adams  and  Haynes  1980)  in  a  comprehensive, 
interregional,  economic  equilibrium  approach.  In  a  dif- 
ferent type  of  supply  and  demand  analysis,  the  Forest 
Policy  Project  of  the  Pacific  Northwest  Regional  Com- 
mission estimated  private  timber  supply  by  maximizing 
present  net  benefits  in  conjunction  with  a  downward- 
sloping  curve  for  regional  timber  demand  (Rahm  1981). 

The  gradual  trend  in  aggregate  timber  supply  analy- 
sis toward  greater  concern  for  economic  variables  was 
preceded  by  several  seminal  papers  in  the  latter  half  of 
this  century.  Vaux  and  Zivnuska  (1952)  stressed  the  use 
of  supply  and  demand  concepts  from  neoclassical 
microeconomic  theory,  thus  focusing  not  on  single  quan- 
tities, but  on  price-associated  functions.  Gregory  (1955) 
undertook  further  theoretical  investigations  concerning 
economic  timber  growth  goals,  and  Vaux's  (1954)  con- 
struction of  cost-efficient  supply  schedules  for  Califor- 
nia sugar  pine  represented  some  of  the  earliest  em- 
pirical work  in  this  area  (Bentley  and  Davis  1976). 

Timber  supply  schedules  based  upon  the  theoretical 
cost-efficiency  concept  have  been  used  in  numerous 

'Problem  analysis:  Methods  for  analyzing  national  timber  supply. 
By  Thomas  J.  Mills,  on  file  at  Rocky  Mountain  Forest  and  Range  Ex- 
periment Station,  Fort  Collins,  Colo.  37  p.  1976. 


analyses  since  Vaux's  1954  study  (Vaux  1973,  Mont- 
gomery et  al.  1975,  Clawson  and  Hyde  1976,  Jackson  and 
McQuillan  1979).  Use  of  statistically  based  modeling  of 
aggregate  timber  supply  trends  has  not  been  as 
widespread,  in  part  because  of  problems  with  data 
availability.  McKillop  (1967),  Adams  (1974),  and  Robin- 
son (1974)  provide  some  of  the  earliest  documented 
investigations  of  statistically  based  timber  supply  sched- 
ules. Use  of  statistical  supply  functions  in  a  national 
timber  assessment  effort  is  documented  by  Adams  and 
Haynes  (1980),  who  also  discuss  integration  of  timber 
supply  techniques  within  an  overall  supply  and  demand 
system  that  projects  stumpage  prices. 

Duerr  (1977)  and  Beuter  (1979),  among  others,56  cri- 
tique past  national  timber  reviews  and  outline  major 
problems  to  be  resolved  in  future  timber  planning 
studies.  Major  criticisms  of  recent  domestic  supply 
projection  methodologies  have  centered  on  problems  in 
adequately  addressing  (1)  growth  and  yield  impacts  of 
shifting  timber  management  strategies  over  time,  (2)  re- 
sponses of  nonindustrial,  private  timber  supplies,  and 

(3)  aggregation  with  respect  to  factors  that  might  affect 
future  timber  availability. 

In  summary,  aggregate  timber  supply  analysis  has 
progressed  from  studies  guided  by  application  of  the 
ratio  of  timber  drain  (removals)  to  growth,  to  com- 
prehensive supply  and  demand  modeling  that  involves 
projection  of  stumpage  prices.  Before  the  1970's,  pro- 
jections of  aggregate  timber  supplies  typically  consisted 
only  of  physical  measures.  The  relatively  recent  model- 
ing of  timber  supply  and  demand  equilibrium  has  in- 
cluded statistically  based  stumpage  supply  schedules, 
with  feedback  mechanisms  linking  stumpage  supply  and 
timber  inventory.  Major  criticisms  of  recent  aggregate 
timber  supply  studies  have  been  directed  at  the  model- 
ing of  the  dynamics  of  timber  investment  opportunities 
and  investment  behavior  on  nonindustrial,  private 
lands,  and  also  at  aggregation  schemes  that  impede  dis- 
aggregation of  results,  provide  uneven  representation, 
and  are  difficult  to  use  as  a  basis  for  policy  formulation. 

COMPONENTS  OF  TIMBER  SUPPLY  ANALYSIS 

Depending  upon  the  objectives  and  degree  of  preci- 
sion sought  in  a  particular  study,  varying  degrees  of 
abstraction  are  involved  in  analyzing  the  complex  inter- 
relationships of  socioeconomic,  physical,  and  biological 
elements  of  timber  supply.  For  expository  purposes,  four 
components  of  aggregate  timber  supply  modeling  are 
discussed  in  this  report:  (1)  land  allocation,  (2)  pro- 
gression of  the  timber  inventory,  (3)  harvest  flows,  and 

(4)  long-term  investment  or  management  strategies. 

"Critique  assessment  topic  no.  1.  Forecasting  and  RPA:  1980 
RPA.  By  Dennis  L.  Schweitzer,  Con  Schallau,  A.  Schuler,  Robert 
Stone,  and  Gregory  Super,  on  file  at  USDA  Forest  Service  Policy 
Analysis  Staff  Office,  Washington,  D.C.  1980. 

*A  critical  analysis  of  timber  supply  research  objectives  for  re- 
gional and  national  assessments.  By  John  H.  Beuter,  on  file  at  the 
Rocky  Mountain  Forest  and  Range  Experiment  Station,  Fort 
Collins,  Colo.  43  p.  1978. 
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Interrelationships  among  the  four  components  are 
quite  important  given  that,  in  theory,  their  dynamic 
interactions  necessitate  a  simultaneous  solution  in 
modeling.  The  dynamic  context  of  a  time  path  of  inter- 
connected components  implies  that  modeling  the  pro- 
gression of  aggregate  timber  inventories  cannot  be 
effectively  achieved  without  considering  the  time  path 
of  the  decision  variables  for  land  use  allocation  and 
management.  Furthermore,  the  composition  of  future 
timber  inventories  and  harvests  is  strongly  related  to 
present  management. 

While  segregation  of  timber  supply  analytical  tools  by 
these  categories  allows  for  a  more  straightforward 
discussion  (corresponding  to  their  typically  separate 
treatment  in  the  literature),  such  a  segregation  should 
not  obscure  the  interdependencies  among  these  compo- 
nents. Furthermore,  considerations  of  the  linkages  of  the 
four  components  to  the  overall  economy  are  important  in 
analyzing  timber  supply  patterns,  although  such  rela- 
tionships have  seldom  been  examined  in  depth  (Adams 
and  Blackwell  1973).  Many  of  the  more  useful  tech- 
niques and  associated  assumptions  for  analyzing  these 
four  timber  supply  components  will  now  be  examined. 


1.  LAND  ALLOCATION 

The  quantity,  quality,  and  location  of  land  that  will  be 
allocated  to  timber  production  have  generally  been 
estimated  exogenously  in  aggregate  timber  supply 
studies,  as  have  the  related  ownership  patterns  and 
objectives  of  management.  For  example,  exogenous  pro- 
jections of  commercial  forest  area  have  been  used  in 
estimating  stumpage  supply  in  the  Timber  Assessment 
Market  Model  (TAMM);  such  projections  are  based  on 
observed  trends  and  subjective  judgment  (Adams  and 
Haynes  1980).  While  Adams  and  Haynes  recognize  that 
the  acreage  devoted  to  timber  production  or  allowed  to 
remain  in  a  forested  state  is  not  independent  of  the 
timber  price  and  growth  projections  made  in  TAMM, 
they  conclude  that  the  decisions  involving  the  many 
possible  uses  of  land  are  too  complex  for  direct  inclusion 
as  endogenous  processes  in  TAMM. 

Some  exogenous  projections  of  land  use,  as  in  the 
TAMM  approach,  treat  land  allocation  to  timber  grow- 
ing as  a  "residual"  use  that  is  calculated  by  sifting  out 
other  uses  of  perceived  higher  economic  value  (e.g.,  agri- 
culture), along  with  institutional  uses  (e.g.,  wilderness), 
from  the  total  potential  base  for  commercial  forest 
acreage  (Wall  1981).  This  approach  does  not  allow  for 
modeling  the  dynamic  interaction  between  changes  in 
the  forest  land  base  and  relative  timber  stumpage 
prices. 


MODELING  APPROACHES 

Burnham  (1973)  used  a  Markov  land  use  simulation 
model  to  depict  intertemporal  land  use  shifts  in 
estimating  future  United  States  cropland  availability. 
Forest  land  was  one  of  six  land  uses  projected,  based 


upon  land  uses  in  the  Southern  Mississippi  alluvial 
valley  between  1950  and  1969.  The  finite  Markov  chain 
process  involved  estimating  probabilities  of  land  use 
shifts  based  on  land  movements  between  use  groups 
over  a  historical  time  period,  and  then  using  these  prob- 
abilities to  project  future  land  uses  at  alternative  points 
in  time.  Projecting  the  future  implications  of  past  land 
use  trends  in  this  manner  does  mask  the  causative 
variables  and,  therefore,  does  not  allow  quantification 
of  the  processes  underlying  land  use  shifts. 

In  a  study  concerned  with  loss  of  wildlife  habitat, 
MacDonald  et  al.  (1979)  also  applied  a  Markov  chain 
model  to  project  bottomland  hardwood  acreages  in  the 
Lower  Mississippi  alluvial  plain  to  1995.  Their  prob- 
abilistic model  was  based  upon  historical  land  use  data 
obtained  primarily  from  interpretation  of  aerial 
photographs  and  other  forms  of  remotely  sensed  im- 
agery. Economic  analysis  in  the  study  showed  that 
agricultural  crop  production  was  more  economically 
attractive  than  investments  in  natural  bottomland  hard- 
wood stands  or  plantations. 

Economic  analysis  of  the  allocation  of  land  as  an  input 
or  factor  of  production  to  timber  growing  has  typically 
used  the  economic  rent  concept,  with  an  unchanging 
alternative  rate  of  return  assumed  for  a  particular 
forest  ownership  (e.g.,  Vaux  1973,  Hyde  1980).  Hyde 
(1980)  applied  the  concept  of  land  rent  (i.e.,  net  annual 
economic  return  per  acre  from  a  land  use)  to  demon- 
strate the  use  of  a  variable  land  "production  function." 
This  function  indicates  land  used  for  timber  production 
at  various  prices  and  the  associated  efficient  harvest 
levels.  Hyde's  land  production  function  is  based  upon 
changing  silvicultural  technologies  for  each  land  class 
brought  into  production.  As  the  various  land  classes 
receive  different  degrees  of  management  intensity,  dif- 
ferent efficient  yields  per  acre  result.  Access  costs,  de- 
pendent on  distance  to  a  site  and  topography,  influence 
the  degree  of  timber  management  applied  to  different 
land  classes.  This  approach  assumes  that  the  landowner 
has  a  profit  maximization  goal,  and  that  higher  stump- 
age prices  will  attract  less  economically  productive 
acres  into  timber  production. 

Access  costs  in  timber  supply  determination  are  also 
addressed  by  Ledyard  and  Moses  (1976)  in  a  conceptual 
model  of  land  use  that  combines  transportation  cost  con- 
siderations (e.g.,  Bradley  1972)  and  capital  theory.  Their 
model  of  land  use  in  long-run,  steady-state,  competitive 
equilibrium  is  based  upon  economic  principles  de- 
scribed by  Von  Thunen  almost  a  century  and  a  half  ago 
(Von  Thunen  1966).  Von  Thunen's  model  allocated  land 
to  alternative  crops  based  upon  yielding  maximum  eco- 
nomic land  rents,  given  crop  prices  and  transportation 
costs  to  market.  Ledyard  and  Moses  extend  Von 
Thunen's  basic  framework  to  model  dynamic  aspects  in- 
volving interactions  of  time  and  transportation  costs, 
which  have  substantial  implications  for  the  use  of 
capital  in  a  forestry  setting.  Their  focus,  like  Hyde's 
(1980),  is  on  a  long-run,  steady-state  equilibrium  rather 
than  on  identifying  the  path  to  equilibrium. 

No  forestry-based  models  currently  exist  with  the 
capability  for  modeling  interregional  linkages  in  land 
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allocation  patterns  involving  forestry.  A  system  for 
analyzing  large-scale,  interregional,  agricultural  land 
allocation,  the  Center  for  Agricultural  and  Rural  Devel- 
opment (CARD)  model,  considers  possible  conversions  of 
forest  to  cropland  when  estimating  interregional  shifts 
in  agricultural  production  based  on  least-cost  agricul- 
tural expansion  to  meet  exogenous  agricultural  com- 
modity requirements  (Huang  et  al.  1981).  However,  the 
CARD  model  is  presently  not  designed  to  estimate  the 
possible  reversion  of  cropland  to  forest. 

GASPLY.— Georgia  Supply  (GASPLY)  is  an  economic 
model  designed  to  analyze  forestry  investments  in 
Georgia  (Montgomery  et  al.  1975,  1976;  Robinson  et  al. 
1978).  The  major  theoretical  bases  of  the  GASPLY 
model— economic  optimization  procedures  for  land  allo- 
cation models— are  similar  to  those  employed  by  Vaux 
(1973)  and  Hyde  (1980).  GASPLY  estimates  long-run 
equilibrium  stumpage  price  and  quantity  and  the  associ- 
ated total  investment  cost.  Also  estimated  at  this 
equilibrium  point  are  the  number  of  acres  and  owners 
involved,  type  of  management  (e.g.,  plantation)  and  loca- 
tion by  planning  district,  and  the  site  quality  of  the 
timber  acreage.  The  criterion  for  GASPLY's  economic 
optimization  is  maximizing  present  net  worth  of  per- 
petual timber  rotations. 

The  economic  modeling  of  land  use  by  GASPLY  is 
based  on  the  relative  profitability  of  agricultural  and 
timber  production  uses,  assigning  to  a  given  land  parcel 
the  use  with  the  highest  economic  rent  (i.e.,  difference 
between  all  discounted  relevant  costs  and  revenues 
associated  with  production  schemes  into  perpetuity). 
Limits  are  set  on  the  amount  of  land  that  can  shift  be- 
tween timber  production  and  agriculture.  Only  the  most 
productive  timberland  found  in  farmer  ownerships  is 
considered  to  be  feasible  for  conversion  to  agriculture, 
while  the  amount  of  farmland  available  for  conversion 
to  timberland  is  restricted  to  idle  farmland. 

Idle  land  eligible  for  conversion  to  forest  is  assigned 
an  opportunity  cost  (when  calculating  soil  expectation 
value)  equal  to  the  economic  returns  that  could  be 
earned  annually  from  alternative  agricultural  use  of  the 
land.  Positive  soil  expectation  values  for  timber  produc- 
tion strategies  cause  such  idle  farmland  to  be  converted 
to  timberland  in  the  GASPLY  model.  The  agricultural  op- 
portunity cost  is  net  of  charges  for  clearing  presently 
forested  land,  and  the  model  allows  for  the  possibility 
that  idle  farmland  can  revert  naturally  to  timberland 
(i.e.,  custodial  management).  Similarly,  land  currently  in 
timber  production  is  shifted  to  agricultural  use  if 
associated  soil  expectation  values  are  negative  when 
alternative  agricultural  returns  are  considered. 

Montgomery  et  al.  (1975)  recommend  that  commercial 
forest  acreage  of  nonindustrial,  private  landowners  be 
reduced  in  GASPLY  input  because  of  economically  in- 
operable tracts  (e.g.,  land  around  homesites).  GASPLY 
adjusts  land  allocation  patterns  for  expected  population 
growth,  with  timberland  shifting  to  urban-related  land 
uses.  The  model  can  also  reflect  additional  ramifications 
of  urbanization,  such  as  property  tax  changes.  Impacted 
acres  that  are  held  in  speculation  (and  will  command 
high  prices  in  the  real  estate  market)  are  assigned 


higher  ad  valorem  property  taxes  that  may  render 
timber  production  uneconomical. 

GASPLY  can  adjust  forest  types  within  the  timberland 
base  to  reflect  the  possibility  that  the  timberland  will 
receive  custodial  care  in  lieu  of  plantation  or  intensive 
natural  stand  management  (Montgomery  et  al.  1975). 
Custodial  management  implies  no  active  timber  manage- 
ment other  than  continued  harvesting  and  protection 
against  fire,  insects,  and  disease.  For  acreage  assumed 
to  receive  custodial  management,  GASPLY  applies  an 
adjustment  to  the  pine  and  oak-pine  types  for  hardwood 
invasion. 

In  summary,  potential  timberland  acreage  subdivi- 
sions are  evaluated  by  GASPLY  under  three  timber 
management  intensities — pine  plantation,  intensively 
managed  natural  stand,  or  custodial— and  the  manage- 
ment intensity  producing  the  highest  present  worth  is 
chosen  for  that  acreage  subdivision  or  cell  if  it  can  ex- 
ceed an  accumulated  agricultural  rent.  The  selected 
timber  management  levels  and  associated  costs  and 
yields  are  assembled  into  a  long-run  timber  supply 
curve,  which  indicates  the  various  amounts  of  timber 
that  can  be  produced  at  different  average  unit  costs. 

ADDITIONAL  MODELING  CONSIDERATIONS 

The  long-run  supply  curve  is  based  upon  the  criterion 
of  economic  efficiency  on  the  part  of  the  timberland 
owner — i.e.,  maximizing  economic  returns  through 
timber  production.  It  also  embodies  major  assumptions 
that  owners  have  (1)  perfect  knowledge  regarding  future 
timber  stumpage  price  levels,  yields,  costs,  etc.; 
(2)  perfectly  competitive  stumpage  and  input  (factor) 
supply  markets;  and  (3)  no  entry  barriers  for  timberland 
investments  (other  than  cost).  Furthermore,  no  dynamic- 
adjustment  paths  (and  associated  costs  and  price 
changes)  and  spatial  detail  are  generally  considered. 
These  assumptions  underlie  most  economic  maximiza- 
tion models  of  timber  supply  (e.g.,  Hoyer  1975;  Jackson 
and  McQuillan  1979,  1980). 

Land  allocation  based  on  neoclassical  economic  effi- 
ciency principles  guided  Davis  (1972),  Dideriksen  et  al. 
(1977),  Hidlebaugh  (1980),  and  Shulstad  and  May  (1980) 
in  estimating  potentials  for  the  conversion  of  forest  land 
to  agricultural  uses.  The  fact  that  many  of  these  land 
use  conversion  or  investment  opportunities  have  not  al- 
ready been  undertaken  suggests  that  the  behavior  of  the 
diverse  class  of  private  landowners  is  not  fully  explain- 
able with  current  economic  efficiency  models  and  data. 
A  host  of  other  studies  of  the  behavior  of  nonindustrial, 
private  forest  landowners  support  this  proposition.7 

Modelers  in  general  have  had  limited  success  in  pro- 
jecting patterns  of  land  use  (Voelker  1975);  this  is  par- 
ticularly true  in  the  case  of  nonindustrial,  private  land 
management.  The  forestry  literature  is  replete  with 
observations  on  a  related  phenomenon  often  termed  the 
"small  forest  owner  problem."  Numerous  proposals  for 

7Problem  analysis:  Improving  the  timber  productivity  of  non- 
industrial  private  forest  lands.  By  J.  E.  deSteiguer,  on  file  at  North- 
eastern Forest  Experiment  Station,  Forestry  Science  Laboratory, 
Princeton,  W.  Va.  37  p.  1980. 
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attaining  more  "efficient"  timber  production  patterns 
on  these  lands  have  been  offered;  some  view  this  class  of 
forest  land  owner  as  the  key  to  significant  increases  in 
future  national  timber  supplies  (USDA  Forest  Service 
1982). 

The  factors  behind  the  apparent  deviation  of  private 
landowner  behavior  from  normative  economic  pathways 
have  not  been  unequivocally  identified,  but  may  be  the 
result  of  several  combined  characteristics  of  the 
market,  owners,  and  analyses  themselves.  Imperfect 
market  information,  uncertainty,  "noneconomic"  goals, 
and  lack  of  technical  skills  (Holley  1979)  are  examples  of 
possible  confounding  market  and  owner-related  factors. 
An  important  analytical  shortcoming  may  be  the  present 
inability  to  account  fully  for  all  the  relevant  returns, 
especially  nonmonetary  benefits,  and  costs  accruing  to 
a  landowner  from  alternative  land  uses. 

This  latter  problem  may  arise  primarily  from  data 
deficiencies  and/or  measurement  errors;  however,  the 
analytical  approach  and  theoretical  foundations  may 
also  be  partially  at  fault.  In  a  study  of  New  England 
landowner  behavior,  Binkley  (1979,  1981)  applied  a 
"utility  maximization"  approach  which  accounted  for 
the  influence  of  some  nontimber  outputs  in  timberland 
management  decisions.  Boyce  (1977,  1978)  has  also 
developed  techniques  for  describing  timber  manage- 
ment activities  in  the  context  of  multiple  benefits.  While 
these  approaches  have  not  been  extensively  tested, 
limitations  to  wider  application  in  the  short  term  appear 
to  center  on  data  deficiencies.  The  lack  of  accessible 
joint  production  data  at  aggregate  levels,  combined  with 
inherent  aggregation  problems  resulting  from  resource 
systems  which  have  both  varying  ecological  and  eco- 
nomic bases,  has  greatly  hampered  aggregate  multi- 
resource  production  analysis.  Problems  of  even  greater 
severity  on  the  demand  side  further  render  integrated 
forest  resource  analysis  a  very  difficult  task  at  ag- 
gregate levels. 

Deviations  of  timber  suppliers  from  norms  prescribed 
in  the  formulation  of  optimization  models  based  on 
financial  maximization  prompted  Marty  (1969)  to  recom- 
mend a  modeling  technique  adapted  to  multigoal  or 
mixed  goal  analysis.  His  recognition  of  the  difficulty  of 
this  task  is  still  relevant  today.  While  some  progress  has 
been  made  over  the  last  decade,  modeling  problems 
related  to  incommensurable  objectives  held  by  many 
timberland  owners  strongly  influence  all  four  dimen- 
sions of  timber  supply  analysis  referred  to  previously. 
This  contributes  substantially  to  uncertainty  in  timber 
supply  modeling  (see  section,  Aggregation  and  Uncer- 
tainty Considerations). 

Projection  of  land  use  allocation  based  on  historical 
relationships  among  economic  and  other  important 
variables  is  one  alternative  tool  to  complement  economic 
efficiency  models.  However,  few  statistically  based 
models  of  forest  acreage  trends  have  been  employed, 
primarily  because  of  substantial  data  problems.  Ade- 
quate time  series  data  for  forest  land  use  on  which  to 
base  statistical  functions  are  difficult  to  construct 
because  of  irregular  measurements  or  estimates  of  com- 
mercial forest  land  in  the  past.  Timberland  acreages  for 


a  particular  state  are  estimated  approximately  every 
ten  years  by  USDA  Forest  Service  Forest  Inventory  and 
Analysis  surveys  (e.g.,  Knight  1973).  Survey  cycles  differ 
among  states.  Attempts  to  augment  periodic  Forest  Serv- 
ice forest  acreage  data  with  other  major  data  sources  is 
hindered  by  differing  criteria  used  in  past  classifica- 
tions of  forest  lands. 

White  and  Fleming  (1980)  relied  upon  farm  woodland 
data  estimated  every  five  years  in  the  Census  of 
Agriculture  (e.g.,  USDC  Bureau  of  the  Census  1977)  in 
analyzing  factors  underlying  competing  land  use  shifts 
in  Georgia.  A  system  of  equations  was  statistically 
estimated  to  explain  simultaneous  changes  in  crop, 
pasture,  and  forest  acreages.  The  most  significant 
determinants  of  forest  acreage  were  the  amount  of  land 
in  farms  and  farm  crop  acreages  for  previous  years 
(Fleming  1980).  Regression  results  were  used  in  a  simu- 
lation model  to  investigate  the  impact  of  government 
land  diversion  programs  on  land  use  patterns  from  1980 
to  2000. 

In  summary,  the  fundamental  question  of  the  alloca- 
tion of  land  as  an  input  to  timber  production  has  been 
analyzed  using  several  distinct  approaches  in  aggregate 
timber  supply  studies.  National  timber  assessments 
have  generally  used  exogenous  estimates  of  timberland 
acreage,  based  on  computations  of  forest  land  as  a 
residual  use  (e.g.,  USDA  Forest  Service  1982).  The 
dynamics  involved  in  the  joint  determination  of  the 
timberland  base  and  relative  timber  stumpage  prices 
have  not  been  modeled,  nor  have  similar  interactions 
with  other  sectors,  particularly  agriculture.  In  short, 
forecasting  the  availability  and  accessibility  of  forest 
acreage  is  a  major  problem  in  aggregate  timber  supply 
analyses. 


2.  TIMBER  INVENTORY  PROJECTION 

A  number  of  methods  have  been  developed  to  project 
future  levels  of  timber  inventory  (growth,  mortality,  ac- 
cumulated stock,  and  structure)  as  it  would  occur  in  an 
unmanaged  forest,  as  well  as  in  response  to  certain 
management  practices  (Avery  and  Burkhart  1983).  Such 
methods  of  projecting  forest  development  can  be  divided 
into  two  broad  classes— direct  and  indirect— based 
upon  whether  the  method  is  applied  locally  to  the  same 
stand  from  which  the  data  used  for  projection  purposes 
are  obtained  (direct  methods),  or  whether  sample  data 
for  numerous  stands  are  synthesized  and  then  extrapo- 
lated to  other  stand  conditions  (indirect  methods). 

There  are  many  models  designed  to  project  timber 
inventory  development,  some  of  which  are  outlined  in 
appendix  B.  These  models  have  been  developed  for 
numerous  species  and  regions;  however,  models  for  a 
single  species  may  vary  considerably,  for  example,  in 
terms  of  the  sample  of  stand  conditions  on  which  the 
model  is  based  or  the  analytical  methods  employed. 
Evaluation  of  these  models  requires  the  synthesis  of 
many  quantitative  and  qualitative  factors  (Buchman  and 
Shifley  1983).  Surveys  of  inventory  projection  models 
are  also  provided  by  Williston  (1975),  Solomon  (1977a), 
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Farrar  (1979a),  Dudek  and  Ek  (1980),  Burkhart  et  al. 
(1981),  Moeur  and  Ek  (1981),  Trimble  and  Shriner  (1981), 
and  Hann  and  Riitters  (1982). 


DIRECT  METHODS 

Direct  projection  methods  are  those  in  which  the 
results  of  operations  performed  upon  a  data  sample  are 
applied  directly  back  to  the  population  from  which  the 
sample  was  obtained,  and  only  to  that  specific  popula- 
tion. The  most  familiar  of  these  methods  is  that  of  stand 
table  projection,  in  which  a  sample  of  growth  rates  is  ob- 
tained from  a  given  stand  and  then  applied  directly  back 
to  that  stand  for  projection  purposes. 

Stand  tables  are  composed  of  frequency  data  showing 
number  of  trees  according  to  such  classification  systems 
as  species,  diameter  at  breast  height  (d.b.h.),  or  height 
classes  (Husch  et  al.  1972).  Stand  tables,  commonly  ex- 
pressed on  a  per  acre  basis,  are  useful  in  depicting  the 
stand  structure  or  distribution  of  tree  sizes  and  species 
in  a  stand.  Stand  table  projection  models  use  estimates 
of  future  diameter  growth,  removals,  mortality,  and 
ingrowth  to  adjust  the  stand  table  in  an  accounting 
fashion  over  time.  Husch  et  al.  (1972)  further  discuss  the 
basis  and  the  advantages  and  disadvantages  of  the  gen- 
eral stand  table  projection  method. 

Although  stand  table  projection  models  may  be  used 
to  address  localized,  individual  stand  development  ques- 
tions, they  have  also  been  used  for  projecting  aggregate 
timber  supplies  when  forest-level  or  aggregate  stand 
questions  are  being  investigated.  To  answer  these 
broader  forest-level  questions,  the  assumption  is  made 
that  a  large,  regional  inventory  can  be  treated  as  one 
all-aged  stand. 

These  generalized  stand  table  projection  models  can 
use  conventional  timber  inventory  data  and  may  be 
relatively  simple  and  inexpensive  to  use.  Some  of  their 
disadvantages  are  that  they  need  a  large  data  base, 
pose  difficulties  when  different  forest  management 
schemes  are  being  modeled,  and  are  conceptually  dif- 
ficult to  interpret.  Another  disadvantage  is  that, 
because  the  forest  is  modeled  in  aggregate,  reliable  tree 
size  and  distribution  information  for  smaller  areas  is 
hard  to  extract.  The  following  example  of  a  stand  table 
projection  model  expands  on  some  of  these  points. 

TRAS.— The  Timber  Resource  Analysis  System 
(TRAS)  developed  by  Larson  and  Goforth  (1970)  has 
been  employed  as  the  inventory  projection  system  for 
many  recent  USDA  Forest  Service  timber  assessment 
studies  at  the  regional  and  national  level  (e.g.,  USDA 
Forest  Service  1982).  TRAS  has  continued  to  evolve 
since  its  inception  in  the  1950's,  and  an  expanded  ver- 
sion, TRAS-1980,  is  described  by  Alig  et  al.  (1982).  TRAS 
is  a  stand  table  projection  model  specifically  designed  to 
solve  three  types  of  problems:  to  reconcile  differences 
among  forest  surveys,  to  update  forest  surveys  com- 
pleted at  different  times  to  a  common  compilation  date, 
and  to  project  long-term  timber  inventory  changes. 

The  TRAS  system  includes  both  an  exponential  size 
class  distribution  model,  the  Q  method,  and  a  parabolic 
size  class  distribution  function,  the  non-linear  interpola- 


tion method  or  NLI  method  (Larson  and  Goforth  1974).  In 
relatively  large  aggregations  of  individual  even-aged 
and  uneven-aged  forest  stands,  the  overall  stand  table  of 
number  of  trees  by  diameter  class  will  approach  an  ex- 
ponential distribution,  for  which  the  Q  method  is  ap- 
propriate. The  NLI  method  was  designed  to  model  the 
parabolic-type  size  class  distribution  of  even-aged 
stands.  Because  collections  of  stands  are  modeled  in  ag- 
gregate, timber  supply  studies  and  such  aggregates  tend 
toward  the  exponential  distribution  (or  the  inverse 
J-shape  as  shown  in  fig.  1);  thus,  only  the  Q  method  or  ex- 
ponential distribution  model  is  discussed  in  this  report. 

The  Q  method  of  TRAS  is  based  on  Meyer's  (1952) 
observation  that  a  graphical  representation  of  the 
diameter  distribution  in  any  large  forest  area,  with  a 
mixture  of  stand  sizes  and  ages,  tends  toward  an  inverse 
J-shaped  curve  (fig.  1).  This  curve  is  expressed  as  an 
exponential  function  that  relates  the  number  of  trees 
per  acre  to  their  respective  d.b.h.  class.  Thus,  an 
average  annual  increase  in  diameter  is  related  to  a 
corresponding  potential  increase  in  the  number  of  trees. 
This  distribution  seems  to  be  generally  true  for  the  ag- 
gregation of  even-aged,  uneven-aged,  or  mixed  stand 
distributions. 

The  aggregation  capabilities  of  TRAS's  inventory  pro- 
jection system  depend  largely  on  the  assumption  that  ag- 
gregate "stands"  constructed  from  a  pooling  of  data 
from  much  smaller  geographic  subdivisions  are  an  ade- 
quate representation  of  the  combined  smaller  units,  as 
gaged  by  the  objectives  and  accuracy  requirements  of 
the  particular  timber  supply  study.  In  practice,  TRAS  is 
often  applied  to  aggregates  that  represent  millions  of 
acres.  Empirical  verification  of  adequate  representation 
by  such  aggregates  is  quite  elusive  because  of  the 
unavailability  of  independent  and  appropriate  data 
bases  upon  which  to  check  regional  abstractions. 
Larson's  test  using  South  Carolina  data  did  indicate  that 
as  few  as  ten  Renewable  Resource  Evaluation  sample 
plots  were  needed  to  contruct  a  "sound"  stand  table  for 
TRAS  projections.4 

The  potential  increase  in  inventory,  due  to  growth  of  a 
stand,  is  computed  as  the  annual  per  acre  change  in  the 
number  of  trees  in  each  2-inch  d.b.h.  class  (fig.  1).  The 
overall  change  in  inventory  is  obtained  by  subtracting 
the  annual  removals  and  mortality  from  the  just- 
computed  potential  increase  in  the  number  of  trees.  This 
resulting  number  of  trees  added  to  the  inventory  is  con- 
verted to  a  total  volume  by  using  the  average  volume  per 
tree  by  d.b.h.  class  and  the  area  of  commercial  forest 
land  (Larson  and  Goforth  1974). 

To  prevent  the  development  of  unrealistically  high  in- 
ventory basal  areas  during  projections,  a  set  of  con- 
straint equations  has  been  developed.  As  described  by 
Larson  and  Goforth  (1970),  this  set  of  equations  reduces 
the  radial  growth  and  ingrowth  and  increases  the  mor- 
tality as  the  total  softwood  and  hardwood  basal  area 
increases.  The  overall  timber  inventory  in  a  region  has 
generally  been  projected  as  two  separate  softwood  and 
hardwood  "average"  stands  by  the  TRAS  model  in 
aggregate  timber  supply  analyses,  although  the  projec- 
tions of  both  stands  are  influenced  by  the  total  stand 
basal  area  constraint. 
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The  mechanics  of  the  TRAS  projection  system  are 
relatively  simple,  and  the  accuracy  of  TRAS  inventory 
projections  essentially  depends  on  the  ability  to  ade- 
quately project  radial  growth,  removals,  and  mortality, 
by  diameter  class,  and  by  ingrowth.8  TRAS  was  designed 
to  project  timber  management  levels  that  are  implicit  in 
the  radial  growth,  mortality,  removals,  and  ingrowth 
data  that  comprise  the  bulk  of  the  input  to  the  stand 
table  projection  process  (i.e.,  TRAS  contains  no  growth, 
mortality,  or  removal  rate  coefficients  internally).  Thus, 
the  implicit  assumption  behind  the  standard  TRAS 
inventory  projection  is  that  timber  management  will  be 
implemented  at  basically  the  same  level  over  the  projec- 
tion period  as  that  which  formed  the  historically  based 
stand  table  input.  This  assumption  regarding  projected 
TRAS  parameters  has  been  strongly  questioned  because 
of  apparent  trends  in  intensifying  timber  management  to 
a  degree  different  from  that  implied  by  standard  TRAS 
projections.8 

TRAS  procedures  were  modified  in  the  1980  RPA 
Assessment  to  incorporate  impacts  of  such  managment 
shifts  (Barber  1980).  The  management  shifts,  based  on 
identified  timber  investment  opportunities,  were  trans- 
lated into  growth  increments  via  Barber's  modification 
of  radial  growth  and  mortality  rates  in  TRAS's  stand 
table.  The  modification  procedure  is  of  an  ad  hoc  nature 
and  was  devised  to  fit  an  immediate  need  of  the  1980 
RPA  Assessment  modeling.  The  cumbersomeness  and 
difficulty  in  monitoring  involved  in  TRAS's  modeling  of 

'Problem  analysis:  Timber  supply  analyses  for  RPA  Assess- 
ments. By  Darius  M.  Adams,  on  file  at  the  Rocky  Mountain  Forest 
and  Range  Experiment  Station,  Fort  Collins,  Colo.  29  p.  1980. 
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Figure  1.— Inverse  J-shaped  form  of  diameter  distribution  assumed 
in  TRAS  model  for  large  forest  areas  with  mixture  of  stand  sizes 
and  ages,  with  diameter  distribution  displacement  shown  for  one 
year's  growth  (Larson  and  Goforth  1974). 


shifts  in  timber  management  intensity  is  one  of  the  major 
reasons  that  the  USDA  Forest  Service  plans  to  supplant 
it,  at  least  in  some  regions  where  management  shifts  are 
significant,  for  the  next  major  national  timber  supply 
assessment.8  This  was  prompted  by  the  need  to  examine 
policy  questions  that  frequently  involve  changes  in 
timber  management  practices  or  intensity. 

The  use  of  TRAS  to  evaluate  certain  changes  from  the 
initial  state  over  long  projection  periods  has  also  been 
questioned.9  For  example,  in  a  shift  from  an  inventory  of 
predominantly  old  growth  to  largely  young  growth,  such 
a  transition  could  cause  the  lower  d.b.h.  classes  to  shift 
from  predominantly  suppressed  to  predominantly  domi- 
nant crown  classes.  The  corresponding  acceleration  in 
growth  over  time  reflected  by  the  TRAS  model  may  not 
be  an  adequate  representation  of  this  change.  The  oc- 
currence of  such  deviations  from  steady-state  systems 
has  led  to  suggestions  that  TRAS  be  limited  to  short-run 
inventory  updates  in  some  parts  of  the  country. 

Other  reported  limitations  in  the  application  of  TRAS 
in  some  timber  supply  analysis  tasks  have  seemingly  in- 
creased with  the  marked  evolution  of  timber  supply 
modeling  and  associated  needs  in  recent  times.  In  part, 
some  applications  of  TRAS  extend  beyond  those  origi- 
nally intended  in  its  design.  For  example,  Row  (1977) 
points  out  severe  limitations  of  TRAS  with  regard  to  the 
previously  mentioned  insensitivity  to  trends  in  manage- 
ment levels  and  investment  in  timber  growing  activities, 
as  well  as  environmental  and  multiple-use  constraints. 
While  noting  that  the  TRAS  system  has  been  validated 
reasonably  well  in  decade-by-decade  biological  projec- 
tion of  timber  inventories  for  which  it  was  designed, 
Row  asserts  that  several  major  timber  supply  investiga- 
tions (Adams  1975,  Holley  et  al.  1975)  have  employed 
TRAS  as  a  timber  inventory  projection  module  primarily 
because  there  was  no  other  model  available  with  ade- 
quate supporting  data. 

With  respect  to  these  criticisms  of  TRAS,  the  use  of 
"average  stands"  to  analyze  broad  regional  areas  has 
been  regarded  as  acceptable  when  detailed  input 
parameters  for  smaller  areas  are  lacking  or  when  there 
is  no  need  to  know  how  small  parts  of  the  total  resource 
are  changing  with  time  (Oswald  1978).  The  basic 
assumption  that  a  large  regional  inventory  can  be 
treated  as  one  large,  all-aged  stand  may  be  reasonable 
in  some  parts  of  the  country.  The  assumption  that  a 
large,  all-aged  stand  represents  a  regional  timber  inven- 
tory structure  is  questioned,  however,  when  even-aged 
management  is  widely  practiced  and  clearcutting  is  the 
accepted  harvest  method. 

Two  methods  have  been  proposed  to  minimize  the 
above  weakness.  One  proposal  is  to  divide  the  inventory 
into  discrete  strata  and  to  use  the  TRAS  model  to  project 
each  stratum  separately  (Oswald  1978).  The  second  pro- 
posal is  to  divide  the  inventory  into  discrete  strata  and 
to  use  the  TRAS  stand  table  projection  technique  for  the 
older  strata  and  to  use  a  yield  table  projection  technique 
for  stands  regenerated  under  more  intensive  even-aged 
management. 

"Personal  communication  with  A.  R.  Stage,  USDA  Forest  Service, 
Moscow,  Idaho. 
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Alternative  approaches  to  the  TRAS  system  for  use  in 
regional  timber  inventory  projections  include  the 
Timber  Resource  Inventory  Model  (TRIM)  and  the  South- 
ern Pine  Age-Class  Timber  Simulator  (SPATS),  models 
that  currently  are  being  developed  and  tested.  The 
TRIM  model  employs  a  more  disaggregated  approach 
than  TRAS,  as  it  stratifies  timberland  acreage  by  age 
classes,  site,  species,  management  intensity,  and  stock- 
ing level  (Tedder  1983).  Acres  can  be  shifted  into  dif- 
ferent levels  of  timber  management  in  two  ways  (1)  pre- 
specified  management  shifts  over  the  entire  planning 
horizon  or  (2)  management  shifts  based  on  economic  in- 
vestment behavior  estimated  internally  by  the  model. 
Yield  table  projections  for  ten-year  age  classes  for  these 
different  levels  of  timber  management  are  made  using 
an  approach  to  normality  technique. 

The  SPATS  5-year  age-class  model  is  intermediate  be- 
tween the  TRAS  and  TRIM  approaches  in  regard  to  level 
of  aggregation.  It  can  model  three  different  forest  types 
separately  for  the  Southeast  and  Midsouth:  (1)  planted 
pine,  (2)  natural  pine,  and  (3)  oak-pine  (Brooks  1983a). 
One  average  site  class  and  stocking  class  are  used  and 
an  approach  to  normality  is  also  utilized  to  adjust  for 
stocking  variations.  Reforestation  practices  can  be 
directly  simulated  in  the  SPATS  model,  including 
estimating  annual  planting  over  a  period  by  using  equa- 
tions developed  using  historical  data  (e.g.,  cost-share 
subsidies). 

INDIRECT  METHODS 

Indirect  timber  inventory  projection  methods  are 
based  on  a  range  of  growth  and/or  yield  observation 
samples  that  are  utilized  in  estimating  the  parameters  of 
projection  models.  The  resultant  models  are  then  ap- 
plied to  subject  stands  larger  in  nature  than  those  from 
which  the  original  growth/yield  data  were  obtained.  In 
essence,  therefore,  these  methods  rely  on  a  two-phase 
process:  the  first  provides  or  describes  the  data  to 
which  the  projection  model  is  to  be  applied  (i.e.,  the  in- 
ventory stands  or  plots);  the  second  sample  provides  the 
data  for  model  construction,  and  for  that  reason  it  is 
more  detailed  than  the  former.  The  model  is  then  ap- 
plied back  to  the  first  sample  for  drawing  inferences 
about  the  entire  forest  studied.  Indirect  models  may  be 
described  in  terms  of  two  broad  classes: 

A.  Whole  stand  yield  tables 

1.  Normal 

2.  Empirical 

3.  Variable  density 

B.  Variable  density  growth  and  yield  models 

1.  Whole  stand 

2.  Size  class  and/or  growth  component 

3.  Individual-tree 

a.  Distance-independent 

b.  Distance-dependent 

Whole  Stand  Yield  Tables 

Yield  tables  describe  stand  development  and  yields  by 
age,  site  quality,  and,  on  occasion,  other  descriptors. 


The  most  familiar  examples  of  such  tables  are  those 
described  as  normal  yield  tables,  which  attempt  to 
describe  the  aggregate  or  whole  stand  yield  of  "fully 
stocked,"  undisturbed  natural  stands  by  site  index  and 
stand  age.  Such  tables  are  constructed  from  measure- 
ments of  many  stands,  but  the  stand  characteristics  are 
observed  only  at  a  single  point  in  time;  thus,  they  provide 
estimates  of  net  production  only  (Curtis  1972).  The  pat- 
tern obtained  by  plotting  yield  or  stand  characteristics 
over  age  is  assumed  to  represent  a  real  growth  pattern 
for  stands  of  this  density.  While  they  have  the  advan- 
tage of  being  relatively  easy  to  construct,  this  is  gener- 
ally outweighed  by  the  fact  that  most  normal  yield  tables 
are  only  applicable  to  relatively  pure  and  even-aged 
stands,  and  they  provide  information  only  on  aggregate 
stand  variables.  Such  tables  do  not  adapt  well  to  man- 
agement purposes,  and  they  are  also  difficult  to  use  for 
"non-normal"  stands. 

A  variation  on  normal  yield  tables  are  empirical  yield 
tables,  which  are  constructed  on  the  basis  of  average 
stand  stocking  conditions,  as  opposed  to  fully  stocked 
stands.  Such  tables  are  constructed  from  sample  plots 
with  varying  degrees  of  stocking,  with  multiple  regres- 
sion techniques  generally  utilized  to  determine  average 
stand  characteristics.  While  these  tables  allow  for 
somewhat  greater  flexibility  in  projecting  stand  yields, 
they  still  do  not  overcome  most  of  the  difficulties  stated 
above  with  respect  to  normal  yield  tables.  Another  im- 
portant disadvantage  of  most  such  yield  tables  is  that 
they  provide  relatively  little  information  on  stand 
growth.  While  yields  are  furnished  for  the  "terminal" 
stand  condition  (i.e.,  end  of  the  projection  period),  little 
information  is  provided  on  initial  stand  conditions  and 
growth  patterns  that  resulted  in  these  terminal  yields. 
This  greatly  reduces  the  utility  of  most  yield  tables  for 
comparing  management  alternatives. 

Variable  density  yield  tables  are  the  result  of  at- 
tempts to  resolve  these  problems  concerning  growth 
representation  in  timber  inventory  projection  modeling. 
Such  tables  provide  yields  for  a  range  of  stand  densities 
for  each  age  and  site  combination.  However,  early 
variable  density  yield  tables  were  based  on  only  single- 
point-in-time  observations,  and  did  not  effectively  pro- 
vide further  information  on  growth  patterns.  More 
recently  developed  variable  density  growth  and  yield 
modeling  techniques  reflect  the  incorporation  of  growth 
data  from  permanent  plot  records  and  corresponding 
equations  to  characterize  forest  growth  as  well  as  its  in- 
tegral forest  yield.  Variable  density  growth  and  yield 
tables  have  been  geared  primarily  to  silvicultural  alter- 
natives in  forest  management  (Ek  and  Dudek  1980). 


Variable  Density  Growth  and  Yield  Models 

Variable  density  growth  and  yield  models  are  con- 
structed to  provide  both  growth  and  yields  correspond- 
ing to  different  stand  densities.  They  generally  include 
specific  growth  functions,  as  well  as  yield  functions 
derived  from  the  integration  and/or  summation  of 
growth  functions.  Such  models  may  be  utilized  for  the 
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construction  of  variable  density  yield  tables,  and  thus 
they  are  of  great  value  in  contributing  to  the  ability  to 
represent  potential  results  of  different  silvicultural 
management  alternatives. 

Whole  Stand  Models 

These  techniques  have  generally  utilized  differential 
equations  to  describe  the  rates  of  change  of  various 
components  of  stand  development.  An  important  initial 
study  in  this  area  was  that  of  Buckman  (1962),  in  which 
the  growth  rate  of  even-aged  stands  was  expressed  as  a 
function  of  age,  site,  and  stocking  in  a  differential 
expression.  With  this  approach,  yields  are  obtained  by 
integrating  the  growth  rate  equations  over  time. 
Depending  upon  the  specific  form  of  the  rate  model,  dif- 
ferentiation may  be  employed  to  derive  stocking  levels 
which  maximize  growth. 

Another  significant  study  that  exemplifies  important 
elements  of  this  approach  was  that  conducted  by  Clutter 
(1963),  in  which  growth  rate  models  were  developed  and 
subsequently  integrated  to  provide  closed-form  ex- 
pressions for  yield  predictions.  Data  for  this  analysis 
were  obtained  from  5-  and  10-year  remeasurements  of 
thinned  permanent  plots  of  loblolly  pine.  The  first  phase 
of  the  analysis  involved  the  development  of  regression 
equations  to  predict  per  acre  cubic  foot  and  basal  area 
yields  from  age,  site,  and  basal  area.  Regressions  for 
basal  area  growth  and  cubic  foot  growth  were  then  used 
to  obtain  prediction  equations  for  future  basal  area  and 
cubic  foot  volume  of  a  given  present  stand;  this  was 
achieved  by  integration  of  the  aforementioned  growth 
equations  with  respect  to  age.  Sullivan  and  Clutter 
(1972)  later  revised  Clutter's  (1963)  loblolly  pine  model 
and  addressed  two  statistical  problems:  (1)  dependent 
parameters  within  a  system  of  equations  and  (2)  obser- 
vations from  remeasured  permanent  plots  that  are  not 
independent. 


Size  Class  and/or  Growth  Component  Models 

Other  techniques  which  utilize  differential  equations 
to  model  forest  development  may  project  growth  com- 
ponents by  size  or  diameter  classes.  Moser's  (1974) 
system  of  nonlinear  differential  equations  for  predicting 
ingrowth,  survivor  growth,  and  mortality  by  diameter 
classes  exemplifies  this  modeling  approach  (see  also 
Moser  et  al.  1979).  An  additional  example  is  the  SHAF 
simulation  model  (Adams  et  al.  1974,  Ek  and  Monserud 
1979),  based  upon  a  stand  table  projection  system 
modified  to  incorporate  stand  growth  component  equa- 
tions given  by  Ek  (1974)  and  later  by  Adams  and  Ek 
(1974). 

Input  for  the  SHAF  model  includes  site  index  and 
number  of  trees  by  2-inch  diameter  classes.  Three  com- 
ponents of  stand  growth  by  diameter  classes  are  then 
estimated  using  nonlinear  difference  equation  models. 
These  components  include  ingrowth,  upgrowth  (i.e., 
number  of  trees  rising  from  one  diameter  class  to  the 
next  higher  class),  and  mortality.  Given  these  three 


growth  component  expressions,  the  number  of  trees  in 
each  diameter  class  at  the  end  of  a  given  growth  period 
may  be  obtained.  By  applying  conversion  factors  to  the 
number  of  trees  in  each  diameter  class,  basal  areas  and 
volumes  for  each  class  may  also  be  calculated.  It  should 
be  noted  that  variations  within  different  models  may  be 
manifest  in  different  forms  for  components  representing 
survivor  growth,  mortality,  and  regeneration. 

Individual-Tree  Simulation  Models 

A  common  characteristic  of  the  preceding  kinds  of 
variable  density  growth  and  yield  models  is  their  reli- 
ance upon  aggregate  stand  or  size  class  characteristics 
such  as  basal  area  or  number  of  trees  for  modeling  of 
forest  development.  The  most  recent  approach  for 
modeling  forest  growth  and  yield  consists  of  individual- 
tree  simulation  models,  in  which  stands  are  described 
on  the  basis  of  characteristics  of  individual  trees,  which 
then  may  be  combined  to  form  stands.  These  models  con- 
sider a  set  or  list  of  individual  trees  and  associated  tree 
characteristics  for  a  plot  and  simulate  the  growth  of 
each  tree  by  explicit  or  implicit  "potential"  growth  func- 
tions; these,  in  turn,  are  modified  by  expressions  reflect- 
ing tree  size  competition  (Ek  1977).  In  this  manner,  tree 
spatial  locations  and  individual  tree  characteristics 
influencing  tree  growth  and  resultant  yields  may  be 
explicitly  considered.  Individual-tree  models  may  be  di- 
vided into  two  general  classes:  (1)  those  that  are 
independent  of  intertree  distances,  and  (2)  those  that  re- 
quire intertree  distances  as  a  necessary  input  (Munro 
1974). 

Distance-Independent  Models. — Individual-tree 
distance-independent  models  do  not  require  intertree 
distances  but  classify  the  competitive  status  of  the  sub- 
ject tree  by  comparing  its  characteristics  (size,  crown 
ratio,  etc.)  in  relation  to  all  other  trees  in  a  sample  plot. 
Each  sample  tree  is  weighted  with  an  expansion  factor 
to  indicate  the  number  of  trees  per  acre  it  represents. 

The  advantages  of  this  type  of  model  are  that  inputs 
consist  primarily  of  conventional  inventory  data  (al- 
though some  modifications  are  required);  it  can  model 
any  age  structure  or  species  mixture;  and  it  provides 
relatively  detailed  information  on  tree  and  stand  param- 
eters and  on  the  effects  of  stand  management.  Disadvan- 
tages of  this  kind  of  model  include  its  relative  complexity 
and  the  inability  to  predict  the  growth  of  specific  single 
trees  with  any  reliability  (Munro  1974).  It  should  be 
noted  that  all  of  the  individual-tree  models  described 
herein  are  similar  in  that  they  involve  some  type  of  tree 
accounting,  product  conversion,  and  summary  capabil- 
ity, with  the  actual  stand  dynamics  treated  by  specific 
program  sections  or  subprograms  for  survivor  tree 
growth,  mortality,  and,  sometimes,  ingrowth  or  regener- 
ation. Somewhat  distinct  from  this  dynamic  part  of  the 
model,  there  may  be  subprograms  for  management  ac- 
tivities, such  as  harvesting  or  thinning.  Other  models 
have  identified  specific  options  for  insect  and  disease 
impact  studies  and  habitat  analyses. 

Individual-tree  distance-independent  models  have 
been  employed  to  simulate  growth  and  yield  for  in- 
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dividual  forest  stands,  small  forested  areas,  and  for  en- 
tire states.  Criticisms  are  sometimes  voiced,  however, 
concerning  the  aggregation  capabilities  of  these  models 
as  well  as  the  potential  costs  of  adapting  them  to  such 
large-scale  assessments.  The  extent  to  which  the 
number  of  plots  needed  to  represent  forest  conditions  of 
interest  would  expand  in  moving  from  substate  to 
regional  or  national  timber  assessments  requires  fur- 
ther examination,  as  does  the  attendant  cost  and  ac- 
curacy of  the  resultant  biological  modeling.  Conversely, 
an  individual-tree  model  capable  of  projections  for  20 
species  may  be  no  more  complex  (and  probably  less  so) 
to  construct  and  use  than  20  "simpler"  models  for  20 
different  species.9  The  difficulties  inherent  in  combining 
diverse  growth  and  yield  models  are  well  known  to 
resource  projection  analysts. 

Individual-tree  distance-independent  models  have 
been  developed  by  Lemmon  and  Schumacher  (1962), 
Botkin  et  al.  (1972),  Goulding  (1972),  Stage  (1973),  and 
USDA  Forest  Service  (1979),  among  others.  Two  ex- 
amples of  this  type  of  model  are  described  below. 

STEMS.— An  individual-tree  distance-independent 
model  considered  potentially  useful  for  regional  or  na- 
tional analyses  is  the  Stand  and  Tree  Evaluation  and 
Modeling  System  (STEMS),  formerly  known  as  the  Forest 
Resources  Evaluation  Program  (FREP)  (USDA  Forest 
Service  1979,  Hahn  et  al.  1979,  Moeur  and  Ek  1981). 
STEMS  is  a  forest  growth  model  designed  to  describe 
stand  development,  with  or  without  management  ac- 
tivities. The  model's  timber  growth  projection  system 
has  been  applied  in  the  Lake  States  (Smith  and  Raile 
1979,  Jakes  and  Smith  1980)  and  is  being  tested  in  the 
Pacific  Northwest;  plans  exist  to  calibrate  it  for  use  in 
the  Central  States  and  Southeast  as  well. 

STEMS  will  handle  a  variety  of  forest  conditions  in- 
cluding pure  or  mixed  stands.  Growth  and  yield  for 
even-  and  uneven-aged  stand  conditions  can  be  simu- 
lated if  the  relevant  biological  response  equations  have 
been  developed  and  tested  for  the  particular  geograph- 
ical area.  STEMS  can  be  used  to  (1)  update  or  project 
plot  tree  lists  to  a  future  time,  or  (2)  evaluate  various 
silvicultural  management  alternatives.  The  model  has 
numerous  timber  management  options  available. 

STEMS  uses  basic  inventory  data  or  data  from  other 
sources  such  as  stand  and  stock  tables  to  construct  an 
input  list  of  individual  trees  representative  of  the  forest 
condition  of  interest.  This  tree  list  contains  species, 
diameter  at  breast  height,  crown  ratio,  and  an  expan- 
sion factor  (per  acre)  for  each  tree.  If  crown  ratio  is 
missing  from  the  input  list,  it  is  approximated  by  using  a 
relationship  with  diameter  and  basal  area. 

The  model's  operational  portion,  which  uses  growth 
equations  for  individual  trees  to  transform  the  input  tree 
list  into  an  updated  output  tree  list,  has  three  parts: 
potential  growth  function,  modifier  function,  and  mor- 
tality function.  The  potential  growth  function  estimates 
how  rapidly  trees  would  grow  in  d.b.h.  if  they  did  not  in- 
teract with  other  stems.  The  modifier  function  reduces 
the  potential  change  to  correspond  with  observed  data 
from  permanent  growth  plots.  The  mortality  rule  com- 
putes a  probability  of  death  for  a  tree  based  on  its 
species  and  growth  rate. 


Output  from  the  model  is  a  list  of  trees  showing 
species,  new  d.b.h.,  and  new  status  if  the  tree  dies,  plus 
associated  product  conversions  and  summary  tables  of 
stand  characteristics.  The  projection  system  does  not 
currently  handle  events  such  as  ownership  changes. 

Prognosis.— In  addition  to  STEMS,  another  individual- 
tree  distance-independent  model  with  potential  for 
aggregate  timber  supply  analyses  is  the  Prognosis  model 
for  stand  development  (Stage  1973,  1979;  Stage  et  al. 
1980).  This  model  has  been  extensively  used  in  the 
northern  Rocky  Mountains.  It  has  been  adapted  to  all 
major  tree  species  found  in  Idaho  for  both  even-aged 
and  uneven-aged  stands  and  may  be  applied  to  similar 
coniferous  types  found  in  the  Pacific  Northwest. 

The  Prognosis  model  predicts  future  development  of 
stands  through  displays  of  long-range  values  of  tree 
diameters,  heights,  crown  ratios,  tree  species  composi- 
tion, and  understory  shrub  composition  and  coverage. 
Like  the  STEMS  model,  Prognosis  uses  species-specific, 
internal  growth  equations  to  simulate  the  growth  and 
yield  of  individual  trees,  which  can  then  be  extrapolated 
to  a  stand  basis.  As  with  most  individual-tree  models,  it 
is  not  restricted  to  specific  cover  types;  rather  it  is  ap- 
plicable to  stands  containing  any  mixture  of  species  or 
age  and  size  classes  that  grow  as  a  community  (Stage 
1973).  Additional  capabilities  include  conversion  of  tree 
dimensions  to  conventional  units  of  timber  yield  in  total 
cubic  and  board  feet,  as  well  as  information  related  to 
wildlife  habitat  changes,  watershed  protection,  and 
recreational  values  that  can  be  derived  from  the  funda- 
mental tree  and  shrub  variables. 

The  Prognosis  model  consists  of  a  set  of  computer  pro- 
grams for  combining  current  silvicultural  knowledge 
with  past  growth  data  to  make  a  forecast  of  stand  devel- 
opment reflecting  the  effects  of  thinning,  fertilization, 
and  regeneration  harvests.  This  program  is  capable  of 
making  projections  for  a  period  as  long  as  a  full  rotation. 

Input  to  the  Prognosis  model  consists  of  habitat  type, 
shrub  cover  and  height,  and  conifer  stocking.  Input  tree 
lists  are  described  by  species,  d.b.h.,  and  crown  ratio. 
Removal  of  individual  trees  from  an  inventory  plot  are 
used  to  reflect  harvests.  The  model  depicts  tree  growth 
through  a  collection  of  submodels  representing  growth 
of  trees  at  different  stages  in  their  life  cycle.  These 
stages  are  as  follows: 

1.  Regeneration  establishment  that  follows  seedlings 
through  the  first  10  years  of  growth. 

2.  Regeneration  development  that  represents  growth 
of  established  seedlings  and  saplings  up  to  the  time 
they  reach  5  inches  d.b.h. 

3.  Poles  and  larger  trees. 

The  regeneration  establishment  submodel  predicts 
the  stocking  characteristics  of  small  (1/300  acre)  plots. 
This  modeling  of  stocking  rate  requires  inputs  describ- 
ing habitat  type,  slope,  aspect,  elevation,  topographic 
position,  time  since  disturbance,  site  preparation 
method,  overstory  composition  and  density,  planting, 
and  distance  from  seed  source. 

The  regeneration  development  submodel  predicts 
height  and  diameter  increments  of  trees  less  than  5 
inches  d.b.h.  Height  increment  is  modeled  using  the 
height  and  species  of  the  subject  tree,  overstory  density, 
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conifer  and  shrub  competition,  habitat  type,  slope, 
aspect,  and  elevation.  The  diameter  increment  is  a  func- 
tion of  species  height  and  predicted  height  increment. 

The  growth  of  pole  and  larger  trees  is  represented  by 
a  set  of  species-specific  equations  that  predict  diameter 
increments  from  values  of  diameter,  crown  ratio,  crown 
competition  factor,  basal  area  in  trees  larger  than  the 
subject  tree,  slope,  aspect,  elevation,  habitat  type,  and 
geographic  region.  The  height  increments  of  these  trees 
are  dependent  on  predicted  diameter  increment,  height, 
diameter,  and  habitat  type. 

The  Prognosis  model  specifies  mortality  as  a  function 
of  species  d.b.h.  and  stand  density  (Hamilton  and  Ed- 
wards 1976).  As  with  the  growth  functions,  these  rela- 
tionships should  be  based  on  data  from  the  area  to 
which  the  model  is  to  be  applied. 

Output  from  the  Prognosis  model  is  separated  into 
individual  tree  summaries  and  forest  stand  summaries. 
Individual  tree  records  consist  of  species,  d.b.h.,  crown 
ratio,  past  diameter  growth,  and  trees  per  acre  associ- 
ated with  each  of  the  original  sample  tree  records. 
Forest  stand  records  consist  of  numbers  of  trees  per 
acre,  distribution  of  d.b.h.,  volume  characteristics  of  the 
total  stand  in  cubic  feet,  surface  areas  of  boles  in 
square  feet,  and  accumulated  tree  height.  In  addition, 
total  cubic  foot  volume  representing  the  growth  of  the  in- 
itial tree  population  and  the  total  cubic  foot  volume 
representing  mortality  are  summarized  (Stage  1973). 

The  Prognosis  model  has  the  capability  to  represent 
effects  of  pest  populations  (e.g.,  Douglas-fir  tussock 
moth)  on  growth  and  mortality  (Crookston  1978,  Monse- 
rud  1978,  Monserud  and  Crookston  1982).  It  may  also  be 
utilized  to  resolve  discrepancies  between  inventory 
growth  rates  and  model-based  estimates,  and  to  rescale 
model-based  growth  and  mortality  rates  to  represent 
nonmodeled  silvicultural  treatments  such  as  fertiliza- 
tion or  reduction  in  shrub  competition.  A  parallel-stand 
growth  processor  permits  interstand  effects  such  as 
harvest  allocation  or  pest  contagion  to  be  represented.9 

Distance-Dependent  Models. — The  second  type  of 
individual-tree  model  requires  intertree  distances  as  a 
necessary  input.  Some  representative  distance- 
dependent  models  have  been  developed  by  Newnham 
(1964),  Mitchell  (1969),  Arney  (1974),  Ek  and  Monserud 
(1974),  and  Daniels  and  Burkhart  (1975). 

In  these  models,  "individual  trees"  on  a  plot  are 
assigned  certain  initial  size  and  spatial  distributions. 
The  trees  are  then  "grown"  according  to  some  function 
of  their  size,  site,  competitive  status,  and  sometimes  a 
random  component  representing  microsite  and/or  gen- 
etic variability.  Competitive  status  for  each  tree  is  quan- 
tified in  terms  of  a  competition  index  that  is  a  function  of 
the  tree's  size  and  the  distance  to  and  size  of  its  neigh- 
bors. Mortality  is  regulated  as  a  function  of  competition 
index,  tree  size,  and/or  growth.  Yield  estimates  are 
made  by  applying  volume  equations  to  the  dimensions  of 
the  trees  (Curtis  1972),  as  in  other  individual-tree 
models. 

Individual-tree  distance-dependent  models  may  pro- 
vide more  detailed  stand  and  tree  measurements  than 
distance-independent  ones.  Several  disadvantages  of 
this  type  of  model  are  that  they  require  (1)  intertree 


distances  as  an  input  (this  is  usually  not  collected  as  in- 
ventory information),  (2)  a  meaningful  competition  index 
for  individual  trees,  and  (3)  more  computing  time. 

Individual-tree  distance-dependent  models  can  be 
considerably  more  sensitive  to  harvest  treatments  and 
reproduction  response  than  diameter  class  models  (Ek 
and  Monserud  1979).  However,  Ek  and  Monserud  also 
suggest  that  whole  stand  models  may  project  aggregate 
stand  characteristics  almost  as  well  as  individual-tree 
distance-dependent  models  at  a  much  lower  cost,  within 
the  range  of  data  used  to  construct  the  former.  Because 
of  their  detail  and  flexibility,  individual-tree  distance- 
dependent  models  may  be  most  useful  in  analyses  of  op- 
timal silvicultural  alternatives  (Adams  and  Ek  1974). 

Although  the  growth  and  yield  estimates  are  con- 
sidered to  be  fairly  accurate,  currently  the  greatest 
obstacle  to  the  use  of  these  models  for  state  and  regional 
analyses  are  the  data  requirements  for  intertree 
distance  coordinates  and  the  expense  of  the  required 
computer  time.  Another  disadvantage  is  their  use  of 
stochastic  methods  for  deducting  mortality.9  As  a  conse- 
quence, the  outcome  of  the  simulation  is  itself  a  random 
process  that  should  be  replicated  to  estimate  the  ex- 
pected value.  For  long-term  projections  or  for  high  mor- 
tality, the  variance  of  the  output  depends  on  the  initial 
numbers  of  tree-records  in  the  simulation. 

PTAEDA.— The  model  PTAEDA,  developed  by  Daniels 
and  Burkhart  (1975),  is  representative  of  the  single-tree 
distance-dependent  models.  PTAEDA  simulates  the 
growth  of  loblolly  pine  under  a  wide  range  of  manage- 
ment alternatives. 

The  two  major  subsystems  in  PTAEDA  deal  with  the 
generation  of  an  initial  precompetitive  stand  and  the 
growth  and  dynamics  of  that  stand.  After  PTAEDA  was 
developed  for  old-field  stands,  management  subroutines 
for  unmanaged  plantations  were  developed  to  simulate 
the  effects  of  site  preparation,  fertilization,  and  thinning 
(Daniels  and  Burkhart  1975).  Trees  are  "grown"  annu- 
ally according  to  their  size,  site  quality,  and  intertree 
competition.  The  probability  of  survival  for  each  tree  is 
calculated  from  a  function  relating  each  tree's  individ- 
ual vigor  and  its  competitive  stress  as  measured  by  esti- 
mates of  photosynthetic  potential.  This  probability  is 
used  to  stochastically  determine  annual  mortality. 

Effects  of  certain  management  practices  are  incor- 
porated into  PTAEDA  by  modifying  the  original  growth 
function.  The  efficiency  of  a  site-preparation  program  is 
expressed  as  the  degree  to  which  a  cutover  site  ap- 
proaches old-field  conditions.  Growth  reductions  on 
cutover  land  are  assumed  to  be  solely  due  to  competing 
vegetation,  because  degradation  in  site  quality  caused 
by  past  management  practices  could  be  described  by  in- 
itially specifying  a  lower  site  index  (Daniels  and 
Burkhart  1975). 

Similarly,  the  response  to  fertilizer  treatments  could 
be  described  as  an  increase  in  site  quality,  as  suggested 
by  Ek  and  Monserud  (1974)  and  Hegyi  (1974).  For  this 
reason,  a  site  adjustment  factor  that  acts  as  a  multiplier 
for  site  index  for  fertilized  stands  is  built  into  the  model. 

Intermediate  thinnings  result  in  a  decrease  in  com- 
petitive stress  for  the  remaining  trees.  These  trees 
become  more  capable  of  competing  for  the  limiting 
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resources  such  as  light,  water,  nutrients,  or  growing 
space.  The  response  to  these  thinnings  is  moderated  by 
a  tree's  own  potential  for  growth  as  measured  by  a  func- 
tion of  crown  size  (Daniels  and  Burkhart  1975). 

Estimates  of  basal  area  per  acre,  trees  per  acre,  total 
stem  cubic  foot  volume,  and  total  aboveground  biomass 
are  also  given  by  PTAEDA.  To  further  describe  stand 
conditions,  the  output  summary  includes  the  mean, 
standard  deviation,  and  range  of  relevant  tree  dimen- 
sions; the  stand  diameter  distribution;  the  average 
height  of  each  diameter  class  for  live  trees;  the  trees 
removed  in  thinning;  and  the  trees  lost  to  mortality. 
Farrar  (1979a)  indicates  predicted  yields  agree  with 
published  reports  for  thinned  and  unthinned  old-field 
plantations. 


SUMMARY 

Direct  inventory  projection  methods  have  been  used 
in  most  large-scale  timber  supply  analyses.  The  utility  of 
the  TRAS  model  (Larson  and  Goforth  1970,  1974)  has 
been  demonstrated  in  the  last  several  USDA  Forest 
Service  timber  assessments.  Although  TRAS  is  fairly  in- 
expensive to  run,  can  represent  any  geographical  level 
by  the  use  of  "average"  acres,  and  has  a  projection 
system  with  fairly  simple  mechanics,  it  has  been  criti- 
cized in  several  respects.  A  major  criticism  is  that  TRAS 
is  not  very  tractable  as  a  tool  for  simulating  changes  in 
timber  management  intensities. 

Among  indirect  methods  for  inventory  projection, 
individual-tree  distance-independent  models  have  often 
been  proposed  as  a  replacement  for  TRAS.  These 
models  have  the  capacity  to  efficiently  simulate  silvicul- 
tural  operations  as  well  as  to  provide  regional  growth 
and  yield  projections.  Models  of  this  type  that  have 
proven  most  useful  are  the  STEMS  model  (USDA  Forest 
Service  1979)  and  the  Prognosis  model  (Stage  1973). 

The  individual-tree  distance-dependent  models  are 
capable  of  providing  accurate  growth  and  yield  projec- 
tions as  well  as  detailed  information  about  stand  struc- 
ture (e.g.,  Daniels  and  Burkhart  1975).  Conversely,  the 
large  data  requirements  of  intertree  coordinates  and 
excessive  computer  time  have  thus  far  prohibited  use  of 
these  models  for  large-scale  analyses. 

Although  the  "inventory  projection"  models  discussed 
were  designed  for  somewhat  different  purposes,  various 
opportunities  exist  for  linking  them  in  regional  or 
national  timber  supply  studies.  For  example,  individual- 
tree-based  models,  such  as  STEMS,  may  be  more  pre- 
cisely classified  as  variable  density  growth  and  yield 
table  generators.  They  could  provide  "harvest  schedul- 
ing" models  with  the  necessary  timber  growth  inlorma- 
tion  to  help  drive  an  overall  timber  inventory  projection 
system.  Expanded  and  improved  regional  growth  and 
yield  modeling  by  such  techniques  will  depend  largely 
on  progress  in  assembling  growth  and  yield  data  in  a 
consistent  manner  for  all  regions,  which  would  facilitate 
modeling  of  alternative  timber  management  possibil- 
ities. In  particular,  it  is  important  that  growth  and  yield 
model  output  forms  be  similar. 


3.  HARVEST  FLOWS 

Harvest  flows  represent  the  timber  volumes  removed 
over  time  from  existing  timber  inventories  or  stock  sup- 
plies of  stumpage.  The  modeling  of  these  flows  is  the 
central  concern  of  aggregate  timber  supply  studies. 
Linkages  among  harvest  flows  and  the  two  components 
of  aggregate  timber  supply  analysis  discussed 
previously — land  allocation  and  timber  growth  and 
yield— involve  dynamic  interrelationships  between  bio- 
logical and  economic  forces.  Furthermore,  because 
timber  inventory  is  both  invested  capital  and  product, 
modeling  of  harvest  flows  is  also  conceptually  linked  to 
timber  investment  processes  that  impact  future  inven- 
tories from  which  harvests  will  be  drawn. 

This  intertemporal  linkage  of  harvest  decisions  and 
longer  term  investment  decisions  may  involve  adjust- 
ments in  timber  inventory  structure  and  volumes  to 
achieve  long-term  objectives,  such  as  profit  maximiza- 
tion, with  investments  (e.g.,  fertilization)  or  disinvest- 
ments (e.g.,  harvest)  based  strongly  on  future  expecta- 
tions of  the  timber  economy.  Timber  market  fluctuations 
in  the  short  term  may  prompt  previously  unforeseen 
revisions  of  long-term  adjustment  paths  so  that  short- 
term  gains  may  be  realized.  Navon  (1982)  observes  that 
the  time  path  of  aggregate  harvest  responds  to  market 
conditions  when  present  net  worth  is  maximized,  as 
would  be  the  case  with  industrial  and  some  nonindus- 
trial  private  ownerships.  For  public  agencies,  however, 
this  path  is  manipulated  to  provide  a  stable  supply  of 
timber  and/or  other  socially  dictated  objectives. 

Thus,  the  timing  and  intensity  of  harvests  chosen 
by  an  owner  are  major  decisions  in  long-run  timber 
management  strategies,  and  their  estimation  greatly 
influences  the  projected  results  of  other  analytical 
components  in  timber  supply  studies.  In  theory,  harvest 
patterns  would  be  chosen  to  maximize  discounted 
economic  returns  over  an  infinite  time  horizon,  if  the 
owner's  goal  were  profit  maximization.  However,  the 
planning  horizons  which  underlie  actual  harvest  deci- 
sions vary  greatly  among  owners  with  respect  to  short- 
run  versus  long-run  perspectives.  These  differences  can 
be  attributed  primarily  to  varying  financial  flow  needs 
and  uncertainty  regarding  financial  returns  from 
various  forest  management  strategies. 

While  marked  differences  among  harvest  patterns 
and  length  of  planning  horizons  exist  for  the  numerous 
types  of  forest  land  owners,  the  flexibility  in  timber  re- 
source adjustment  and  harvest  over  time  is  constrained 
to  a  large  degree  by  underlying  biological  and  economic 
processes  of  timber  production.  Distinguishing  between 
short-run  and  long-run  situations  in  aggregate  timber 
supply  analysis  is  somewhat  problematic,  especially  in 
view  of  the  long  production  periods  inherent  in  timber 
production.  In  an  economic  context,  the  short  run  is 
characterized  by  the  presence  of  a  fixed  production  fac- 
tor, while  in  the  long  run  all  factors  are  considered 
variable.  Lyon  and  Sedjo  (1983)  observe  that  in  forestry 
the  short  run  embodies  a  perspective  of  timber  produc- 
tion in  which  timber  manifests  many  characteristics  of  a 
nonrenewable  resource.  In  this  context,  the  question  of 
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short-run  supply  relates  to  the  rate  at  which  the  existing 
stock  is  harvested. 

Duerr  (1960)  also  proposes  temporal  characteri- 
zations or  conditions  for  timber  supply:  stock  supply, 
short-run  supply,  and  long-run  supply.  Stock  supply  of 
stumpage  arises  from  existing  timber  inventory  and, 
because  the  existing  growing  stock  inventory  is  also  the 
stumpage  factory,  stock-supply  responses  tend  to  over- 
shadow the  other  two  supply  responses.  Stock  supply  is 
theoretically  governed  by  the  opportunity  costs  of 
holding  the  stock  or  inventory  in  relation  to  alternative 
rates  of  return.  Short-run  supply  depends  upon  the  rate 
of  growth  of  the  growing  stock  and,  consequently,  upon 
the  cost  of  variable  inputs  that  might  be  used  during  a 
timber  rotation.  Land  is  the  main  fixed  input,  and  the 
short  run  may  consequently  span  several  decades.  The 
long  run  corresponds  to  the  period  when  all  inputs  are 
variable,  including  the  decision  whether  to  invest  in 
timberland. 

Although  this  outline  of  Duerr's  (1960)  timber  supply 
theory  is  somewhat  simplified,  it  still  suggests  the 
complexity  of  the  dynamics  of  timber  supply  over  time. 
With  aggregation  across  many  timber  management 
situations,  including  various  forest  stand  structures  and 
conditions  and  diverse  ownership  classes  with  varying 
time  frames  and  differing  management  perspectives,  it 
becomes  exceedingly  difficult  to  distinguish  between 
short  run  and  long  run  for  analytical  purposes.  The 
criterion  of  land  as  a  fixed  factor  of  production  over  a 
particular  time  period  may  apply  to  some  owners  but  not 
to  others.  Thus  the  precise  time  period  for  which  this 
factor  remains  fixed  may  be  difficult  to  identify  from  an 
aggregate  standpoint.  Aggregate  analysis  in  effect 
synthesizes  a  broad  range  of  rotations  which  vary 
according  to  a  wide  variety  of  species  and  management 
conditions. 

The  concept  of  the  short  run  as  the  transition  period, 
during  which  the  gradual  depletion  of  old-growth  timber 
and  the  merging  of  various  management  patterns  results 
in  a  "steady  state"  forest  situation  (Lyon  and  Sedjo 
1983),  represents  one  conceptual  approach  to  this  prob- 
lem. In  this  context,  the  long  run  represents  harvest 
levels  after  the  steady  state  has  been  achieved.  Such 
estimates  are  of  obvious  importance  with  respect  to  the 
long-term  status  of  the  timber  resource;  however,  this 
framework  is  difficult  to  apply  in  empirical  supply 
analyses.  Such  an  approach  is  addressed  further  in  the 
following  section  on  timber  investment  modeling.  The 
present  focus  is  upon  the  major  techniques  for  modeling 
harvest  flows,  with  less  emphasis  placed  upon  explicit 
criteria  for  long-term  investment  strategies. 

Modeling  of  aggregate  harvest  flows  has  generally 
been  addressed  through  two  distinct  analytical  ap- 
proaches: (1)  simulation  of  harvest  flows  based  upon 
historical  relationships  among  current  period  harvest, 
stumpage  prices,  and  some  proxy  for  the  opportunity 
costs  of  holding  timber  inventory  (Adams  et  al.  1982);  or 
(2)  optimization  techniques  such  as  mathematical  pro- 
gramming that  schedule  harvest  flows  so  as  to  result  in 
the  maximization  of  some  ownership  objective. 

Simulation  models  predict  or  display  the  conse- 
quences of  a  selected  range  of  management  alterna- 


tives; they  generally  provide  information  for  selection  of 
a  desirable  strategy  by  some  method  external  to  the 
simulation  model.  Duerr  et  al.  (1975)  classify  optimiza- 
tion or  "analytical"  models  as  those  that  evaluate  alter- 
natives in  order  to  find  an  optimum.  The  modifier 
"optimization"  is  used  in  place  of  Duerr  et  al.'s 
"analytical"  to  conform  to  popular  usage  of  these  terms 
in  the  literature. 

Simulation  techniques  are  often  used  in  positive 
analyses,  i.e.,  those  in  which  responses  of  the  timber 
resource  to  management  in  the  future  are  projected  as  a 
likely  continuation  of  response  patterns  that  occurred  in 
the  past.  Optimization  techniques  are  generally  associ- 
ated with  the  normative  modeling  concept  in  which 
goals  or  objective  functions  with  associated  constraints 
are  specified  or  prescribed.  The  distinction  between 
simulation  and  optimization  techniques  is  not  always 
straightforward  (Field  1978);  some  modeling  approaches 
incorporate  both  as,  for  example,  in  the  case  of  a 
mathematical  programming  technique  utilized  within 
the  broader  framework  of  a  simulation  model.  Despite 
their  simultaneous  use  for  some  modeling  problems, 
simulation  and  optimization  are  discussed  separately  as 
the  two  most  commonly  used  approaches  for  the  model- 
ing of  harvest  flows  in  aggregate  timber  supply 
analyses. 


SIMULATION  TECHNIQUES 

Simulation  has  been  used  to  represent  the  dynamic 
processes  for  a  wide  variety  of  situations  in  forest 
management  (Bare  1971).  In  simulating  large-scale 
forestry  systems,  key  parameters  and  variables  may  be 
systematically  varied  to  allow  monitoring  of  system  per- 
formance under  a  range  of  conditions.  A  dynamic  simu- 
lation model  is  a  representation  of  a  system  as  it  evolves 
over  time.  A  deterministic  simulation  model  does  not 
contain  variables,  otherwise  it  is  stochastic.  Difficulties 
inherent  in  model  validation,  including  statistical  valida- 
tion problems  for  stochastic  simulation  models,  are 
discussed  by  Reynolds  et  al.  (1981). 

With  respect  to  the  modeling  of  aggregate  harvest 
flows,  simulation  techniques  have  employed  various 
assumptions  regarding  harvest  trajectories.  For  ex- 
ample, baseline  simulations  have  been  obtained  by 
assuming  that  harvest  levels  would  increase  until  they 
equal  annual  growth  (USDA  Forest  Service  1973).  While 
earlier  simulation  efforts  tended  to  focus  upon  biological 
and/or  silvicultural  aspects  of  timber  management,  the 
emphasis  in  simulating  harvest  flows  appears  to  be 
shifting  more  toward  incorporation  of  economic 
variables  or  proxies  (Adams  1977,  McKillop  1967, 
Robinson  1974).  Simulation  of  the  response  of  private 
timber  harvest  to  changes  in  production  costs  and 
prices  was  used  in  the  1980  RPA  Assessment  (Mills  and 
Alig  1979,  Adams  and  Haynes  1980).  Prior  to  this, 
harvest  responses  to  economic  variables  had  generally 
been  determined  in  a  judgmental  fashion  for  major  na- 
tional timber  assessments. 

TAMM- The  Timber  Assessment  Market  Model 
(TAMM),  used  in  the  1980  RPA  Assessment,  is  a  compre- 
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hensive  system  of  analytical  techniques  that  models 
lumber,  plywood,  and  stumpage  markets  in  which  sup- 
ply and  demand  interact  to  determine  equilibrium  prices 
(Adams  and  Haynes  1980).  Its  stumpage  supply  modeling 
will  be  discussed  as  an  example  of  simulating  harvest 
flows  using  econometric  techniques,  which  are  statis- 
tical functions  estimated  by  using  economic  precepts  to 
guide  model  form. 

TAMM's  several  analytical  components  estimate 
future  prices,  consumption,  and  production  for  softwood 
lumber  and  plywood,  hardwood  lumber,  and  both  soft- 
wood and  hardwood  stumpage.  Annual  harvest  flows 
are  estimated  by  short-term  supply  functions,  for  which 
stumpage  prices  act  as  independent  variables.  These 
prices  are  obtained  through  a  linkage  of  the  stumpage 
and  product  markets;  this  is  achieved  by  a  complex 
simultaneous  solution  procedure  for  equations  involving 
equilibrium  prices  and  quantities  (Adams  and  Haynes 
1980).  In  effect,  stumpage  prices  are  estimated  each 
model  year  for  use  in  the  short-run  supply  functions,  and 
these  supply  functions  shift  over  time  as  other  independ- 
ent variables  change  (e.g.,  timber  inventory). 

After  the  annual  harvest  flows  are  estimated  in 
TAMM,  they  are  then  incorporated  into  the  TRAS  model 
(discussed  previously)  that  has  served  as  an  inventory 
projection  module  for  TAMM.10  The  TRAS  model  then 
calculates  the  timber  inventory  progression  for  the  next 
annual  cycle.  The  timber  inventory  at  the  end  of  that 
simulated  cycle  is  then  fed  into  the  short-term  supply 
functions  as  one  determinant  of  the  next  year's  harvest 
volumes.  This  cycle  then  repeats  over  the  time  period  of 
analysis. 

TAMM  is  one  of  the  most  comprehensive  models 
available  for  projecting  timber  market  activities,  partly 
because  of  the  attention  given  to  the  spatial  juxtaposi- 
tion of  the  various  supply  and  demand  regions  and  the 
linkage  of  the  stumpage  and  product  markets.  Regional 
competitive  advantages  attributable  to  differences 
across  supply  regions  in  the  cost  of  timber  production 
are  an  example  of  economic  behavior  that  can  be  simu- 
lated by  such  spatial  modeling.  Also,  stumpage  prices 
computed  by  a  simple  simultaneous  equations  pro- 
cedure using  the  equilibrium  values  from  the  product 
market  are  used  in  the  stumpage  supply  functions  to 
estimate  harvest  volumes. 

Domestic  supply  region  boundaries  used  in  TAMM 
modeling  for  the  1980  RPA  Assessment  (fig.  2)  were 
largely  dictated  by  availability  of  aggregate  forest 
inventory  data  (Adams  and  Haynes  1980).  Short-run 
supply  or  harvest  functions  used  in  the  1980  RPA  As- 
sessment were  estimated  for  the  eight  domestic  supply 
regions  separately  for  softwood  and  hardwoods  and 
four  ownership  categories.  For  example,  eight  supply 
functions  were  estimated  for  softwood  and  hardwood 
stumpage  supplies  of  four  owner  classes  in  the  South- 
Central  region,  one  of  the  domestic  supply  regions. 

"The  TRAS  model  itself  has  been  modified  to  estimate  harvest 
flows  via  stumpage-supply  equations  (Alig  et  al.  1982).  Up  to  six  in- 
dependent variables  can  be  specified  for  the  equations  (e.g.,  stump- 
age price)  with  the  user  providing  both  the  appropriate  variables  and 
coefficients. 


State  (1) 


Groups  of  timbersheds 
(512  possibilities) 

_\_  Timbershed  (10) 


Groups  of  owner  class 
(16  possibilities) 

Owner  class 

(5  per  timbershed) 

Administrative  units 
(64  for  entire  state) 

.  _\_  Grouped  resource  units 
( 1  20  for  entire  state) 

Basic  resource  units 
(200  for  entire  state) 


Figure  2.— Domestic  timber  supply  regions  used  in  the  1982  RPA 
Assessment  (US0A  Forest  Service  1982). 

The  TAMM  system  has  the  flexibility  for  harvest  flow 
estimations  for  more  disaggregated  supply  regions,  the 
major  limitation  being  the  existing  levels  of  statistical 
reliability  for  collection  and  organization  of  data  on 
timber  inventories  and  owner  supply  behavior.  This 
data  limitation  varies  for  different  geographical  areas, 
and  augmenting  data  sources  tends  to  be  a  long-term 
undertaking. 

A  typical  (50-year)  TAMM  projection  is  relatively  in- 
expensive because  of  TAMM's  computerized  nature, 
especially  given  the  relatively  low  cost  of  long-term  pro- 
jections by  the  TRAS  inventory  projection  module.  The 
speed  and  low  cost  of  TAMM  projections  make  the 
system  especially  useful  for  sensitivity  analyses.  The  in- 
fluence of  underlying  assumptions  can  be  tested  through 
several  projections  by  systematically  varying  input 
related  to  those  assumptions. 

The  short-term  harvest  model  embodied  in  TAMM  is 
based  on  the  following  precepts:  (1)  private  stumpage 
owners  will  vary  harvest  directly  with  stumpage  price, 
and  (2)  private  harvest  will  increase  with  available 
timber  inventory  levels  (Adams  and  Haynes  1980). 
Adams  and  Haynes  (1980)  reported  that  TAMM  ex- 
plained historical  harvest  flow  behavior  reasonably  well 
for  the  forest  industry  ownership  in  1980  RPA  Assess- 
ment analyses,  but  was  generally  not  as  successful  in 
modeling  private,  nonindustrial  harvest  flows.  TAMM's 
modeling  of  public  stumpage  supply  behavior  involves 
exogenously  determined  upper  harvest  bounds  arising 
from  allowable  cut  levels,  established  by  legislated 
policies. 

The  modeling  of  private  timber  harvest  by  TAMM 
employs  essentially  the  same  functional  forms  for  both 
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the  industrial  and  nonindustrial  owners.  The  relatively 
limited  success  in  modeling  nonindustrial  timber  harvest 
may  be  attributable  to  this  parallel  approach,  which 
does  not  directly  recognize  the  importance  of  non- 
monetary forest  returns  to  nonindustrial  owners.  Also, 
industry-related  considerations,  such  as  aspects  of  in- 
tegration in  industrial  firms,  are  ignored  by  TAMM's 
short-term  harvest  modeling,  which  is  based  directly  on 
price  and  timber  inventory  variables. 

TAMM's  difficulties  in  successfully  modeling  harvest 
behavior  of  the  diverse  class  of  nonindustrial,  private 
landowners  have  been  preceded  by  numerous  studies 
that  suggest  few  generalizations  can  be  made  regarding 
timber  management  patterns  by  nonindustrial  land- 
owners (e.g.,  Stoddard  1961,  Sizemore  1973,  McComb 
1975,  Clawson  1979).  Binkley  (1979,  1981)  has  proposed 
a  model  of  timber  harvest  behavior  of  nonindustrial, 
private  landowners  that  addresses  multiple  objectives 
of  ownership,  in  the  context  of  balancing  harvest  income 
and  personal  consumption  of  the  forest's  nontimber  out- 
puts. Binkley's  theoretical  model  is  based  on  utility  max- 
imization through  consumption  of  (1)  goods  and  services 
purchased  through  income  from  timber  harvests  and 
(2)  goods  based  on  nontimber  outputs  from  the  forest 
(e.g.,  recreation). 

Data  was  not  available  for  estimating  the  theoretical 
timber  supply  model  outlined  by  Binkley  (1981).  For  ex- 
ample, existing  data  for  even  the  relatively  small  study 
area  was  inadequate  for  estimating  technical  tradeoffs 
between  timber  and  nontimber  outputs.  A  "stochastic 
utility"  model  of  choice  was  developed  using  the  avail- 
able data,  with  a  maximum  likelihood  logit  estimator 
used  to  estimate  the  probability  of  harvest.  Because 
probability  of  harvest  was  analyzed  and  other  study 
conditions  were  variable,  Binkley's  (1981)  empirical 
results  cannot  be  directly  compared  to  those  for 
TAMM's  modeling  of  private  harvest  behavior.  How- 
ever, results  of  Binkley's  work  indicate  that  probability 
of  timber  harvest  is  strongly  influenced  by  stumpage 
price,  while  income  is  negatively  related  to  harvest 
probability.  Farmers  are  shown  to  be  more  likely  to 
harvest  timber  than  nonfarmers,  and  the  farmer  group 
also  appears  to  be  more  price  responsive. 

OPTIMIZATION  TECHNIQUES 

Optimization  models  represent  another  class  of  quan- 
titative analytical  techniques  frequently  applied  to  the 
modeling  of  harvest  flows.  While  simulation  procedures 
are  often  used  to  develop  a  range  of  alternatives  without 
specifying  a  single  most  desirable  solution  to  the  model- 
ing problem,  optimization  techniques  are  designed  to 
identify  a  particular  "optimal"  solution  that  satisfies  a 
set  of  predetermined  output  criteria. 

The  specific  aim  of  optimization  is  the  maximization 
(or  minimization)  of  the  value  of  an  objective  function 
that  incorporates  the  goal  of  the  modeling  effort  (e.g., 
maximizing  present  net  value  or  minimizing  costs).  For 
some  problems,  the  model  is  allowed  to  arrive  at  an  op- 
timal solution  in  an  unconstrained  fashion;  in  many  in- 
stances, however,  a  variety  of  constraints  are  specified 


as  conditions  for  achievement  of  the  optimal  solution. 
Optimization  models  for  timber  management  have 
generally  been  of  this  latter  variety.  Variables  in  the  ob- 
jective function  and  constraint  equations,  whose  value 
in  the  optimal  solution  are  to  be  determined,  are  re- 
ferred to  as  decision  variables. 

Optimization  techniques  applied  to  estimate  harvest 
flows  in  timber  supply  analyses  have  used  a  variety  of 
biological,  economic,  and  other  specified  objective  func- 
tion criteria,  as  well  as  a  wide  range  of  model  con- 
straints (e.g.,  Navon  1971).  Johnson  and  Scheurman 
(1977),  Bell  (1977),  Field  (1978),  Marty  (1979),  and  Hann 
and  Brodie  (1980)  review  optimization  tools  for  harvest 
estimation  in  even-aged  forests,  and  Hann  and  Bare 
(1979)  discuss  techniques  for  uneven-aged  cases.  Op- 
timization techniques  evaluate  harvest  decisions  as  part 
of  a  total  timber  management  package  that  includes 
other  management  practices;  these  management  alter- 
natives are  prespecified  external  to  the  model.  This  im- 
portant interdependency  of  the  management  variables 
in  terms  of  joint  optimization  will  be  discussed  further  in 
the  following  section. 

Among  these  optimization  approaches,  linear  pro- 
gramming (LP)  has  been  utilized  most  frequently  for  the 
modeling  of  harvest  flows.  This  in  part  reflects  the  abil- 
ity of  LP  to  address  large-scale  modeling  problems  and 
its  degree  of  computational  efficiency  relative  to  other 
optimization  procedures.  A  primary  assumption  of  any 
LP  analysis  is  that  the  problem  being  modeled  is  capable 
of  being  formulated  in  such  a  way  that  an  optimal  solu- 
tion to  that  problem  may  be  found.  A  recent  application 
of  LP  utilized  by  the  USDA  Forest  Service  for  National 
Forest  lands  is  discussed  below  as  an  example  of  the  op- 
timization framework  for  modeling  harvest  flows. 

FORPLAN.— The  Forest  Planning  (FORPLAN)  model 
is  a  linear  programming  (LP)  package  that  has  been  ap- 
plied to  the  complex  optimization  problems  associated 
with  National  Forest  planning  mandated  by  the  National 
Forest  Management  Act  of  1976  (Johnson  et  al.  1980). 
FORPLAN  is  actually  a  software  package  that  serves  as 
a  linear  programming  matrix  generator,  an  interface  to 
the  UNIVAC  1100  FMPS  LP  solution  algorithm,  and  a 
report  writer  for  describing  the  LP  solution.  The  LP 
matrix  generator  converts  data  input  into  a  matrix  of 
rows  and  columns.  This  matrix  is  then  used  in  the 
simultaneous  solution  of  land  allocation,  land  manage- 
ment or  investment  scheduling  activities,  and  output 
mix.  The  report  writer  can  produce  final  output  in  a 
variety  of  specified  formats.  FORPLAN  essentially 
evolved  from  a  timber  management  model  known  as  the 
Multiple  Use-Sustained  Yield  Calculator  (MUSYC)  and 
retains  emphasis  on  timber  analysis  capabilities; 
however,  prescriptions  and  output  scheduling  for  a 
variety  of  resources  can  be  examined  as  well,  if  relevant 
biological  and  economic  information  is  available. 

FORPLAN  can  construct  an  LP  based  on  either  Model 
I  or  Model  II  structures  for  timber  harvest  analysis 
(Johnson  and  Schuerman  1977).  Basically,  these  two 
model  structures  differ  in  the  manner  in  which  they 
define  timber  management  variables  and  handle  multi- 
ple harvests  within  the  planning  horizon.  In  the  Model  I 
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structure,  a  variety  of  management  activities  can  occur 
on  a  particular  land  area  over  the  entire  planning 
horizon,  with  the  specified  area  retaining  its  identity 
throughout  the  process.  In  the  Model  II  approach,  in- 
dividual parcels  of  land  are  not  kept  intact  through  time 
because  harvested  acres  are  recombined  into  new  man- 
agement units.  Model  I  is  more  effective  in  keeping  track 
of  location  on  the  ground,  while  Model  II  is  more  con- 
ducive to  minimizing  model  size.  In  terms  of  earlier 
timber  management  models,  FORPLAN  can  employ  the 
Model  I  structure  of  the  Timber  Resources  Allocation 
Model  (Timber  RAM)  developed  by  Navon  (1971)  or  the 
Model  II  structure  of  the  MUSYC  model  noted  above. 

Joyce  et  al.  (1983)  review  the  history  and  foundations 
of  linear  programming  applications  in  the  natural 
resources  area,  including  the  use  of  FORPLAN  as  a 
multiresource  analysis  tool  for  relatively  small  areas. 
The  timber  harvest  analysis  capabilities  of  FORPLAN 
are  similar  to  other  LP  models  used  for  optimization 
analyses  (e.g.,  Ware  and  Clutter  1971),  and  they  are 
discussed  only  briefly  in  this  report.  Johnson  and 
Schuerman  (1977)  compare  different  types  of  LP  models 
that  utilize  Model  I  or  Model  II  frameworks,  discuss  the 
different  solution  techniques  that  have  been  applied, 
and  also  review  quadratic  versions  of  Model  I  and 
Model  II  (e.g.,  Walker  1971). 

Basically,  an  LP  model  such  as  FORPLAN  can  max- 
imize or  minimize  a  linear  objective  function  subject  to 
linear  constraints.  Ten  different  possible  forms  of  objec- 
tive functions  for  FORPLAN  are  listed  by  Joyce  et  al. 
(1983).  FORPLAN  can  also  generate  a  piecewise  approx- 
imated downward-sloping  demand  curve  for  timber. 
Fixed  prices  are  assumed  for  all  other  outputs.  At  this 
time,  costs  can  be  assigned  per  acre  for  timber-related 
management  prescriptions,  and  per  unit  output  for  all 
outputs  (including  timber).  It  is  anticipated  that  the 
capability  to  assign  costs  per  acre  for  all  management 
prescriptions  will  soon  be  available  in  FORPLAN.  An  op- 
tion to  include  fixed  costs  is  also  available  in  the  model. 

Output  constraints  or  output  targets  are  an  important 
part  of  the  LP  formulation  within  FORPLAN.  These 
targets  set  minimum  or  maximum  levels  of  outputs  to  be 
obtained  in  the  LP  solution  and  actually  drive  the  model 
in  some  instances.  Inclusion  of  certain  specific  targets  is 
directly  mandated  by  the  NFMA  regulations.  The  fifteen 
kinds  of  tables  and  graphs  which  constitute  the  outputs 
of  the  FORPLAN  model  include  the  objective  function 
value,  constraint  output,  land  allocation  areas,  harvest 
reports  (area  and  volume),  economic  reports  (costs, 
prices,  revenues,  present  net  worth,  present  net 
benefit),  regeneration  levels,  and  cultural  treatments 
(Johnson  et  al.  1980). 

FORPLAN  allows  considerable  flexibility  in  the  user 
designation  of  analysis  areas  and  management  prescrip- 
tions. Analysis  areas  are  assumed  to  respond  the  same 
way  to  a  given  management  prescription,  regardless  of 
where  they  occur.  Analysis  areas  are  defined  by  six 
levels  of  identification:  levels  1-3,  user-specified  land 
characteristics  (60  categories  for  each  level);  level  4, 
working  groups  (9  categories);  level  5,  land  classes  (15 
categories);   level  6,   existing  vegetation  classes  (60 


categories).  Every  alternative  management  prescription 
specified  in  FORPLAN  must  apply  to  a  given  analysis 
area.  In  addition,  alternative  prescriptions  can  be  in- 
cluded to  apply  to  the  regeneration  classes  created  by 
prescriptions  with  harvest  practices. 

Production  coefficients  associated  with  specific  treat- 
ments and  land  types  are  input  into  the  model  as  yield 
tables.  They  generally  represent  the  number  of  units  per 
acre  of  output  resulting  from  a  given  treatment  applied 
to  a  given  analysis  area.  Prices  and  costs  are  entered  in 
the  same  manner  as  the  yield  tables  noted  above.  The 
time  frame  that  applies  to  management  prescriptions  as 
well  as  to  resultant  outputs  is  quite  variable;  it  may  con- 
sist of  up  to  30  time  periods  of  from  1  to  20  years  each. 

FORPLAN  model  formulations  may  involve  con- 
straints relating  to  any  or  all  of  the  following:  area  and 
volume  control,  harvest  flow  (e.g.,  nondeclining  yield), 
ending  inventory,  management  emphasis  and  intensity, 
and  cultural  treatments.  These  constraints  may  be 
represented  by  numerous  configurations,  and  generali- 
zations are  inappropriate  in  this  report.  For  example, 
the  determination  of  which  management  alternative  to 
examine  within  the  FORPLAN  model,  and  the  derivation 
of  the  associated  biological  response  (e.g.,  yield  tables) 
are  for  the  most  part  forest-specific. 

Apple  (1982)  surveyed  FORPLAN  users  in  7  of  the  9 
USDA  Forest  Service  regions  regarding  the  utility  of 
FORPLAN  for  planning  activities  on  National  Forest 
lands.  Most  respondents,  who  consisted  primarily  of 
forest  planners  and  operations  research  analysts,  indi- 
cated that  between  one  and  twelve  initial  model  runs 
were  required  to  obtain  a  usable  solution.  Approxi- 
mately 60%  of  respondents  were  utilizing  Model  I  for 
analytical  runs;  the  remaining  40%  were  employing 
Model  II.  Model  I  was  regarded  as  less  costly  and  easier 
to  track,  while  Model  II  was  considered  to  be  better 
suited  for  forest  management  data  and  a  more  flexible 
analytical  tool. 

Among  some  of  the  unexpected  benefits  arising  from 
the  use  of  FORPLAN,  planners  and  analysts  felt  that  it 
provided  a  better  evaluation  of  inventory,  an  increased 
understanding  of  the  role  of  economics  in  forest  manage- 
ment, better  coordination  of  management  functions,  and 
a  better  overall  understanding  of  the  forest  management 
system.  Problems  listed  most  frequently  included  the 
need  for  more  direction  in  model  use,  the  lack  of  or  poor 
quality  of  input  data,  and  lack  of  time  for  analysis,  par- 
ticularly with  respect  to  the  performance  of  sensitivity 
analysis  for  model  results.  Users  indicated  that  the 
major  areas  of  information  needs  included  economic 
data  and  yield  data  for  timber  and  other  resources. 

Although  the  FORPLAN  model  has  primarily  been 
applied  at  the  individual  National  Forest  level,  the  rela- 
tionship of  forest-level  FORPLAN  implementation,  as 
well  as  other  optimization  techniques,  to  aggregate 
timber  supply  analysis  is  of  major  importance  within  the 
context  of  this  report.  Aggregate  timber  supply  analysis 
for  national  assessments  necessarily  involves  a  critical 
linkage  of  potential  harvest  flows  from  public  and  pri- 
vate sources;  however,  these  two  sources  are  typically 
modeled  with  varying  detail.  Public  timber  harvest  flows 
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in  the  most  recent  national  timber  assessment  were 
essentially  exogenous  estimates,  while  private  harvests 
were  estimated  using  TAMM's  simulation  modeling. 
Questions  have  arisen  regarding  the  possible  expanded 
use  of  FORPLAN  output  from  forest-level  modeling  for 
national  assessments,  or  use  of  a  modeling  framework 
similar  to  FORPLAN's  but  at  a  broader  regional  level. 

These  concerns  regarding  the  potential  applicability 
of  optimization  techniques  for  modeling  public  harvest 
flows  at  aggregate  levels  (e.g.,  National  Forest  region) 
are  also  relevant  to  the  modeling  of  aggregate  private 
harvest  flows  in  some  regions.  The  multilevel  or  hierar- 
chical decision  problem  of  different  planning  levels  (e.g., 
national,  regional,  and  Forestflocal)  is  difficult  to  solve 
because  of  the  need  to  recognize  (1)  that  variables  under 
control  of  policy  makers  are  distinct  from  those  under 
control  of  the  micro  unit  or  behavioral  decision  maker 
and  (2)  conflicting  objective  functions  may  exist  for  dif- 
ferent owners  throughout  the  hierarchy  (Candler  et  al. 
1981).  The  goals  of  economic  units  can  also  often  be 
stated  more  explicitly  and  unambiguously  at  the  micro 
level  than  at  the  macro  level.  Specification  of  appropri- 
ate or  realistic  objective  functions  and  constraints  at  ag- 
gregate levels  for  the  various  owner  classes  could  be  dif- 
ficult tasks,  particularly  for  the  diverse  nonindustrial 
class.  Aggregation  across  such  diverse  owner  and 
timber  resource  conditions  may  also  lead  to  relatively 
large  model  size  and  cost  if  needed  modeling  detail  is  to 
be  achieved  (Ashton  et  al.  1980).  Little  has  been 
reported  in  the  literature  regarding  the  feasibility  of 
modeling  aggregate  timber  harvest  flows  solely  within 
an  optimization  framework;  however,  mathematical  pro- 
gramming has  been  utilized  in  somewhat  analogous 
roles  to  model  agricultural  supply  activities  (e.g., 
Meister  and  Nicol  1975). 

Another  potential  difficulty  with  respect  to  expanded 
use  of  optimization  modeling  involves  questions  regard- 
ing computational  efficiency  that  stem  from  the  neces- 
sity of  arriving  at  a  single  optimal  level  in  order  to 
satisfy  the  production  of  a  desired  output  (e.g.,  timber). 
As  an  alternative  to  addressing  the  large  computational 
requirements  of  the  LP  process,  Hoganson  (1981)  de- 
scribes a  method  whereby  the  LP  problem  is  formulated 
as  a  simulation.  This  simulation  model  is  driven  by 
shadow  prices  (dual  variables  in  an  LP  formulation)  and 
the  iterative  process  continues  until  harvest  levels  are 
driven  to  a  level  close  to  that  specified  in  the  regional  LP 
problem.  The  technique  has  the  potential  for  significant 
cost  reductions  in  the  analysis  of  large  LP  problems. 


OTHER  APPROACHES 

Goal  programming  (GP)  is  a  particular  form  of  LP  used 
to  find  a  solution  to  a  resource  allocation  problem 
involving  several  objectives,  subject  to  a  set  of  linear 
constraints.  Depending  on  the  type  of  formulation,  all 
goals  may  be  considered  simultaneously  in  a  composite 
(and  single)  objective  function  or  sequentially  in  a  series 
of  objective  functions  (Bell  1976,  Field  1978,  Field  et 
al.  1980).  The  choice  variables  are  deviational  ones- 


showing  overachievement  or  underachievement  of  the 
specified  goal  levels  of  output.  As  with  FORPLAN,  most 
GP  models  have  been  applied  to  relatively  small  areas. 
Schuler  et  al.  (1977)  suggest  that  the  largest  application 
problem  involves  determination  of  the  technical  coeffi- 
cients (i.e.,  the  amount  of  output  that  can  be  anticipated 
from  a  given  input). 

Another  approach  based  on  mathematical  program- 
ming is  Berck's  (1979)  empirically  based  model  of  long- 
run  timber  supply  behavior.  Berck  tested  the  hypothesis 
that  private  owners  cut  their  timber  at  a  rate  exceeding 
an  optimal  one  based  on  alternative  rates  of  returns. 
This  analysis  involved  estimating  the  rate  of  time 
preference  indicated  by  the  interest  or  discount  rate  im- 
plied in  historical  management  of  Douglas-fir  by  private 
landowners,  assuming  they  maximized  present  value  of 
profits.  Berck  found  that  cutting  rates  for  private  land- 
owners indicated  they  were  discounting  the  future  at  a 
rate  much  lower  than  the  rate  of  return  available  for 
other  private  investments  and  much  lower  than  that  pro- 
posed by  Arrow  (1976)  as  a  social  discount  rate.  Berck's 
assumption  of  rational  price  expectations  by  private 
landowners  in  order  to  maximize  profits  over  a  175-year 
period  implies  knowledge  and  behavior  that  may  not  be 
representative  of  many  private  forest  land  owners.  Non- 
timber  benefits  from  the  private  forests  were  not  in- 
cluded in  this  modeling  of  private  forest  owner  behavior. 

TREES.— An  example  of  a  comprehensive  approach  to 
modeling  aggregate  harvest  flows  is  the  Timber 
Resources  Economic  Estimation  System  (TREES)  (John- 
son et  al.  1975,  Tedder  et  al.  1980)  model.  This  model 
contains  a  number  of  harvest  estimation  options,  includ- 
ing variable  harvest  scheduling  methods  that  use  a 
multiple  iteration  binary  search.  The  aggregation 
scheme  of  the  TREES  model  is  noteworthy  for  its  sophis- 
tication and  flexibility. 

The  TREES  model  was  originally  developed  for  a 
comprehensive  study  of  future  timber  supply  for  the 
state  of  Oregon  (Beuter  et  al.  1976).  The  model  projects 
different  harvest  levels  in  response  to  varying  assump- 
tions about  land  use  changes,  growth  rates,  trends  in 
silvicultural  intensity,  harvest  control  policies,  and 
utilization  efficiencies.  TREES  is  capable  of  (1)  project- 
ing timber  volumes  resulting  from  either  natural  stand 
development  or  in  response  to  specified  treatment  activ- 
ities, (2)  simulating  harvests  at  numerous  administrative 
levels,  and  (3)  considering  multiple  stocking  levels  in  in- 
ventory projections.  The  TREES  model  can  also  project 
timber  inventories  for  both  even-aged  and  uneven-aged 
forest  structures.  The  key  to  these  capabilities  is  the 
availability  of  yield  tables.  TREES  does  not  provide  such 
tables;  it  uses  them  to  make  timber  inventory  projec- 
tions. Input  to  the  TREES  model  includes  base  growth 
and  yield  information  via  growth  functions,  tables,  or  a 
combination  of  the  two.  These  are  keyed  according  to 
forest  type,  site,  and  acreage  of  each  class  by  owner  and 
administrative  unit. 

Seven  harvest  scheduling  options  allow  the  user  to 
simulate  harvest  levels  under  current  public  and  private 
management  practices  and  to  simulate  new  harvest 
schedules  to  evaluate  policy  alternatives.  These  harvest 
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schedules  are  developed  through  either  single  (fixed 
harvest  schedule)  or  multiple  (variable  harvest  sched- 
ule) iterative  computer  runs.  The  three  fixed  methods  do 
not  use  an  iterative  search  to  find  the  harvest  level,  so 
no  optimization  occurs  (Schmidt  and  Tedder  1980). 

The  three  fixed  methods  can  schedule  harvest  in  ac- 
cordance with  the  following  rules: 

1.  Area  control:  Number  of  acres  to  be  harvested  is 
determined  by  dividing  total  available  acres  by  a 
prespecified  rotation  length. 

2.  Percent  of  inventory:  A  proportion  (fixed  or 
variable)  of  the  available  inventory  is  to  be 
harvested  each  period  according  to  the  harvest 
priority  chosen  (e.g.,  oldest  age  first). 

3.  Absolute  amount:  A  specified  timber  volume, 
selected  on  the  basis  of  the  harvest  priority  rule,  is 
to  be  harvested  each  period. 

The  four  variable  or  multiple  iteration  procedures  em- 
body principles  of  several  other  harvest  scheduling 
models.  In  development  of  these  timber  flows,  a  binary 
search  technique  estimates  the  following: 

1.  An  even  flow  of  timber  volume  sustainable  over  an 
entire  specified  projection  interval,  similar  to  the 
approach  used  in  the  SIMAC  model  (Sassaman  et 
al.  1972).  Alternatively,  an  even  flow  of  timber 
volume  is  estimated  for  each  harvest  period  within 
the  projection  interval.  When  this  technique  is 
used,  the  allowable  cut  is  recalculated  at  the 
beginning  of  each  harvest  period  throughout  the 
entire  projection  period.  This  approach  embodies 
the  ideas  of  the  SORAC  model  developed  by  Chap- 
pelle  and  Sassaman  (1968). 

2.  An  even-flow  as  in  item  1,  except  that  the  variable 
of  interest  is  a  linear  function  of  volume  instead  of 
volume  itself.  Such  volume-dependent  variables  in- 
clude employment  or  gross  revenue. 

3.  Harvest  rates  based  on  the  economic  criterion  of 
present  net  worth  (PNW)  in  a  constrained  maximi- 
zation framework.  This  approach  is  similar  to  the 
ECHO  model  proposed  by  Walker  (1971),  where  the 
volume  harvested  is  determined  for  a  timberland 
owner  facing  a  negatively  sloping  demand  curve 
that  reflects  quantity  and  price  interactions  and 
market  power. 

4.  Harvest  schedules  to  maximize  present  net  benefit 
(PNB),  which  is  the  net  discounted  difference  be- 
tween a  measure  of  the  consumer's  "willingness  to 
pay"  and  total  cost.  This  harvest  scheduling 
method  is  analogous  to  item  3,  except  that  it  max- 
imizes PNB  rather  than  PNW,  and  the  former  is 
relatively  more  difficult  to  estimate.  The  output 
projections  are  given  by  timbershed  and  decade 
for  timber  volume  and  growth,  volume  and  value  of 
harvests,  volume  and  value  by  products,  and  pay- 
ments in  lieu  of  taxes. 

One  of  the  strengths  of  the  TREES  model  is  its 
elaborate  control  structure  when  it  is  used  to  aggregate 
the  timber  harvests.  This  control  structure  organizes 
resource  units  in  such  a  way  that  the  model  can  keep 
separate  and  monitor  a  large  number  of  resource  units 


(fig.  3).  The  hierarchical  aggregation  scheme  of  TREES 
can  currently  aggregate  information  for  up  to  eight  dif- 
ferent levels. 

The  basic  resource  unit  for  TREES  is  the  smallest 
piece  of  land  for  which  a  beginning  timber  inventory 
may  be  entered.  Examples  of  resource  units  from  the 
Oregon  study  are  even-aged  units  represented  by 
volume  per  acre  by  age  class,  site  class,  and  species  for 
western  Oregon,  and  uneven-aged  units  represented  by 
number  of  trees  per  acre  by  diameter  class  in  eastern 
Oregon.  The  next  level  of  aggregation  is  the  grouped 
resource  unit,  a  collection  of  basic  resource  units  all 
possessing  a  similar  resource  quality  (e.g.,  a  group  of 
forest  stands  of  single  species  and  site  class  and  all 
being  managed  under  the  same  harvest  method).  The 
grouped  resource  units  or  management  units  are  the 
lowest  level  at  which  timber  harvests  and  projections 
are  estimated. 

Grouped  resource  units  can  be  combined  to  form  ad- 
ministrative units.  Administrative  units  are  sensitive  to 
institutional  constraints  and  closely  follow  the  level  at 
which  federal  harvests  are  regulated.  Harvests  are  gen- 
erally calculated  at  this  level,  and  each  projection  has 
its  own  assumptions  about  land  use  changes,  manage- 
ment intensities,  and  harvest  control  procedures. 

Administrative  units  may  be  combined  within  an 
owner  class,  and  owner  classes  can  be  grouped  within  a 
timbershed.  The  tirribersheds  can  then  be  aggregated 
into  groups  within  a  state  or  a  multistate  region.  For  the 
Oregon  study,  the  state  was  divided  into  a  collection  of 
timbersheds  that  were  defined  as  the  geographic  area 
from  which  the  primary  needs  of  a  wood  processing 
center  were  supplied.  A  major  advantage  of  the  flexible 
aggregation  scheme  of  the  TREES  model  is  that  harvest 
quantities  can  be  determined  at  any  level  from  the 
grouped  resource  unit  up  to  the  state  level. 

The  TREES  model  was  also  utilized  in  a  Pacific  North- 
west Regional  Commission  Study  of  timber  supply  in 
Washington,  Oregon,  and  Idaho  (Rahm  1981).  In  the 
baseline  analysis,  the  harvests  of  public  sector  owners 
are  assumed  to  continue  along  the  current  planned  path 
through  the  next  century.  Private  timber  owners  are 
assumed  to  respond  fully  to  market  forces,  harvesting 
timber  at  a  rate  that  maximizes  the  present  net  value  of 
their  timber  income  stream  in  every  decade  for  the  next 
100  years.  In  an  "optimizing  simulation,"  the  TREES 
model  selects  harvest  schedules  in  conjunction  with 
downward-sloping  demand  curves  (net  of  planned  public 
harvest)  for  delivered  logs  in  two  sub  regions.  The 
assumed  objective  function  is  the  maximization  of  pres- 
ent net  benefit;  harvests  are  selected  in  each  time 
period  so  that  the  discounted  net  revenue  from  the  last 
(marginal)  unit  of  wood  harvested  is  just  equal  to  what 
that  unit  would  yield  if  allowed  to  grow  and  be  har- 
vested in  the  next  time  period. 

Because  the  flow  of  delivered  logs  is  the  focus  of  this 
equilibrium  analysis,  the  costs  of  removing  and  hauling 
the  timber  determine  the  corresponding  stumpage 
sources.  Per  unit  costs  by  age  class  for  felling,  bucking, 
skidding,  and  hauling  are  required  for  final  harvests 
and  thinnings.  Price  adjustment  factors  are  utilized  to 
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Figure  3.— Control  structure  of  the  TREES  model  used  in  Oregon  study  to  aggregate  timber 
supply  information  (adapted  from  Johnson  etal.  1975). 


reflect  the  differing  values  of  logs  harvested  from  stands 
of  different  ages. 

Following  this  study,  Rahm  (1981)  pointed  out  the  fol- 
lowing model  limitations  with  regard  to  timber  inventory 
projection  capabilities  of  the  TREES  model:  (a)  certain 
combinations  of  timber  management  activities  cannot  be 
simulated,  necessitating  the  use  of  broad  average  man- 
agement practices;  (b)  the  time  path  of  changes  in  timber 
management  intensity  seems  unrealistic  because  the 
proportion  of  total  acres  in  a  class  rather  than  the  pro- 
portion harvested  is  shifted  after  harvest;  (c)  a  manage- 
ment intensity  change  (e.g.,  fertilization  in  middle-aged 
natural  stands)  cannot  begin  until  acres  are  harvested; 
and  (d)  acres  shifted  from  nontimber  management  to 
timber  management  are  assumed  to  have  the  same 
average  age  structure  and  volume  per  acre  as  acres 
already  in  the  management  base. 

While  the  TREES  model  considers  the  cost  of  remov- 
ing timber  from  the  stump  and  transporting  it  to  market, 
Rahm  (1981)  indicates  that  the  model  does  not  properly 
specify  the  changes  in  costs  associated  with  shifts  in  the 
geographic  concentration  of  harvests  over  time.  The 
model  selects  timber  for  harvest  in  areas  that  yield  the 


highest  average  net  value  based  on  the  harvest  costs. 
This  constrains  timber  harvests  solely  to  a  subset  of 
analysis  areas  (i.e.,  grouped  resource  units)  in  a 
stepwise  manner,  rather  than  having  them  more  widely 
distributed  geographically.  Because  of  this  inter- 
subregional  harvest  optimization  procedure,  the  result- 
ing harvest  flows  may  not  necessarily  be  a  realistic  in- 
dicator of  the  actual  location  of  future  timber  harvests 
within  a  region. 

Rahm  (1981)  further  points  out  that  the  TREES  model 
tends  to  be  expensive  to  operate  because  its  "all- 
purpose"  orientation  results  in  a  size  and  complexity  of 
the  model  that  tend  to  make  individual  runs  costly.  Also, 
because  it  is  an  extremely  open-ended  tool,  a  consider- 
able amount  of  time  is  required  for  the  user  to  become 
familiar  with  its  operation.  However,  once  basic  model 
runs  are  set  up,  revisions  can  be  made  relatively  quickly 
and  large  data  sets  can  be  accommodated  satisfactorily. 

The  TREES  model  also  does  not  necessarily  calculate 
economic  optima,  but  produces  approximations  because 
of  assumptions  such  as  those  listed  by  Schmidt  and  Ted- 
der (1980).  For  example,  one  assumption  is  equivalent  to 
ignoring  the  opportunity  cost  of  delaying  future  harvest. 
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Rahm  (1981)  suggests  that  such  deviations  from 
economic  optima  may  not  be  significantin  some  applica- 
tions, especially  in  view  of  the  degree  of  error  inherent 
in  some  data  that  is  used. 


INTEGRATION  OF  HARVEST  FLOWS  AND  TIMBER 
INVENTORY  PROJECTION  MODELING 

Linkage  of  the  four  analysis  components  in  aggregate 
timber  supply  studies  has  been  improved  in  recent 
years,  particularly  with  respect  to  interactions  between 
biological  and  economic  processes.  However,  the  mod- 
eling of  harvest  flows  and  investment  has  seldom  been 
as  detailed  as  that  for  timber  growth  in  past  studies. 
Future  harvest  patterns  were  sometimes  estimated  on 
the  basis  of  simplified  assumptions  or  scenarios,  with  no 
explicit  consideration  given  to  linkages  with  trends  in 
important  economic  variables  (e.g.,  prices).  Despite  re- 
cent efforts  to  devote  more  attention  to  economic  proc- 
esses, questions  still  remain  regarding  the  appropriate 
balance  of  modeling  detail  for  biological  and  economic 
components  in  aggregate  timber  supply  studies. 

While  more  effective  integration  of  biological  and 
economic  models  has  slowly  occurred,  the  flexibility  of 
such  analytical  packages  for  examining  timber  supply 
policy  questions  is  still  rather  limited.  For  example,  the 
TAMM  harvest  flow  modeling  discussed  previously 
relied  on  timber  inventory  projections  by  the  TRAS 
model,  which  is  relatively  inflexible  in  its  capabilities  to 
simulate  timber  management  alternatives.  This  rigidity 
in  the  biological  simulation  component  diminishes  the 
usefulness  of  the  overall  analytical  package  for  timber 
supply  analysis. 

Models  that  are  independently  developed  and  later 
linked  together  often  require  major  adjustments  in  order 
to  permit  efficient  and  effective  joint  operation.  A  major 
problem  may  occur  when  biological  and  economic 
models  must  interface  at  a  common  level  of  aggregation. 
For  example,  the  integrity  of  TRAS's  biological  modeling 
may  be  diminished  if  it  is  applied  to  large  geographical 
regions  which  serve  as  the  bases  for  harvest  flow  equa- 
tions or  other  economic  modeling.  Aggregation  across 
many  silvicultural  conditions,  such  as  site  classes,  stock- 
ing levels,  and  stand  structures,  may  more  adversely 
impact  model  validity  for  the  biological  modeling  com- 
ponents than  for  those  involved  with  modeling  economic 
processes. 

SUMMARY 

The  TAMM  model  (Adams  and  Haynes  1980)  repre- 
sents the  most  comprehensive  approach  for  integrating 
short-run  private  harvest  behavior  based  on  historical 
patterns  with  long-term  timber  management  strategies. 
Econometrically  derived  short-run  harvest  flow  equa- 
tions are  shifted  through  time  in  a  TAMM  simulation,  as 
in  the  1980  RPA  Assessment.  Stumpage  price  and  timber 
inventory  levels  are  the  independent  variables  in  the 
harvest  flow  equations,  which  are  somewhat  ad  hoc  in 


nature.  The  TAMM  approach  has  the  strongest  ex- 
planatory power  with  regard  to  forest  industry 
behavior,  but  it  generally  encounters  major  difficulties 
in  attempting  to  explain  the  harvesting  behavior  of 
nonindustrial,  private  landowners.  Other  empirical 
studies  of  nonindustrial  harvesting  behavior  include 
Binkley's  (1981)  dichotomous  procedure  for  estimating 
the  probability  of  harvest  in  a  certain  year.  Although 
substantial  data  problems  were  encountered,  stumpage 
prices  were  shown  to  strongly  influence  the  probability 
of  harvest. 

The  Forest  Planning  (FORPLAN)  model  (Johnson  et  al. 
1980)  represents  the  application  of  optimization  tech- 
niques to  the  modeling  of  harvest  flows.  This  linear 
programming  (LP)  model  evolved  from  previous  LP  ap- 
proaches for  timber  harvest  scheduling,  but  FORPLAN 
differs  from  its  predecessors  in  its  multiresource  dimen- 
sions. Applicability  of  the  FORPLAN  model  and  other 
optimization  approaches  to  the  modeling  of  aggregate 
harvest  flows  is  hindered  in  part  by  difficulties  in 
specifying  objective  functions  and  constraint  equations 
for  such  analyses.  Other  mathematical  programming  ap- 
proaches include  Berck's  (1979)  empirically  based 
model  of  long-term  timber  harvesting  patterns.  This 
modeling  of  the  intertemporal  dependency  of  harvest 
and  management  decisions  involves  estimation  of 
private  landowner  rate-of-time  preference  as  a  behav- 
ioral parameter. 

The  Timber  Resources  Economic  Estimation  System 
(TREES)  Qohnson  et  al.  1975,  Tedder  et  al.  1980)  is  a 
relatively  elaborate  harvest  scheduling  model  that  has 
seven  harvest  scheduling  options.  The  variable  harvest 
scheduling  methods  (e.g.,  present  net  worth  algorithm) 
use  a  multiple  iteration  binary  search,  and  several  op- 
tions can  be  applied  to  either  even-aged  or  uneven-aged 
stands.  TREE's  aggregation  scheme  has  considerable 
potential  for  expanded  adoption  in  aggregate  timber 
supply  studies,  as  it  utilizes  a  comprehensive,  hier- 
archial  control  structure.  The  aggregation  control  struc- 
ture allows  harvest  quantities  to  be  determined  at  any 
level  from  a  grouping  of  basic  resource  units  up  to  the 
state  level. 


4.  TIMBER  INVESTMENTS 

Investments  to  enhance  the  productive  capacity  of 
forests  involve  utilization  of  inputs  which  influence  the 
growth  of  the  inventory  over  time  and  future  harvest 
levels.  The  process  of  timber  production  may  require  in- 
put expenditures  (e.g.,  artificial  regeneration)  well  in 
advance  of  when  the  effects  of  those  inputs  are  realized 
(i.e.,  gestation  period).  Capital  theory  deals  with  the 
time  path  of  capital  accumulation,  and  harvest  flows 
discussed  in  the  preceding  section  are  part  of  the 
capital  formation  process  in  forestry,  representing 
disinvestment  or  liquidation  of  accumulated  timber 
capital  stocks. 

Although  timber  investments  are  important  for  long- 
run  timber  supply,  analyses  of  aggregate  timber  invest- 
ment patterns  have  been  rather  limited.  Decisions  to 
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undertake  intertemporally  linked  timber  investments 
such  as  planting,  thinning,  or  fertilization  theoretically 
depend  on  anticipated  physical  yields,  opportunity  costs 
of  such  investments,  and  stumpage  prices.  However, 
with  the  diversity  of  owner  objectives  (even  within 
owner  classes)  and  the  widespread  deficiencies  in 
technical  production  and  behavioral  data,  realistic 
modeling  of  long-term  timber  management  behavior 
remains  one  of  the  more  difficult  tasks  facing  timber 
supply  analysts. 

The  theory  of  aggregate  forest  investment  is  not  well 
developed,  and  almost  all  applied  timber  investment 
analyses  have  been  conducted  at  the  stand  level.  While 
there  has  been  considerable  statistical  testing  of  ag- 
gregate investment  theories  with  regard  to  the  general 
economy,  results  from  these  efforts  have  generally  not 
been  very  satisfactory  and  no  widely  accepted  theory  of 
investment  exists  (Ackley  1978,  Clark  1979). 

Aggregate  timber  supply  analyses  have  employed  a 
variety  of  assumptions  and  approaches  for  modeling 
timber  investment  patterns  in  long-run  supply  studies. 
The  importance  of  transitional,  as  opposed  to  steady- 
state,  aspects  of  aggregate  timber  supply  has  been  em- 
phasized to  different  degrees  in  various  studies  of 
timber  supply.  Transitional  patterns  involve  states  of 
disequilibrium  and  change  in  the  stumpage  sector.  The 
disequilibrium  pertains  to  both  stock  variables  that  have 
no  time  dimension  (e.g.,  timber  inventory),  flow  variables 
(e.g.,  harvest)  that  are  time  dependent,  and  associated 
ratios  of  the  two.  Movements  in  the  "system"  compris- 
ing the  timber  economy  may  result  in  an  essentially  con- 
tinual disequilibrium  over  time  because  of  continuous 
changes  in  external  circumstances  such  as  productive 
techniques,  population,  consumer  preferences,  and  gov- 
ernment actions. 

Conversely,  some  studies  emphasize  steady-state 
behavior  for  at  least  some  aspects  of  a  timber  supply 
system  over  time.  Examples  are  long-run  harvest  levels 
formulated  according  to  constant  economic  parameters 
over  time,  or  repetitive  cycles  of  identical  forest 
management  strategies  on  land  assumed  to  be  initially 
bare.  Steady-state  conditions  imply  that  relevant 
variables  all  grow  at  an  identical  rate,  which  is  a 
generalization  of  the  concept  of  the  stationary  state  for 
which  the  relevant  variables  all  remain  constant  (i.e., 
grow  at  a  zero  rate).  The  stationary  equilibrium  state  is 
one  in  which  the  system  operates  in  a  similar  manner 
year  after  year  unless  it  is  disturbed  by  some  outside 
force.  Furthermore,  in  steady-state  analysis  production 
functions  or  technology  are  assumed  to  be  invariant 
over  time.  Some  analyses  may  also  assume  asymptotic 
behavior  in  relation  to  inherent  limits  for  the  system. 

The  theory  of  long-term  timber  investment  based  pri- 
marily on  steady-state  or  static  profit  maximization  prin- 
ciples is  discussed  by  Hyde  (1980)  and  Jackson  (1980). 
Hyde  applies  analytical  tools  based  on  economic  effi- 
ciency norms  in  a  Pacific  Northwest  timber  supply 
analysis,  and  devotes  special  attention  to  public  forest 
land  activities.  As  previously  mentioned,  Hyde  also  gives 
an  extended  theoretical  treatment  of  the  economically 
efficient  allocation  of  forest  land  among  market  and 


nonmarket  uses,  which  encapsulate  long-term  invest- 
ment strategies. 

Hyde's  static  approach  focuses  on  the  economic  effi- 
ciency of  public  land  allocation  and  management,  and 
he  concludes  that  increased  levels  of  both  timber  and 
nontimber  outputs  are  possible.  This  approach  is  similar 
to  that  of  Clawson  and  Hyde  (1976)  and  it  does  not 
address  some  timber  market  ramifications  of  increased 
public  supply.  These  include  the  responses  of  private 
producers  to  any  changes  in  public  supply  behavior  and 
interregional  linkages  of  product  and  stumpage 
markets.  The  important  question  in  the  Pacific  North- 
west with  respect  to  old-growth  timber  and  conversion 
policies  is  handled  in  a  very  general  manner,  with  no 
examination  of  dynamic-adjustment  paths.  Feasibility  of 
attaining  a  long-run  annual  harvest  level  that  would  not 
significantly  increase  relative  stumpage  prices  is  in- 
vestigated, rather  than  the  actual  path  of  an  optimal 
conversion  policy  and  the  resulting  market  implications. 

Jackson  (1980)  orients  his  theoretical  discussion 
toward  the  microeconomic  framework  of  the  private 
timber  firm,  developing  a  theory  of  timber  supply  under 
conditions  in  which  production  is  time  dependent.  For 
example,  under  a  particular  set  of  assumptions,  regen- 
eration inputs  and  rotation  ages  as  decision  variables 
are  shown  to  be  substitutes,  rather  than  complements, 
in  the  productive  relationships.  These  types  of  timber 
production  relations  are  important  for  long-term  timber 
management  analysis,  and  Jackson  also  points  out  infor- 
mation gaps  that  exist  in  the  area  of  timber  production 
processes. 

Two  different  types  of  timber  investment  analysis 
techniques  are  discussed  next  which  address  dynamic 
aspects  of  long-term  timber  supplies— optimal  control 
theory  and  analysis  of  investment  opportunities  in  an 
overall  simulation  framework.  These  two  analytical 
systems  reflect  a  contrast  between  theoretical  and  ap- 
plied timber  investment  modeling.  They  also  focus  on  dif- 
ferent aspects  of  the  multidimensional  timber  supply 
question,  particularly  the  distinction  between  short- 
term  transitional  supply  responses  and  long-run  sta- 
tionary state  behavior. 

OPTIMAL  CONTROL  THEORY 

Optimal  control  models  describe  the  evolvement  of  a 
system  over  time  and  determine  optimal  levels  of  deci- 
sion variables  over  time,  relying  on  the  maximum  princi- 
ple of  optimality  (Clark  1976).  State  equations  describe 
the  evolution  of  the  system  from  an  initial  state, 
resulting  from  the  application  of  a  given  control.  The 
fundamental  problem  in  optimal  control  theory  is  to 
determine  a  feasible  control  that  maximizes  an  objective 
functional  (i.e.,  an  optimal  control)  in  order  to  determine 
an  optimal  path  or  trajectory  through  time.  The  max- 
imum principle  gives  certain  necessary  conditions  that 
must  be  satisfied  by  an  optimal  control.  Techniques  of 
optimal  control  theory  can  handle  nonlinear  optimiza- 
tion problems  and  inequality  constraints. 

SPOC— The  Supply  Potential  Optimal  Control  (SPOC) 
model  is  an  optimal  control  theory  model  for  estimating 
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regional  long-term  supply  of  timber  (Lyon  and  Sedjo 
1983).  A  recent  application  using  hypothetical  timber 
supply  data  addresses  the  constrained  maximization 
problem  of  optimal  timber  harvesting  and  forest  in- 
vestments or  management  practices.  A  long-run  supply 
potential  based  on  economic  considerations  is  investi- 
gated, using  an  endogenous  determination  of  the 
economically  optimal  intertemporal  levels  of  silvicul- 
tural  inputs.  Optimal  rate  of  drawdown  or  conversion  of 
existing  old-growth  stands  is  examined,  and  output 
levels  are  projected  for  steady-state  conditions  after  the 
transition.  Timber  strategies  that  are  examined  pertain 
to  the  rate  at  which  the  timber  stock  is  changed  by  rates 
of  reforestation  and  afforestation,  and  also  by  other 
management  that  alters  biological  growth  rates. 

In  the  SPOC  model,  composition  of  the  long-run  supply 
curve  changes  as  the  transition  progresses  from  short- 
run  fixed  timber  stocks  to  those  in  the  long  run,  where 
the  stock  is  variable.  Once  the  transition  is  complete,  a 
steady-state  forest  situation  ensues.  The  length  of  the 
transition  period  is  specified  beforehand,  and  the  SPOC 
model  does  not  solve  for  it  using  a  simultaneous  optimi- 
zation approach. 

Optimal  long-term  timber  supply  potential  is  analyzed 
by  the  SPOC  model  as  a  discrete  optimal  control  theory 
problem,  with  the  large  multiperiod  problem  decom- 
posed into  a  series  of  smaller,  single  time  period  prob- 
lems. The  model  optimizes  the  level  of  management  input 
and  rotation  length,  and  provides  projections  of 
regeneration,  growth,  and  harvest  for  each  class  of 
forest  lands.  Important  real  world  considerations  such 
as  ownership  patterns  and  institutional  constraints  are 
considered  to  be  noneconomic  factors  and  are  not  ad- 
dressed, given  that  the  intent  of  the  SPOC  modeling  is  to 
focus  on  basic  long-term  biological  and  economic 
processes. 

Biological  variables  are  embodied  in  production  func- 
tions in  the  SPOC  model;  initial  conditions  and  "laws  of 
motion"  comprise  necessary  background  information. 
Initial  conditions  include  such  items  as  acres  of  forest  by 
age  group  and  land  class  and  composite  regeneration  in- 
put for  each  of  the  acres.  Laws  of  motion  are  rules  that 
govern  the  system,  including  those  that  redefine  acres  of 
trees  in  an  age  group  and  the  regeneration  input  from 
one  year  to  the  next. 

Economic  optimization  by  the  SPOC  model  proceeds 
by  land  class,  with  optimal  regeneration  input  and  rota- 
tion period  calculated  for  the  stationary  state  and  then 
for  the  desired  transition  period.  Control  variables  in 
the  transition  period  phase  are  harvest  levels  and  the 
level  of  regeneration  input  in  each  year  by  land  class. 
The  stationary  state  variables  are  the  acres  of  trees  by 
age  and  land  class  and  the  associated  regeneration  in- 
put for  each  of  those  acres.  The  model  maximizes  the 
discounted  present  value  of  the  intertemporal  net  sum  of 
producer  and  consumer  surplus  (net  surplus)  subject  to 
an  appropriate  set  of  technological  and  cost  constraints, 
with  the  optimal  control  theory  algorithm  used  as  the 
solution  technique. 

Optimal  control  variables  are  solved  for  by  the  use  of 
gradient  vector  techniques,  in  which  time  paths  are 


based  on  the  direction  of  maximum  ascent  of  the  present 
value  of  a  "net  surplus  hill."  Total  costs  are  the  sum  of 
regeneration  and  harvesting  costs.  Regeneration  input 
expenditures  depend  on  the  acres  harvested,  with  no 
assumed  regeneration  lag.  Harvesting  involves  a  fixed 
entry  cost  per  acre  and  a  variable  cost  that  depends 
upon  total  volume  harvested. 

Values  for  the  objective  function  represent  present 
value  of  the  net  surplus  stream,  which  is  maximized  sub- 
ject to  a  set  of  constraints.  These  are  constraints  on  the 
values  of  the  control  variables  and  the  laws  of  the  mo- 
tion of  the  system,  e.g.,  the  portion  of  acres  harvested 
has  to  be  nonnegative  and  less  than  or  equal  to  one.  Con- 
strained maximization  of  the  objective  function  pro- 
ceeds by  decomposing  the  overall  problem  into  a  series 
of  subproblems  (i.e.,  maximum  principle).  Each  iteration 
in  the  gradient  method  moves  the  time  path  or  trajectory 
of  the  variables  closer  to  the  optimum. 

The  gradient  vector  which  determines  the  direction  of 
maximum  ascent  of  the  net  surplus  hill  is  the  derivative 
of  the  objective  function  with  respect  to  the  control 
variables.  Thus,  it  is  the  vector  of  shadow  prices  or 
values  of  the  control  variables,  i.e.,  rate  of  change  in  the 
objective  function  per  unit  change  in  a  control  variable. 
These  shadow  prices  indicate  whether  an  increase  in  a 
particular  control  variable  will  increase  the  value  of  the 
objective  function. 

Applications  of  the  SPOC  models  have  used  hypothet- 
ical data  to  project  the  economically  optimum  time  path 
of  prices  and  harvest  volumes.  Variants  of  the  model 
can  be  applied  to  examine  regions  where  data  are  avail- 
able; however,  data  availability  in  the  short  term  may  be 
a  major  constraint  on  actual  application  and  testing  of 
the  embodied  techniques.  Lynn  and  Sedjo  (1983)  propose 
that  several  regions  could  be  incorporated,  if  data  were 
available,  by  introducing  each  region  as  a  separate  land 
class  or  groups  of  land  classes. 

The  potential  usefulness  of  optimal  control  models  in 
aggregate  timber  supply  analyses  is  difficult  to  evaluate 
because  it  is  a  relatively  new  analytical  tool.  Zilberman 
(1982)  suggests  that  optimal  control  analysis  is  a  power- 
ful tool  which  expands  the  range  of  issues  that  agricul- 
tural economists  can  deal  with,  as  well  as  increasing 
their  effectiveness.  Wide-scale  application  of  this 
approach  may  be  impractical  at  present,  however, 
because  of  the  data  requirements  for  solving  dynamic 
optimization  problems.  Optimal  control  models  with 
stochastic  elements  are  much  more  complex  than  deter- 
ministic ones  (e.g.,  Lyon  and  Sedjo  1983),  including  the 
most  advanced  technique,  named  adaptive  control,  that 
takes  into  account  the  expected  gains  from  future  learn- 
ing in  determining  optimal  control  levels  (Zilberman 
1982). 


RELATED  APPROACHES 

Cohan  (1982)  also  utilized  optimal  control  theory  in 
constructing  a  detailed  theoretical  model  of  private  sec- 
tor forest  management  and  timber  supply.  This  long-run 
model  determines  optimal  investments  in  timber  man- 
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agement  activities  (e.g.,  fertilization),  the  timing  and 
levels  of  harvest,  and  cost-minimizing  harvest  patterns. 
This  dynamic  model  of  forest  management  and  timber 
supply  is  based  on  the  assumption  of  profit  maximizing 
behavior  and  is  linked  to  a  set  of  simpler  models  of  mills, 
market  interactions,  and  forest  products  demand.  This 
network  is  designed  to  be  a  flexible  system  of  market 
equilibrium  models  of  the  timber  industry.  This  optimal 
control  model  has  also  not  been  applied  extensively  to 
date. 

Mathematical  programming  techniques,  including 
linear  programming,  goal  programming,  dynamic  pro- 
gramming, integer  and  mixed  integer  programming,  and 
quadratic  programming,  have  also  been  used  in  timber 
investment  analyses  (Bare  1971,  Martin  and  Sendak 
1973).  They  are  designed  to  select  an  optimal  solution 
for  a  set  of  variables,  often  called  activities.  The  optimal 
outcome  is  the  numerical  maximum  or  minimum  of 
some  specified  performance  criterion  or  objective 
function. 

Dynamic  programming  (DP)  has  been  used  at  the 
stand  level  to  solve  for  the  simultaneous  determination 
of  rotation  and  intermediate  treatments  in  timber  invest- 
ment analyses.  The  DP  algorithm  requires  discrete  time 
and  stocking  interval  specification,  in  contrast  to  the  use 
of  continuous  time  or  removal  variables  in  some  optimal 
control  theory  formulations.  DP  and  its  applications  to 
forestry  are  discussed  by  Schreuder  (1968)  and  Brodie 
et  al.  (1978).  One  of  the  earliest  reported  applications 
for  analyzing  multistage,  sequential  decisionmaking 
regarding  timber  investment  paths  was  that  of  Amidon 
and  Akin  (1968).  This  is  essentially  a  recasting  of  Chap- 
pelle  and  Nelson's  (1964)  marginal  analysis  of  optimal 
loblolly  pine  stocking  levels  at  the  stand  level  (based  on 
an  iterative  computer  program)  into  a  dynamic  program- 
ming mode.  More  recent  applications  of  DP  are  dis- 
cussed by  Brodie  et  al.  (1978),  Brodie  and  Kao  (1979), 
Kao  and  Brodie  (1979),  Martin  and  Ek  (1981),  and  Riit- 
ters  et  al.  (1982). 

ANALYSIS  OF  INVESTMENT  OPPORTUNITIES 

TAMM.— The  TAMM  model  (Adams  and  Haynes 
1980)  examined  in  the  previous  section  represents  the 
applied  state  of  the  art  in  aggregate  timber  investment 
modeling.  TAMM  embodies  a  unique  interregional  link- 
age among  long-term  management  or  investment  model- 
ing, short-term  supply  (harvest)  estimation,  and  stump- 
age  price  determination.  Investment  modeling  by  TAMM 
is  based  on  normative  analyses  that  involve  economic 
maximization  calculations,  but  are  later  linked  to 
positive  harvest  flows  modeling.  Application  of  TAMM 
timber-investment  modeling  in  an  aggregate  timber- 
supply  analysis  is  discussed  by  Adams  et  al.  (1982). 

The  TAMM  investment  model  is  an  iterative  process 
that  determines  the  level  of  timber  investment  based  on 
present  net  worth  calculations  and  the  associated  ex- 
tent of  growth  augmentation  in  the  TRAS  timber  inven- 
tory projection  module  (Adams  et  al.  1982).  Stumpage 
price  impacts  of  those  added  investments  are  then 


simulated,  with  only  a  part  of  the  eligible  investments 
undertaken  each  year.  A  new  stumpage  price  is  then  fed 
back  to  recalculate  the  level  of  timber  investment  and 
growth  augmentation.  This  entire  cycle  is  repeated  until 
there  is  a  sufficient  convergence  of  stumpage  price.  Con- 
vergence of  price  implies  that  price  expectations  at  time 
of  investment  are  actually  realized  when  the  investment 
matures  and  the  resulting  stumpage  increment  is 
marketed  (i.e.,  perfect  price  expectations). 

Forest  investments  influence  short-term  supply  rela- 
tions in  the  TAMM  model,  changing  future  harvests  and 
stumpage  prices.  Owners  are  assumed  to  recognize  this 
interrelation  between  investments  and  future  prices,  in- 
cluding all  future  investments,  under  the  assumption  of 
perfect  price  expectations.  Constraints  on  rates  of  in- 
vestments are  imposed  by  limited  availability  of  capital, 
with  economic  gain  from  timber  production  assumed  to 
be  the  sole  objective  of  all  owners. 

Constant  discount  or  interest  rates  are  used  in  the 
timber  investment  analysis  in  TAMM  to  account  for 
opportunity  costs  of  investing  in  timber  growing; 
however,  forestry  is  the  only  major  land  use  considered, 
and  changing  relative  prices  for  products  from  other 
uses  of  the  land  are  not  considered.  Timber  growers  are 
assumed  to  have  perfect  price  expectations  (i.e.,  ex- 
pected prices  are  actually  realized  when  investments 
mature  several  decades  into  the  future)  and  to  have  the 
sole  objective  of  maximizing  their  financial  position. 
Timber  production  yields  and  returns  are  also  usually 
modeled  under  certainty,  with  no  explicit  allowance  for 
land  use  shifts  or  productivity  setbacks. 

Reasonableness  of  assumptions  such  as  economic  ra- 
tionality and  perfect  price  expectations  by  private 
owners  are  difficult  to  evaluate  fully  because  of  the  data 
deficiencies  surrounding  actual  private  timber  invest- 
ment behavior.  Assuming  profit  maximization  as  the  sole 
objective  of  the  selection  and  implementation  of  timber 
management  practices  on  private  lands  ignores  other 
owner  objectives  and  concerns  for  amenity  and  other 
nonmarket  forest  outputs.  The  annual  rate  at  which  ad- 
missible investments  are  undertaken  is  arbitrary 
because  little  is  known  regarding  actual  financial  and 
other  constraints  that  influence  private  timber  invest- 
ment behavior.7 

The  TAMM  investment  analysis  techniques  were  first 
employed  in  the  1980  RPA  Assessment  (USDA  Forest 
Service  1982).  The  development  of  the  SPATS  model  for 
the  South  has  also  expanded  the  analytical  capability  of 
the  TAMM  system  in  terms  of  examining  the  timber  sup- 
ply impacts  from  government  expenditures  to  subsidize 
timber  management  on  nonindustrial  lands  (Brooks 
1983b).  Continuing  research  on  the  TAMM  analysis 
framework  includes  development  of  a  Timber  Resource 
Inventory  Model  (TRIM)  for  major  timber  supply  regions 
that  incorporates  modeling  of  management  intensifica- 
tion practices.  The  general  format  of  the  TRIM  model  is 
based  on  attributes  of  the  aggregation  structure  of  the 
TREES  model  (Tedder  et  al.  1980)  discussed  earlier  and 
also  on  TAMM's  investment  modeling  (Tedder  1983). 

Modeling  of  management  intensification  for  even- 
aged  stands  by  the  TRIM  model  relies  on  the  yield  table 
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projection  approach.  This  age  class  approach,  similar  to 
that  for  the  TREES  model,  allows  the  timber  inventory  to 
be  segregated  so  that  management  practices  can  be  tai- 
lored according  to  site,  species,  stocking  levels,  and  age. 
This  is  more  specific  than  the  1980  TAMM  approach, 
which  relied  on  an  augmented  TRAS  model  (Barber 
1980),  where  all  acres  in  an  aggregate  were  assumed  to 
receive  the  same  basic  package  of  management  prac- 
tices over  time. 

The  economic  process  that  drives  the  TRIM  invest- 
ment model  is  based  on  perfect  price  expectations,  as  in 
the  1980  TAMM  approach.  If  the  estimated  rate  of 
return  is  equal  to  or  exceeds  a  specified  alternative  rate 
of  return,  investments  are  implemented  on  a  perpetual 
rotation  basis.  Annual  prices  used  in  the  calculations  of 
soil  expectation  values  are  taken  from  the  TAMM 
model,  while  costs  will  be  obtained  from  regional  studies 
(e.g.,  Rahm  1981  and  Dutrow  et  al.  1982).  This  timber 
resources  inventory  projection  model  would  make  pro- 
jections by  ten-year  intervals,  while  linked  to  the  annual 
TAMM  market  model. 

OTHER  MODELING  CONSIDERATIONS 

Applications  of  techniques  for  analyzing  timber  in- 
vestment strategies  are  numerous,  and  many  variations 
have  been  applied  to  the  broad  range  of  timber  species 
across  the  nation.  A  cross  section  of  these  timber  invest- 
ment analysis  techniques  are  listed  by  region  in  appen- 
dix C,  which  also  includes  sources  of  cost  information 
useful  for  investment  analyses.  These  techniques  have 
frequently  been  constructed  to  analyze  stand  level  treat- 
ment possibilities,  and  they  are  typically  designed  for 
analyses  at  substate  levels  (e.g.,  Hepp  1982).  Thus, 
analysis  of  the  biological  and  economic  characteristics 
of  alternative  timber  management  schemes  at  aggregate 
levels  would  probably  require  considerable  adjustments 
in  associated  timber  growth  and  yield  modeling  to  better 
address  the  less  detailed  but  more  extensive  data  needs 
of  aggregate  analyses. 

All  of  the  major  timber  investment  analysis  tech- 
niques currently  available  are  based  on  normative  deci- 
sion rules,  which  assume  that  owner  objectives  are 
generally  based  on  maximization  of  economic  gain  from 
the  timber  resource  alone.  No  positive  aggregate  timber 
investment  model  has  been  reported,  and  timber  supply 
analysts  face  the  question  of  how  to  represent  ade- 
quately realistic  future  price  expectations  mechanisms 
and  likely  supply  responses  by  stumpage  producers. 
Because  planned  supply  responses  or  investment  deci- 
sions involve  estimates  of  the  future  based  on  imperfect 
knowledge  regarding  stumpage  prices,  technology, 
yields,  institutions,  etc.,  those  future  estimates  are  sub- 
ject to  error.  Those  errors  give  rise  to  differences  be- 
tween planned  and  realized  supply  response  (Jensen  and 
Day  1961,  Stone  1970),  and  little  is  known  regarding  the 
extent  of  this  gap  for  timber  investments.  Most  timber 
investment  analyses  utilized  in  examining  future  timber 
supply  levels  also  usually  assume  that  all  timber  invest- 
ments actually  implemented  will  be  carried  to  maturity 
at  a  certain  productive  level. 


This  latter  assumption  may  be  too  optimistic  in  some 
cases  because  of  evidence  of  investment  failures  and 
land  use  changes  cited  in  several  studies.4  Studies  of  the 
retention  and  condition  of  timber  investments  initiated 
through  public  programs  are  listed  in  appendix  D.  In 
general,  an  inconsistent  and  rather  meager  information 
base  exists  to  support  the  appropriate  degrees  of  adjust- 
ment of  physical  and  financial  yields  of  timber  invest- 
ments to  account  for  unforeseen  productivity  reductions 
and  land  use  changes.  This  is  particularly  true  for  non- 
industrial,  private  timber  investments,  which  often  are 
managed  by  several  different  landowners  over  a  matu- 
ration period  that  typically  spans  at  least  several 
decades. 

The  various  approaches  to  modeling  aggregate  timber 
investment  have  generally  been  normative  in  nature. 
They  have  emphasized  different  biological  and  economic 
dimensions,  especially  as  they  relate  to  the  dynamics  of 
timber  supply.  Static  optimization  techniques  have  been 
utilized  to  investigate  displacement  from  equilibrium 
positions  based  on  competitive  market  conditions,  speci- 
fied constraints,  and  the  objective  of  profit  maximiza- 
tion, but  the  explicit  time  path  of  variable  changes  is  not 
predicted.  A  hybrid  framework  that  uses  short-run,  em- 
pirically derived  harvest  relations  in  conjunction  with 
normative  long-term  timber  investment  modeling  in  a 
dynamic  analysis  supports  the  first-generation  TAMM 
model. 

The  proliferation  in  the  application  of  operations 
research  techniques  in  forestry  over  the  past  several 
decades  has  also  led  to  a  large  body  of  literature  on 
timber  investment  analyses  using  mathematical  pro- 
gramming models  suited  to  electronic  data  processing 
(e.g.,  Brodie  et  al.  1978).  The  mathematical  basis  of  op- 
timization techniques  necessarily  implies  abstract 
representations  of  biological  and  social  systems,  with 
approximations  and  simplifying  assumptions  generally 
required  if  the  model  is  to  be  tractable.  The  proper 
criterion  for  judging  whether  the  fomulated  model  is 
then  a  valid  representation  is  whether  it  predicts  the 
relative  effects  of  alternative  courses  of  action  with  suf- 
ficient accuracy  to  permit  a  sound  decision  (Hillier  and 
Lieberman  1980).  Such  determination  is  quite  subjective 
with  respect  to  timber  investment  because  long  matura- 
tion periods  hinder  timely  objective  validation  (and  con- 
sequent modification,  if  needed). 


AGGREGATION  AND  UNCERTAINTY 
CONSIDERATIONS 

Estimating  the  validity  of  aggregate  models  and  the 
concurrent  degree  of  uncertainty  associated  with  aggre- 
gate timber  supply  studies  is  a  major  problem.  The 
degree  of  uncertainty  inherent  in  estimating  future 
components  of  timber  supply  is  strongly  related  to  the 
aggregation  schemes  employed,  which  by  their  nature 
condense  diverse  and  complex  relationships  into  a 
relatively  small  number  of  essential  characteristics.  The 
broad  geographical  range  of  most  aggregate  analyses 
invariably  includes  diverse  owner  classes  with  differing 
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timber  management  motivations,  institutional  con- 
straints, resource  and  market  knowledge,  etc.  Not  only 
does  this  imply  that  predicting  future  timber  supply  is 
difficult,  but  in  many  situations  substantial  uncertainty 
exists  with  regard  to  current  areas,  volumes,  ecological 
capabilities,  prices,  etc.  (Bentley  1981). 

Aggregation  and  uncertainty  in  timber  supply  model- 
ing are  related  in  part  to  the  quantity  and  quality  of 
historical  data  available  as  a  basis  for  projection.  Data 
deficiencies  vary  with  respect  to  biological  and  eco- 
nomic dimensions  of  timber  supply  modeling,  and  also 
within  these  two  broad  categories.  The  levels  of  timber 
growth  and  yield  data  have  been  augmented  signifi- 
cantly during  the  last  several  decades;  substantial  prob- 
lems remain,  however,  regarding  the  suitability  of  such 
data  for  purposes  of  modeling  responses  to  an  array  of 
management  treatments. 

Availability  and  quality  of  data  for  economic  charac- 
teristics of  timber  management  is  also  highly  variable; 
data  are  particularly  deficient  in  certain  areas.  A  com- 
prehensive data  base  for  modeling  landowner  response 
to  alternative  market  conditions  does  not  exist;  there- 
fore, consideration  of  the  dynamic  nature  of  human  in- 
teractions with  biological  and  economic  systems  in  the 
real  world  is  quite  limited.  These  deficiencies  have 
forced  analysts  to  employ  less  data-demanding  models, 
often  limiting  the  expansion  of  analyses  to  adequately 
include  other  resource  interactions  and  impeding  the 
construction  of  timber  supply  analyses  at  less  aggre- 
gated levels.  Strategies  for  dealing  with  data  limitations 
are  offered  by  Rose  et  al.  (1981),  while  Hamilton  (1978) 
and  Lewis  and  Ellefson  (1983)  discuss  related  aspects  of 
required  levels  of  precision  with  respect  to  information 
needs  for  forest  management. 

Even  with  ideal  data  sources,  however,  determination 
of  the  mode  of  data  aggregation  for  regional  and  na- 
tional timber  supply  studies  would  likely  not  involve  a 
straightforward  decision  process.  Appropriate  geo- 
graphic bases  sometimes  vary  greatly  for  analyzing  eco- 
logical, social,  and  economic  concerns  in  supply  studies 
that  include  multiresource  and  institutional  considera- 
tions. Although  aggregation  patterns  may  have  a  sub- 
stantial bearing  on  the  quality  of  supply  analysis  results, 
a  set  of  useful  aggregation  guidelines  based  on  adequate 
quantitative  considerations  has  not  emerged.  At  least 
initially,  aggregation  schemes  are  dictated  by  the  objec- 
tives of  a  given  study,  but  the  direction  of  later  adjust- 
ments in  response  to  data  availability  is  often  not 
predictable. 

Teeguarden  (1977),  in  discussing  the  complex  of  fac- 
tors involved  in  aggregation  considerations  for  RPA 
Assessments,  suggests  that  relatively  homogeneous  eco- 
systems be  the  lowest  level  of  analysis,  with  subsequent 
aggregation  to  regions  and  to  the  nation.  Marty  (1969) 
and  Beuter  (1979)  have  suggested  the  need  for  less  ag- 
gregated geographic  and  ownership  groupings  than 
have  been  utilized  in  national  timber  appraisals.  While 
other  suggestions  can  be  found  in  the  literature,  no  com- 
plete set  of  readily  transferable  guidelines  based  on 
quantitative  or  economic  principles  is  available  to  aid  in 
devising  aggregation  schemes  for  timber  supply  studies. 


The  question  of  appropriate  level  of  aggregation  also 
needs  to  be  examined  regarding  timber  demand  as  well. 

Two  additional  aspects  of  data  aggregation  that  may 
cause  problems  in  timber  supply  modeling  pertain  to 
nonlinear  systems  and  variance  of  aggregates.9  In  the 
first  instance,  aggregation  can  cause  bias  in  systems 
that  are  nonlinear,  particularly  for  systems  with  many 
step  functions  (e.g.,  decisions  to  harvest,  depending  on 
volume  or  costs).  The  second  case  pertains  to  whether  a 
particular  model  adequately  represents  the  effects  of 
variance  in  aggregates.  For  systems  with  strong  positive 
feedback,  ignoring  variance  may  introduce  bias  (e.g.,  ig- 
noring the  variability  of  growth  rates  from  the  mean  for 
large,  dominant  trees  may  produce  underestimation  of 
growth). 

Aggregate  timber  supply  results  are  generally  based 
on  analytical  techniques  dependent  upon  a  mixture  of 
experimental  and  nonexperimental  data  sources,  with 
little  (if  any)  indication  of  the  overall  statistical  reliabili- 
ty of  the  ultimate  results.  The  manner  in  which  such 
statistical  problems  might  interact  or  even  be  com- 
pounded over  relatively  long  timber  supply  projection 
periods  has  generally  not  been  addressed.  While  the 
problem  of  combining  techniques  or  data  of  varying 
statistical  reliability  is  not  unique  to  timber  supply 
analyses,  the  particular  mixture  of  specific  data  prob- 
lems and  long  time  periods  involved  in  timber  supply 
studies  create  some  special  statistical  concerns.  Propa- 
gation of  error  does  appear  to  limit  the  usefulness  of 
large,  relatively  complex  models. 

One  general  approach  used  to  address  these  prob- 
lems has  been  the  use  of  sensitivity  analyses,  the  testing 
of  the  influence  of  changes  in  input  items  or  assumptions 
on  model  results.  Sensitivity  analyses  are  especially 
relevant  for  mathematical  programming  techniques  that 
are  widely  used  in  forest  land  planning,  given  the  long 
planning  horizons  and  the  assumption  of  certainty  for 
model  imputs. 

Sensitivity  analyses  have  been  applied  in  various 
fashions  for  timber  supply  analyses.  For  example, 
Frayer  and  Jones  (1970)  use  a  Monte  Carlo  technique  to 
investigate  the  effects  of  each  input  item  on  timber 
stand  projection  variability.  Results  indicate  that  sam- 
pling errors  for  growth,  mortality,  removal,  and  in- 
growth parameters  largely  determine  projection  vari- 
ability. Goforth  and  Mills  (1975)  devised  a  technique  to 
determine  degree  of  data  error  required  to  influence  the 
outcome  of  timber  investment  analyses. 

One  proposal  for  improving  credibility  of  natural 
resource  supply  forecasts  involves  the  comparison  of  a 
series  of  relatively  short-term  projections  with  one  con- 
tinuous projection.5  In  this  approach,  reexamination 
(and  readjustment  if  necessary)  of  assumptions  and 
techniques  are  possible  at  each  breakpbint.  Another 
suggestion  involves  definition  of  a  systematic  procedure 
to  ensure  that  each  component  of  the  projection  process 
and  the  overall  process  itself  is  examined  for  validity  by 
specialists  other  than  the  forecaster/analyst  and  the 
potential  users.5 

Reducing  the  degree  of  uncertainty  in  timber  supply 
analyses,  or  at  least  not  unduly  obscuring  it  in  describ- 


26 


ing  results  of  such  studies,  is  a  primary  consideration  in 
devising  aggregation  schemes.  Much  has  been  written 
on  the  general  topic  of  uncertainty  associated  with 
investment  analyses  and  planning  (e.g.,  Borch  1968, 
Arrow  and  Lind  1970),  with  applications  to  forest  land 
use  planning  discussed  by  Dowdle  (1962),  Marty  (1964), 
Thompson  (1968),  Teeguarden  (1969),  Lundgren  and 
Thompson  (1972),  and  Fight  and  Bell  (1977),  among 
others.  As  noted  in  the  previous  section,  Hoganson 
(1981)  offers  a  specific  approach  for  addressing  the 
problem  of  uncertainty  and  accountability  in  stand-level 
optimization  (i.e.,  linear  programming)  modeling  by 
reformulating  the  linear  programming  problem  within  a 
simulation  framework.  However,  there  appear  to  be  no 
standard  pathways  to  follow  in  evaluating  and  tem- 
pering aggregate  forest  investment  analysis  or  supply 
projections  with  regard  to  uncertainty.  Zones  of  uncer- 
tainty, usually  expanding  greatly  with  time,  are  some- 
times assigned  to  forecasted  estimates  in  timber  supply 
studies.  However,  their  estimation  and  use  appear  to  be 
rather  subjective  in  nature. 

A  relatively  large  degree  of  uncertainty  in  timber  sup- 
ply analyses  is  related  to  the  estimation  of  future 
technological  changes.  This  refers  to  both  innovation 
and  implementation  of  technologies  across  diverse 
owner  settings.  According  to  economic  theory,  produc- 
tion function  frontiers  are  based  on  application  of  the 
most  efficient  technology,  which  therefore  implies  a 
fixed  technology  base  over  the  production  period;  how- 
ever, production  processes  that  stretch  over  several 
decades  at  a  minimum  suggest  that  this  assumption  may 
not  be  valid.  While  many  recognize  that  substitution  of 
other  inputs  (e.g.,  genetically  improved  trees)  for  land  in 
timber  production  will  increase  in  the  future,  precise 
forecasts  of  the  rate  and  actual  composition  of  the 
technical  change  are  difficult  to  obtain.  Heady  (1952) 
has  outlined  some  basic  principles  to  consider  in 
attempting  to  approximate  the  direction  of  the  impact  of 
technical  change;  these  are  based  on  the  nature  of  the 
technical  change  and  demand-supply  characteristics  of 
products  (Castle  1977). 

Aggregate  timber  supply  studies  can  provide  critical 
clues  to  issues  regarding  natural  resources  availability, 
but  often  such  studies  are  questioned  because  adequate 
measures  of  risk  or  uncertainty  are  not  provided.  Most 
aggregate  studies  have  employed  deterministic  models 
that  provide  only  a  point  estimate  of  the  outcome  for  a 
specified  set  of  assumptions  and  timber  management 
strategies.  This  is  again  related  in  part  to  data  prob- 
lems. It  is  not  uncommon  for  timber  growth  and  yield 
projections  to  be  published  without  any  meaningful  in- 
dication of  changes  in  the  level  of  reliability  over  the 
projection  period. 

In  summary,  given  that  substantial  uncertainty  and 
aggregation  bias  might  be  associated  with  aggregate 
timber  supply  studies,  the  degree  of  sophistication  war- 
ranted in  a  study's  framework  and  underlying  modeling 
techniques  is  largely  a  judgmental  matter.  Standard 
statistical  tests  or  other  objective  guidelines  are  very 
difficult  to  apply  in  a  comprehensive  manner,  largely 
because  of  the  quality  and  mixture  of  available  data.  A 


key  question  in  the  planning  stages  of  a  timber  supply 
study  is  whether  the  study  approach  is  consistent  and 
efficient  in  terms  of  the  stated  objective.6  However,  cur- 
rently there  are  difficulties  in  determining  changes  in 
the  quality  of  projection  results  needed  to  improve  deci- 
sionmaking to  selected  degrees.  This  has  also  impeded 
rigorous  and  objective  evaluations  of  model  adequacy. 

OVERVIEW  OF  METHODOLOGY 

Techniques  and  supporting  data  bases  for  analyzing 
aggregate  timber  supplies  have  evolved  at  different 
rates  for  the  four  major  analytical  components  re- 
viewed: forest  land  use  allocation,  growth  and  yield  of 
the  timber  resource,  removals  or  harvest,  and  timber  in- 
vestments over  time.  Aggregate  timber  supply  analyses 
prior  to  1970  typically  depicted  timber  supplies  in 
physical  or  biological  dimensions,  with  relatively  little 
analysis  of  their  relationship  to  market  forces.  Applied 
analytical  techniques  have  generally  evolved,  however, 
from  relatively  simple  ratio  and  trend  analyses  to  the 
use  of  complex  econometric  techniques  in  modeling  the 
equilibrium  of  supply  and  demand. 

The  allocation  of  land  to  timber  production  has  been 
exogenously  determined  in  national  timber  assessments, 
with  adjustments  in  a  residual  timberland  base  reflect- 
ing changes  in  higher  valued  land  uses.  Systematic 
methodologies  for  projecting  timberland  acreages  have 
been  based  largely  on  economic  optimization  criteria 
(e.g.,  Hyde  1980),  with  notable  exceptions  being  Markov 
land  use  simulation  (e.g.,  Burnham  1973)  and  econo- 
metric approaches  (e.g.,  White  and  Fleming  1980). 

The  accuracy  of  results  from  projecting  timberland 
acreages  with  techniques  that  embody  optimization 
criteria  has  not  been  satisfactory  in  some  cases, 
especially  with  regard  to  nonindustrial,  private  lands. 
Simulation  techniques  applied  in  attempting  to  remedy 
this  situation  have  also  met  with  limited  success, 
primarily  because  of  land  use  data  limitations. 
Documented  information  for  historical  timberland 
acreages  is  characterized  by  irregular  time  frames  and 
nonstandard  classification  criteria  among  major  data 
sources. 

Assuming  that  a  method  for  characterizing  the  forest 
resource  has  been  selected  and  the  resource  has  been 
stratified  on  the  basis  of  such  a  characterization,  timber 
inventory  projection  models  are  used  to  project  repre- 
sentative units  of  each  stratum  (or  the  entire  stratum) 
over  time  in  timber  supply  analyses.  For  example,  simul- 
taneous growth  and  yield  equations  (e.g.,  Sullivan  and 
Clutter  1972)  can  be  utilized  for  stand-level  projections 
for  a  given  initial  density,  age,  site,  etc.,  stratum.  In  an 
alternative  approach,  the  STEMS  model  (USDA  Forest 
Service  1979)  generates  projections  on  an  individual- 
tree  basis  for  a  particular  stratum. 

Given  the  diversity  of  species  types  and  physical  set- 
tings, timber  inventory  modeling  has  evolved  according 
to  two  approaches:  direct  and  indirect  methods.  Direct 
methods  involve  modeling  growth  and  yield  for  forest 
areas  using  only  sampling  information  for  those  areas 
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(i.e.,  applying  sample  data  directly  back  to  the  popula- 
tion from  which  the  sample  was  obtained).  An  example 
is  a  stand  table  projection  technique  such  as  the  TRAS 
model,  which  applies  rates  of  change  (e.g.,  mortality 
rate)  to  a  set  of  initial  stand  conditions  (Alig  et  al.  1982). 
TRAS  has  been  used  in  several  regional  and  national 
timber  supply  appraisals,  but  its  design  makes  it  some- 
what inflexible  for  modeling  impacts  of  alternative 
timber  management  strategies. 

Indirect  timber  inventory  projection  methods  are  ap- 
plied to  subject  stands  larger  than  those  from  which  the 
original  growth  and  yield  data  are  obtained,  and  essen- 
tially involve  a  two-phase  process.  Whole  stand  yield 
tables  (McArdle  et  al.  1961)  and  variable  density  growth 
and  yield  tables  or  models  (Buckman  1962)  comprise  the 
two  broad  classes  of  indirect  methods.  Variable  density 
growth  and  yield  models  may  be  used  to  generate  yield 
tables,  and  there  are  three  major  classes  of  these 
variable  density  models:  whole  stand  (Buckman  1962), 
size  class  and/or  growth  component  (Moser  1974),  and 
individual-tree  simulation  models  (Stage  1973).  The  first 
two  classes  rely  upon  aggregate  stand  characteristics 
such  as  number  of  trees  or  basal  area  in  modeling  forest 
development.  The  more  recently  developed  individual- 
tree  simulation  models  describe  stands  on  the  basis  of 
characteristics  of  individual  trees. 

Individual-tree  models  are  classified  according  to 
whether  or  not  they  incorporate  intertree  distance 
within  the  projection  framework.  Intertree  distance- 
independent  models  are  represented  by  the  STEMS 
model  (USDA  Forest  Service  1979)  and  the  Prognosis 
model  (Stage  1973).  These  models  do  not  consider  spatial 
patterns  or  individual-tree  locations.  Consequently,  they 
are  operationally  faster  than  distance-dependent 
models  (e.g.,  Daniels  and  Burkhart  1975),  presumably 
without  a  significant  trade-off  in  terms  of  modeling 
precision  for  aggregate  timber  supply  studies.  Distance- 
dependent  models  can  provide  relatively  detailed 
evaluations  of  timber  management  alternatives,  but  are 
relatively  taxing  in  terms  of  data  requirements  and  com- 
puter time. 

The  design  and  implementation  of  an  inventory  pro- 
jection system  that  would  improve  upon  current  capabil- 
ities for  aggregate  timber  supply  analyses  are  buffeted 
by  multifaceted  considerations.  Data  availability, 
accuracy  requirements,  cost  effectiveness,  time  and 
personnel  constraints,  and  manageability  are  examples 
of  interrelated  (and  often  conflicting)  considerations 
that  are  difficult  to  weigh  because  of  quantification 
problems.  While  the  growing  sophistication  of  computer 
programming  technology  facilitates  construction  of  in- 
dividual models  with  larger  storage  and  computation 
capabilities,  the  reliance  upon  a  single  model  (e.g., 
TRAS)  in  future  aggregate  timber  supply  modeling  may 
diminish  as  the  state  of  the  art  evolves.  It  is  unlikely  that 
a  single  growth  model  could  provide  the  most  efficient 
growth  estimates  for  all  species  and  conditions.  Also, 
the  incorporation  or  adaptation  of  all  of  the  component 
relationships  needed  to  apply  any  one  model  to  all 
species  and  conditions  would  likely  be  impeded  by  trans- 
lation and  data  availability  problems.  Linkage  of  models, 


or  the  aggregation  of  results  from  distinct  models  (using 
statistical  reliability  standards),  appears  to  offer  sub- 
stantial promise,  though  coordination  problems  will  like- 
ly be  substantial. 

Timber  inventory  projection  models  suitable  for  large- 
scale  analyses  have  generally  been  developed  sepa- 
rately from  economic  analysis  techniques,  the  latter 
having  been  subsequently  linked  to  them  for  the  purpose 
of  modeling  the  interrelationships  between  biological 
development  of  timberland  and  human  intervention. 
This  has  somewhat  hindered  supply  analysts  from  ex- 
amining impacts  arising  from  shifts  in  timber  manage- 
ment intensity  suggested  by  investment  analyses. 

Harvest  flow  models  interface  with  timber  inventory 
projection  models  in  modeling  the  impacts  on  inventory 
development  that  result  from  adjusting  timber  inventory 
levels  in  response  to  harvest  decisions.  Harvest  flows 
are  often  calculated  as  part  of  an  overall,  optimal, 
timber  management  package  based  upon  long-run  con- 
siderations (e.g.,  Hyde  1980).  However,  in  the  1980  RPA 
Assessment,  such  harvest  flows  were  modeled  as  short- 
run  functions  of  stumpage  price  and  timber  inventory 
levels  (Adams  and  Haynes  1980). 

Simulation  and  optimization  techniques  have  been 
utilized  separately  or  in  various  mixtures  in  modeling 
aggregate  harvest  flows.  These  approaches  differ  with 
respect  to  the  assumed  basis  of  stumpage  supply 
responses  of  producers,  including  the  posited  intertem- 
poral linkage  between  harvests  and  long-range  in- 
vestments. The  TAMM  model  (Adams  and  Haynes  1980) 
simulates  future  harvest  flows  by  using  econometrically 
estimated  short-run  equations  based  on  observed  values 
from  a  historical  period  approximately  several  decades 
in  length.  An  optimization  model  such  as  FORPLAN 
(Johnson  et  al.  1980)  seeks  to  maximize  (or  minimize)  the 
value  of  an  objective  function  that  includes  values  for 
decision  variables  representing  management  practices 
for  specific  areas.  Thus,  simulation  models  have  been 
employed  to  project  timber  supply  trends  based  on 
historical  patterns  or  other  assumptions,  while  optimiza- 
tion models  can  estimate  potential  supply  responses 
linked  to  efficient  shifts  in  timber  management  on  some 
lands  (assuming  certain  goals).  The  basic  distinction  is 
between  simulating  how  the  system  tends  to  operate,  in 
contrast  to  how  it  could  operate  according  to  some  norm 
or  optimality  rule. 

Simulations  of  timber  supply  activities  in  aggregate 
timber  supply  studies  have  typically  been  made  with 
econometric  techniques,  which  use  statistical  methods 
to  quantify  economic  relationships.  Econometric  models 
for  aggregate  timber  supply  processes  often  include 
several  variables  that  are  exogenous  to  the  model,  thus 
leading  to  a  partially  open  modeling  structure.  In  par- 
ticular, macroeconomic  variables  (e.g.,  interest  rates) 
may  exert  large  influences  on  timber  supply  activities. 
These  are  quite  difficult  to  forecast  for  future  condi- 
tions, particularly  when  the  long  time  horizons  charac- 
teristic of  forestry  planning  are  considered.  Other 
simulation  techniques  that  might  be  used  in  place  of  or 
in  conjunction  with  econometric  techniques  are  sto- 
chastic simulation  and  system  dynamics.  In  particular, 
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system  dynamics  may  allow  a  more  explicit  method  for 
modeling  the  complex  interactions  that  characterize 
timber  supply  processes,  while  providing  a  general  level 
of  understanding  regarding  the  associated  dynamic 
tendencies  over  the  long  time  periods  (Meadows  1980). 

While  simulation  models  are  readily  adaptable  to  the 
modeling  of  aggregate  harvest  flows,  optimization 
techniques  such  as  linear  programming  encounter  dif- 
ficulties regarding  specification  of  appropriate  objec- 
tive functions  and  constraint  equations  for  regional  or 
national  analyses  (McCarl  and  Spreen  1980).  This  prob- 
lem is  significant  for  both  private  and  public  owner- 
ships. The  extent  of  deviation  of  actual  harvesting  pat- 
terns from  "optimal"  timber  harvest  and  management 
strategies  assumed  or  portrayed  in  optimization  analy- 
ses is  presently  unknown,  but  some  evidence  exists 
concerning  management  activities  applied  to  cost- 
shared  timber  investments  on  nonindustrial  private 
lands  (e.g.,  Kurtz  et  al.  1980). 

Modeling  of  aggregate  timber  investment  behavior 
has  been  integrated  with  that  for  harvest  flows  in 
several  different  ways.  An  example  of  incorporating 
timber  investment  patterns  based  on  optimization 
criteria  within  an  overall  timber  supply  simulation  ap- 
proach is  the  TAMM  modeling  applied  for  the  1980  RPA 
Assessment  (Adams  et  al.  1982).  A  different  approach  is 
exemplified  by  the  long-term  timber  supply  potential 
model  developed  by  Lyon  and  Sedjo  (1983).  This  is  based 
on  optimal  control  theory  and  projects  harvest  and  in- 
vestment levels  during  both  transitional  and  steady- 
state  phases.  Such  an  approach  allows  simultaneous 
determination  of  harvest  and  investment  levels  over  an 
analysis  period;  it  has  not  been  extensively  applied  and 
tested  to  date. 

The  current  array  of  timber  investment  analysis  tech- 
niques at  the  stand  level  varies  considerably  across  the 
country,  primarily  because  the  degree  of  sophistication 
of  the  underlying  growth  and  yield  modeling  is  not 
uniform  across  regions.  Timber  growth  and  yield  char- 
acteristics differ  so  greatly  across  some  geographical 
locations  in  terms  of  species,  physiography,  climate,  and 
other  factors,  that  separate  growth  and  yield  models 
must  be  linked  to  relatively  standardized  economic 
analysis  tools  for  particular  regions.  As  a  consequence, 
differences  in  the  evolution  of  data  collection  and  biolog- 
ical response  modeling  efforts  for  various  packages  of 
timber  management  treatments  have  resulted  in  varying 
capabilities  for  evaluating  timber  investment  oppor- 
tunities in  different  parts  of  the  country. 

Specific  areas  of  improvement  needed  in  long-term 
timber  management  modeling  include  positive  timber 
investment  modeling  and  more  in-depth  examination  of 
the  dynamics  of  timber  inventory  development  in  re- 
sponse to  varying  degrees  and  timing  of  investments.  At 
present,  no  positive  models  of  aggregate  timber  invest- 
ment behavior  exist.  Empirical  work  devoted  to  the  sta- 
tistical testing  of  timber  investment  theories  has  been 
quite  limited,  and  no  fully  satisfactory  or  widely  ac- 
cepted theory  of  aggregate,  long-term  timber  investment 
has  been  established.  A  better  understanding  of  private 
timber  investment  decisionmaking  is  required,  par- 


ticularly with  respect  to  the  processes  of  economic 
expectation  formation  (e.g.,  as  for  stumpage  prices)  for 
conditions  that  may  not  actually  materialize  for  several 
decades  at  a  minimum.  Appropriate  levels  of  detail  and 
abstraction  for  aggregate  timber  investment  models  also 
require  careful  scrutiny,  including  the  establishment 
of  useful  measures  of  general  aggregate  investment 
behavior. 

Investment  modeling  based  on  the  assumption  of 
timber  profit  maximization  by  owners  may  represent  an 
upper  limit  or  overestimation  of  actual  investment,  in 
light  of  the  currently  deficient  data  and  analyses  related 
to  other  objectives  for  various  owner  groups.  From  an 
overall  perspective,  generalizations  regarding  appro- 
priate levels  of  abstraction  and  aggregation  in  aggre- 
gate timber  supply  studies  are  tenuous;  this  is  due  in 
part  to  the  varying  objectives  of  such  studies  and  to  diffi- 
culties in  verifying  whether  harvest  and  investment  pat- 
terns depicted  by  these  models  adequately  reflect  actual 
occurrences. 

The  degree  and  timing  of  future  improvements  in  ag- 
gregate timber  supply  analysis  may  depend  upon  meas- 
ures taken  to  strengthen  existing  data  sources.  Data 
deficiencies  have  been  alluded  to  within  numerous 
timber  supply  studies,  which  themselves  have  generally 
drawn  on  a  variety  of  data  sources.  With  respect  to 
basic  timber  supply  data,  timber  resource  data  is  some- 
times criticized  as  being  inadequate  for  modeling  pur- 
poses; however,  economic  data  may  be  even  less  ade- 
quate. This  has  particularly  hampered  econometric 
analyses,  which  require  reasonably  long  series  of 
consistent  price  and  production  data.  For  example,  data 
on  timber  production  and  transportation  costs  based  on 
location  and  accessibility  considerations  have  not  been 
assembled  or  organized  in  readily  retrievable  forms  for 
use  in  extending  supply  modeling  capabilities. 

The  quality  dimensions  of  timber  supply  represent 
another  instance  in  which  data  deficiencies  have  im- 
peded modeling.  Timber  supply  at  aggregate  levels  is 
comprised  of  many  species,  sizes,  and  grades  of  trees  or 
logs,  each  of  which  may  have  a  distinct  price.  Existing 
data  do  not  presently  allow  these  differences  to  be  ex- 
tracted or  monitored. 

Some  of  the  primary  improvements  needed  in  aggre- 
gate timber  supply  data  include: 

1.  Levels  of  private  timber  harvest  by  region  and 
owner  group. 

2.  Timber  inventory  by  region  and  owner  group. 

3.  Socioeconomic  characteristics  of  forest  owners 
(e.g.,  income) 

4.  Actual   timber   investment   behavior  by  forest 
owners. 

5.  Forest  land  management  opportunities  by  region 
and  owner  group  (collected  on  a  periodic  basis).8 

The  bulk  of  the  empirical  aggregate  timber  supply 
analyses  can  be  described  without  reference  to  any- 
thing but  the  most  basic  microeconomic  theory.  While 
the  transition  from  theoretical  abstraction  to  empirical 
reality  has  been  hindered  by  data  problems,  numerous 
opportunities  exist  to  build  upon  basic  timber  supply 
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theory  a  structure  that  is  more  flexible  and  realistic 
with  regard  to  the  diverse  ownership  settings. 

Ultimately,  research  may  provide  an  overall  model 
capable  of  representing  aggregate  timber  supply 
according  to  the  theory  of  the  rational  forest  owner.8 
This  will  entail  the  simultaneous  determination  of  the 
levels  of  timber  harvest,  intensity  of  timber  manage- 
ment, and  land  area  devoted  to  the  timber  growing 
enterprise  over  time  (along  with  supporting  timber 
inventory  projection  capabilities).  It  will  require  more 
attention  to  the  systems  structure  of  the  analytical  prob- 
lem (e.g.,  feedback  loops  and  performance  criteria) 
which  in  turn  may  provide  a  framework  for  integrating 
the  components  of  aggregate  timber  supply  analysis.  An 
additional  requirement  is  a  better  understanding  of 
owner  objectives,  including  the  nature  and  extent  of  in- 
terest in  market  and  nonmarket  forest  products  and  the 
"role"  of  the  forest  enterprise  in  the  overall  operations 
of  a  forest  owner,  as  well  as  the  role  of  forestry  in  local, 
regional,  and  national  economies.  However,  given  the 
state  of  the  art  and  data  deficiencies,  analysts  will  be 
forced  in  the  interim  to  separate  the  overall  problem 
into  smaller,  more  tractable  components,  as  outlined  in 
this  report. 
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Appendix  A.— The  date,  primary  analytical  approach,  and  scope  for 
major  aggregate  timber  supply  studies  in  the  United  States 


Study 


Year 


Analytical  approach 


Scope 


USDA  Forest  Service 
(Capper  Report) 

USDA  Forest  Service 
(Copeland  Report) 

USDA  Forest  Service 
(Reappraisal) 

USDA  Forest  Service  (Timber 
Resources  Review) 

USDA  Forest  Service 
(Timber  Trends) 

Gedney  et  al.  (PNW  Economic 
Base  Study) 

USDA  Forest  Service  (Douglas- 
Fir  Supply  Study) 

Vaux 


U.S.  President's  Advisory 
Panel  on  Timber  and  the 
Environment 

USDA  Forest  Service 
(Outlook) 


Gedney  et  al. 
Clawson  and  Hyde 


USDA  Forest  Service 
(Assessment) 


USDA  Forest  Service  (Timber 
Harvest  Scheduling  Issues 
Study) 

Rahm 


USDA  Forest  Service  (Analysis 
of  the  Timber  Situation) 

Brooks 


1920        Growth-drain  ratio  analysis 


1933        Growth-drain  ratio  analysis 


1948        Growth-drain  ratio  analysis 


1958        Growth-drain  ratio  analysis 


1965        Growth-drain  ratio  analysis 


1966        Stand  table  projection 


1969  Investment  analysis  and  even- 

flow  modeling 

1970  Investment  analysis  and  even- 

flow  modeling 

1973  TRAS  stand  table  projection 
and  supply/demand  "gap" 
analysis 

1973  TRAS  stand  table  projection 
and  supply/demand  "gap" 
analysis 

1975  TRAS  stand  table  projection 

1976  Yield-table  projection  with 

investment  analysis 

1976a  TRAS  stand  table  projection 
and  supply/demand  "gap" 
analysis 

1976b       Investment  analysis  and  even- 
flow  modeling,  with  price 
projections 

1981  TREES  inventory  projection, 

with  economic  simulation  of 
timber  harvests  utilizing  a 
downward  sloping,  regional 
demand 

1982  TAMM  economic  equilibrium 

modeling 

1983b       Yield  table  projection,  with 
investment  functions  and 
linkage  to  the  TAMM  model 


U.S. 
U.S. 
U.S. 
U.S. 
U.S. 


Pacific  Northwest,  Inland 
Empire 

National  Forests  in  Douglas- 
fir  region 

Douglas-fir  region 


U.S. 
U.S. 

Pacific  Coast  States 
Coastal  Pacific  Northwest 

U.S. 

National  Forests 
Washington,  Oregon,  Idaho 

U.S. 

Private  lands  in  the  South 


42 


Appendix  B.— Selected  recent  studies  of  growth  and  yield  modeling 
identified  by  principal  species  and  region 


Author  Year  Principal  species  Geographic  area 


Whole  stand  models 


Adams  and  Ek 

Alexander  et  al. 

Alig  et  al. 

Alimi  and  Barrett 

1974, 
1975 
1975 
1982 
1977 

Northern  hardwoods 

Engelmann  spruce 

U.S.  species 

Mixed  conifer  hardwood 

Lake  States 

Colorado,  Wyoming 
U.S. 

Northeast 

Barnes 

1955 

Slash  pine 

Florida 

Beck  and 
Oella-Bianca 

1970, 
1972 

Yellow  poplar 

Southern  Appalachians 

Bennett  et  al. 

1959 

Slash  pine 

Southeast 

Bennett 

1970a, 
1970b 

Slash  pine 

Southeast 

Bennett  and  Clutter 

1968 

Slash  pine 

Southeast 

Brodie  and  Rose 

1975 

Jack  pine 

Wisconsin 

Bruce  et  al. 

1977 

Douglas-fir 

Pacific  Northwest 

Buckman 

1962 

Red  pine 

Lake  States 

Buongiorno  and  Michie 

1980 

North  Central  hardwoods 

North  Central  States 

Burkhart  et  al. 

1972a 

Loblolly  pine 

Southeast 

Burkhart  et  al. 

1972b 

Loblolly  pine 

Southeast 

Burkhart  and  Strub 

1974 

Loblolly  pine 

Southeast 

Cao  et  al. 

1982 

Loblolly  pine 

Southeast 

Chambers  and  Wilson 

1  !ul  1  lUvl  \J    HI  IVJ    WWII  O \-/  1  1 

1972 

Douglas-fir 

Pacific  Northwest 

Chambers  and  Wilson 

1978 

Western  hemlock 

Pacific  Northwest 

Clutter 

1963 

Loblolly  pine 

Southeast 

Clutter  and  Allison 

1974 

Monterrey  pine 

New  Zealand 

Coile  and  Schumacher 

1964 

Loblolly  and  and  slash  pine 

South 

Curtis 

1967 

Coastal  Onunla^-f  ir 

Coastal  Pacific  Northwpst 

\y  VCiO  lul     1    tAVrf  Mils    1  1VI  U  lilCul 

r^iirtjc  pt  al 

VUI  llO  CI  CM. 

1981 

Coastal  noiinla^-fir 

KJ CI 0 1  CI  1    LfUU^jluO  III 

Coastal  Pacific  Northwest 

nahm^ 
Lycxi  II  i  io 

1Q64 
1975 

I  nrinpnnlp  ninp 

Orpnnn 

Dale 

1972 

Mixed  oaks 

Iowa,  Missouri,  Kentucky, 
Ohio 

Dell  et  al. 

1979 

Slash  pine 

South 

Dieraul  and  Mader 

1965 

Loblolly  pine 

Southeast 

Edminster 

1978 

Ponderosa  pine,  lodgepole 
pine,  spruce-fir 

Rocky  Mountains 

Ek 

1974 

Northern  hardwoods 

Lake  States 

Ek  and  Brodie 

1975 

Aspen 

Lake  States 
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Appendix  B.— Selected  recent  studies  of  growth  and  yield  modeling 
identified  by  principal  species  and  region  (continued) 


Author 

Year 

Principal  species 

Geographic  area 

Farr 

1967 

White  spruce 

Interior  of  Alaska 

Farrar 

1979 

Longleaf  pine 

South 

Fedducia  et  al. 

1979 

Loblolly  pine 

South 

Gardner  et  al. 

1982 

Mixed  hardwoods 

South 

Goebel  and  Shipman 

1969 

Loblolly  pine 

South  Carolina 

UUCUCI  CI  al. 

1974 

I  nhlollv  ninp 

South  Carolina 

OIUouMUW  dllU  LVdllo 

OlOOl  1   [Jll  IO 

Cni  1 1 hp^s  i 

OVJU  1 1  ICOo  l 

Hann 

1980 

Pnndprn^a  ninp 

Southwest 

Hoyer 

1975 

rjmmla^-fir 

L/UUUIQw  III 

Parifir  Northwp^t 
(Washington) 

1  nhlnllv  ninp 

Texas 

Lenhart  and  Clutter 

1971 

Loblolly  pine 

Georgia 

Lohrey  and  Bailey 

1977 

Longleaf  pine 

Louisiana  and  Texas 

Lynch 

1958 

Ponderosa  pine 

Northeast  Washington, 
Idaho  West  Montana 

Lynch 

1  SOL 

MiypH  ^npnip^*  a^npn-\A/hitp 

IVIIACU   OjJCOICo.  dO)JCI  1  WIIIIC 

pine 

MacKinney  et  al. 

1937 

Loblolly  pine 

Southeast 

McArdle  et  al. 

1961 

Coastal  Douglas-fir 

Western  Oregon,  Washington 

Meyer 

1938 

Ponderosa  pine 

West 

Moser  et  al. 

1979 

Mixed  hardwoods 

Lake  States 

Murphy 

1982 

Shortleaf  pine 

South  (West  Gulf) 

Murphy  and  Sternitzke 

1979 

Loblolly  pine 

South  (West  Gulf) 

Myers 

1971 

Ponderosa  pine 

Rocky  Mountains 

Myers  et  al. 

1976 

Ponderosa  pine 

Southwest 

Nelson  et  al. 

1961a 

Virginia  pine 

Southeast 

Nelson  et  al. 

1961b 

Loblolly  pine 

Southeast 

Payandeh 

1973 

Northern  species 

Southeastern  Canada 

Schlaegel 

1971 

Aspen 

North  Central  States 

Schnur 

1937 

Oak-hickory 

East 

Schreuder  et  al. 

1979 

Slash  pine 

South 

Schumacher  and  Coile 

1960 

Southern  pines 

Southeast 

Smalley  and  Bailey 

1974a 

Loblolly  pine 

Southeast 

Smalley  and  Bailey 

1974b 

Shortleaf  pine 

Southeast 

Smith  et  al. 

1975 

Mixed  hardwoods 

South 
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Appendix  B.— Selected  recent  studies  of  growth  and  yield  modeling 
identified  by  principal  species  and  region  (continued) 


Author 

Year 

Principal  species 

Geographic  area 

Solomon 

1977a, 
1977b 

Northern  hardwoods 

New  England 

Staebler 

1955 

Coastal  Douglas-fir 

Western  Oregon,  Washington 

Sullivan  and  Clutter 

1972 

Loblolly  pine 

South 

USDA  Forest  Service 

1976c 

Southern  pines 

South 

Warrack 

1959 

Douglas-fir 

British  Columbia 

Wiley  and  Murray 

1974 

Coastal  Douglas-fir 

Pacific  Northwest 

Worthington  et  al. 

1960 

Red  alder 

Pacific  Northwest 

Individual-tree  distance-independent  models 

Botkin  et  al. 

1972 

Northeast  mixed  forest 

Northeast 

Ek  et  al. 

1980 

Pure  or  mixed  species  stands 

Lake  States 

Goulding 

1972 

Douglas-fir 

British  Columbia 

Lemmon  and  Schumacher 

1962, 
1963 

Ponderosa  pine 

Rocky  Mountain  Region 

Martin 

1978 

Red  pine 

Wisconsin 

Oliver  and  Powers 

1978 

Ponderosa  pine 

California 

Stage 

1973, 
1979 

Ponderosa  pine,  Douglas-fir, 
lodgepole  pine,  mixed 
western  species 

Northern  Rocky  Mountains, 
Pacific  Northwest 

USDA  Forest  Service 

1979 

Lake  States  species 

Lake  States 

Individual-tree  distance-dependent  models 

Arney 

Daniels  and  Burkhart 
Daniels  et  al. 

1974 
1975 
1979 

Douglas-fir 
Loblolly  pine 
Loblolly  pine 

Pacific  Northwest 

Southeast 

South 

Dress 

1970 

Even-aged  pure 
species  stands 

Not  geographically 
limited 

Ek  and  Monserud 

1974 

Northern  hardwoods 

Lake  States 

Goldsmith 

1976 

Eastern  white  pine 

Northeast 

Hatch 

1971 

Red  pine 

Lake  States 

Hegyi 

1974 

Jack  pine 

Ontario 

Lee 

1967 

Lodgepole  pine 

British  Columbia,  Alberta 

Lin 

1974 

Douglas-fir, 
western  hemlock 

Pacific  Northwest 

Mitchell 

1969 

White  spruce 

New  England 

Mitchell 

1975 

Douglas-fir 

Pacific  Northwest 

Newnham 

1964 

Douglas-fir 

British  Columbia 

Newnham  and  Smith 

1964 

Douglas-fir, 
lodgepole  pine 

British  Columbia 
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Appendix  C— Selected  recent  timber  investment  and  cost  studies  by  region 
of  the  United  States  and  principal  species 


Author 

Year 

Principal  species 

Geographic  area 

Timber  investment  studies 

Adams  and  Ek 

1974 

Northern  hardwoods 

Lake  States 

Adams  and  Ek 

1975 

Northern  hardwoods 

Lake  States 

Amidon  and  Akin 

1968 

Loblolly  pine 

South 

Anderson 

1968 

Southern  pine 

Georgia 

Anderson  and 
Guttenberg 

1971 

Oak-pine  conversion  to 
loblolly  and  slash  pine 

South 

Beuter  and  Handy 

1974 

Douglas-fir 

Western  Oregon 

Broderick  et  al. 

1982 

Loblolly  pine 

Virginia 

Brodie  et  al. 

1978 

Douglas-fir 

Pacific  Northwest 

Brodie  and  Kao 

1979 

Douglas-fir 

Pacific  Northwest 

Callahan  and  Smith 

1974 

Black  walnut 

Midwest 

Chappelle  and  Nelson 

1964 

Loblolly  pine 

South 

Clawson  and  Hyde 

1976 

Douglas-fir 

Coastal  Pacific  Northwest 

Ek  and  Brodie 

1975 

Aspen 

Lake  States 

Fight  and  Gedney 

1973 

Douglas-fir 

Pacific  Northwest 

Flick  et  al. 

1980 

Loblolly  pine 

South  (highlands) 

Flora 

1966 

Ponderosa  pine 

Pacific  Northwest 

Forest  Industries  Council 

1980 

U.S.  species 

U.S. 

Gansner  and  Herrick 

1973 

Upland  oak 

Ohio 

Gedney  et  al. 

1975 

Pacific  Northwest  softwoods 

Pacific  Coast  States 

Hardie 

1977 

Loblolly  pine 

Mid-Atlantic 

Herrick  and  Morse 

1968 

Appalachian  forests 

Virginia,  W.  Virginia 

Hyde 

1980 

Douglas-fir 

Pacific  Northwest 

Jackson  and  McQuillan 

1979 

Northern  Rocky 
Mountain  species 

Montana 

Knight  and  McClure 

1974 

Southern  pine 

Southeast 

Koss  and  Scott 

1978 

Douglas-fir 

Pacific  Northwest 

Leak 

1980 

Northern  hardwoods 

New  England 

Lewis  and  Chappelle 

1964 

Southern  species 

Virginia 

Lundgren 

1966 

Red  pine 

Lake  States 

Manthy 

1970 

Softwoods/hardwoods 

Pennsylvania 

Marty 

1973 

Softwoods 

U.S. 

McCauley  and  Marquis 

1972 

Northern  hardwoods 

North 

Mills 

1976 

U.S.  species 

U.S. 
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Appendix  C— Selected  recent  timber  investment  and  cost  studies  by  region 
of  the  United  States  and  principal  species  (continued) 


Author  Year        Principal  species  Geographic  area 


Timber  investment  studies— continued 


Randall 

1977 

Douglas-fir 

Pacific  Northwf 

Row 

1973 

Southern  pine 

South 

Row 

1978 

U.S.  species 

U.S. 

Utz  and  Sims 

1981 

Upland  oak 

East 

Vaux 

1954 

Sugar  pine 

California 

Vaux 

1973 

California  conifers 

California 

Webster 

1960 

White  pine,  Norway  spruce 

Pennsylvania 

Cost  studies 

C^rtrtle  i  n 
VAJI  imi  1 

1Q71 

INUI  11  1 

1QRD 

u.o. 

Hilliker  e>t  al 
niltmct  cl  dl. 

1  sLca  Qtotac  cnci.oi.cic 
Lctrvc  Olalco  opcuico 

Ldl\c  oldlco 

Mill<;  ft  al 

IVIIIIO  CI  dl* 

1Qft? 

fviudix  diiu  r\uL»t;ici 

C /-» i  (thorn  cnfirioc 
wUUlllclll  optJUIcc 

Cr,i  it  h 

MOaR  el  al. 

1Q77 

oouinern  species 

Qai  ilh 

oouin 

Moak  et  al. 

1980 

Southern  species 

South 

Moak 

1982 

Southern  species 

South 

Somberg  et  al. 

1963 

Southern  species 

South 

Sunda  and  Lowry 

1975 

Loblolly  pine 

South 

Weaver  and  Osterhaus 

1976 

Loblolly  pine 

South 

Wikstrom  and  Alley 

1967 

Lodgepole  pine 

Montana,  Idaho 

Yoho 

1961 

Southern  species 

South 

Yoho  et  al. 

1969 

Southern  species 

South 
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Appendix  D.— Studies  of  retention  and  condition  of  timber  investments 
on  private  lands 


Author 


Year 


Program 


Geographic 
area 


Practice 
type 


Performance 
measure 


Alig  et  al. 


1980 


Kingsley  and  Mayer  1972 


Kurtz  et  al. 


Mills  and  Cain 


Nodine  and 
Marsinko 

Risbrudt  and 
Ellefson 

Stone 

Tennessee  Valley 
Authority  (TV A) 

Williston 


Williston  and  Dell 


1980 

1978 

1979 

1983 

1970 
1962 

1972 
1974 


Soil  Bank 

Mixture  of  cost 
share  and  private 

Agricultural 
Conservation 
Program  (ACP) 

Forestry  Incentives 
Program  (FIP) 

Soil  Bank 


Forestry  Incentives 
Program  (FIP) 


Soil  Bank,  ACP, 
CCC,  private 

Yazoo-Little 
Tallahatchie  Flood 
Prevention  project 

Civilian  Conservation 
Corps  (CCC) 


South 
North 
East 

U.S. 

South  Carolina 
U.S. 

Michigan 
Tennessee  Valley 


Pine  plantations 
Conifer  plantations 
Conifer  plantations 


Hardwood  and  conifer 
practices 

Pine  plantations 


Hardwood  and  conifer 
practices 


Pine  and  hardwood 
plantations 


Northern  Mississippi        Pine  plantations 


Northern  Mississippi       Pine  plantations 


Retention  and 
condition 

Condition 


Retention  and 
condition 


Timber  yield  and 
financial  return 

Retention  and 
condition 

Timber  yield  and 
financial  return 


Retention  and 
condition 

Retention  and 
condition 


Retention  and 
condition 
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GLOSSARY 


DIAMETER  CLASS:  A  classification  of  trees  based  on 
tree  diameter  (including  bark)  measured  at  breast 
height  (4-5  feet  above  the  ground).  D.b.h.  is  the  com- 
mon abbreviation  for  diameter  at  breast  height,  and 
2-inch  intervals  or  diameter  classes  are  commonly 
used. 

DIRECT  GROWTH  AND  YIELD  MODEL:  Models  in 
which  the  underlying  sampling  information  is  used 
directly  as  a  basis  for  projecting  forest  development 
for  a  particular  area  or  stand. 

EMPIRICAL  YIELD  TABLES:  Tabular  presentation  of 
the  yield  of  stands  or  trees  for  average  stand  condi- 
tions under  existing  management  practices. 

ENDOGENOUS  VARIABLES:  Variables  whose  values 
are  simultaneously  determined  by  the  model  and 
which  the  model  is  designed  to  explain. 

ENGINEERING  FUNCTION:  A  function  (or  model) 
based  on  parameters  estimated  according  to  tech- 
nical or  engineering  efficiency  considerations. 

EXOGENOUS  VARIABLES:  Variables  originating  from 
external  causes,  whose  values  are  determined  out- 
side the  model  but  influence  the  model. 

HARVEST  FLOWS:  Annual  or  periodic  estimates  of 
timber  that  will  be  cut  over  time  on  a  particular 
aggregate. 

HARVEST  SCHEDULING:  A  method  for  estimating  a 
planned  sequence  of  cutting  and  reforestation  activi- 
ties that  are  scheduled  according  to  a  prescribed 
norm,  including  the  acres  and  volume  harvested 
through  each  activity. 

INDIRECT  GROWTH  AND  YIELD  MODEL:  Models 
that  are  applied  to  subject  stands  or  areas  larger  in 
nature  than  those  from  which  the  underlying  sample 
growth  and  yield  data  was  obtained,  with  two  dif- 
ferent sets  of  data  used  for  model  construction  and 
application. 

INGROWTH:  Number  of  trees  that  grow  into  a  partic- 
ular timber-size  class  during  a  specified  interval  of 
time. 

INVENTORY:  Quantity  of  stumpage  existing  at  a  certain 
point  in  time  for  a  specified  geographical  area  (this 
term  is  also  often  used  to  refer  to  the  activity  of  col- 
lecting data  on  the  stock  of  timber). 

LONG  RUN  or  LONG  TERM:  A  period  of  time  in  which 
all  factors  involved  in  the  production  process  are  con- 
sidered variable. 

MARGINAL  ANALYSIS:  Economic  optimization  princi- 
ple that  contribution  to  total  revenue  of  an  additional 
unit  produced  (i.e.,  marginal  revenue)  should  equal 
the  contribution  to  total  cost  of  producing  an  addi- 
tional unit  (i.e.,  marginal  cost). 

NORMAL  YIELD  TABLES:  Tabular  presentation  of  the 
yield  of  "fully  stocked,"  undisturbed  natural  stands 
by  site-index  and  stand-age  categories. 


NORMATIVE  MODELING  APPROACH:  Analytical 
framework  designed  to  show  what  ought  to  be,  given 
certain  conditions  and  assumptions  (i.e.,  prescrip- 
tive). 

POSITIVE  MODELING  APPROACH:  Analytical  framed 
work  designed  to  describe  things  as  they  do  exist  or 
are  likely  to  exist  based  on  empirical  or  historical  evi- 
dence (i.e.,  descriptive  or  predictive). 

PRECOMMERCIAL  THINNING:  Removal  of  some  un- 
merchantable trees  from  a  forest  stand  to  promote  en- 
hanced growth  of  the  remaining  trees. 

PRODUCTION  FUNCTION:  Technical  relationship  be- 
tween inputs  and  outputs  in  a  production  process;  it  is 
often  expressed  as  an  equation,  as  when  the  quantity 
of  stumpage  (output)  is  defined  as  a  function  of  the 
relevant  inputs  (e.g.,  land). 

RESOURCES  PLANNING  ACT  (RPA):  The  Forest  and 
Rangeland  Resources  Planning  Act  of  1974,  which  re- 
quires integrated  planning  between  levels  of  state 
and  federal  agencies,  including  assessments  at  each 
level  that  culminate  in  a  national  assessment  of  the 
renewable  resources. 

SAWTIMBER:  Live  trees  that  are  above  a  certain  d.b.h. 
and  contain  merchantable  lengths;  for  example,  soft- 
wood sawtimber  trees  in  the  Pacific  Northwest  must 
be  11.0  inches  in  d.b.h.  or  larger  and  contain  at  least 
a  12-foot  sawlog  or  two  noncontiguous  8-foot  sawlogs. 

SENSITIVITY  ANALYSIS:  Procedure  by  which  values 
for  coefficients  or  variables  are  changed  and  the  im- 
pact on  a  target  value  or  solution  observed.  This 
technique  is  often  used,  for  example,  to  see  if  an  op- 
timal solution  is  highly  sensitive  to  the  values  used  for 
variables  that  are  not  known  with  certainty. 

SHORT  RUN:  Period  of  time  in  which  one  or  more 
factors  involved  in  the  production  process  is  con- 
sidered fixed. 

STOCK:  Entire  amount  of  standing  timber  existing  at  a 
particular  point  in  time;  it  is  also  termed  physical  sup- 
ply or  inventory. 

STOCKING:  Degree  of  occupancy  of  land  by  trees, 
measured  by  basal  area  and/or  number  of  trees  by 
size  or  age  and  spacing,  compared  to  a  standard. 

TIMBER  DEMAND:  Schedule  of  the  quantities  of  stump- 
age that  consumers  are  willing  to  purchase  at  differ- 
ent prices  in  a  market  specified  as  to  time  and  place. 
Demand  for  stumpage  is  "derived"  from  the  demand 
for  other  wood  products  such  as  lumber  and  plywood. 

TIMBER  SUPPLY:  Schedule  of  the  quantities  of  stump- 
age that  producers  are  willing  to  offer  for  sale  at  dif- 
ferent prices  in  a  specific  market  area  in  a  given  time 
period. 

TIMBER  SUPPLY  ANALYSIS:  Examination  and  investi- 
gation of  a  timber  supply  complex,  its  elements,  and 
their  relations,  in  the  context  of  a  functional  relation- 
ship between  quantity  of  stumpage  that  would  be  sup- 
plied (or  produced)  and  price. 
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U.S.  Department  of  Agriculture 
Forest  Service 

Rocky  Mountain  Forest  and 
Range  Experiment  Station 


The  Rocky  Mountain  Station  is  one  of  eight 
regional  experiment  stations,  plus  the  Forest 
Products  Laboratory  and  the  Washington  Office 
Staff,  that  make  up  the  Forest  Service  research 
organization. 

RESEARCH  FOCUS 

Research  programs  at  the  Rocky  Mountain 
Station  are  coordinated  with  area  universities  and 
with  other  institutions.  Many  studies  are 
conducted  on  a  cooperative  basis  to  accelerate 
solutions  to  problems  involving  range,  water, 
wildlife  and  fish  habitat,  human  and  community 
development,  timber,  recreation,  protection,  and 
multiresource  evaluation. 


RESEARCH  LOCATIONS 

Research  Work  Units  of  the  Rocky  Mountain 
Station  are  operated  in  cooperation  with 
universities  in  the  following  cities: 


Albuquerque,  New  Mexico 
Flagstaff,  Arizona 
Fort  Collins,  Colorado" 
Laramie,  Wyoming 
Lincoln,  Nebraska 
Rapid  City,  South  Dakota 
Tempe,  Arizona 


•Station  Headquarters:  240  W.  Prospect  St.,  Fort  Collins,  CO  80526 


